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Abstract 



Probing the Frontiers in QCD 

William A. Horowitz 

With the energy scales opened up by Rhic and Lhc the age of high-pr physics 
is upon us. This has created new opportunities and novel mysteries, both of 
which will be explored in this thesis. The possibility now exists experimentally 
to exploit these high momentum particles to uniquely probe the unprecedented 
state of matter produced in heavy ion collisions. At the same time naive 
theoretical expectations have been dashed by data. 

The first puzzle we confront is that of the enormous intermediate-pr az- 
imuthal anisotropy, or v 2 , of jets observed at Rhic. Typical lines of reasoning 
lead to an anticorrelation between v 2 and the overall jet suppression, Raa] the 
larger the v 2 the smaller the Raa- By simultaneously plotting the two this re- 
lationship becomes manifest and it is clear that usual energy loss mechanisms 
cannot reproduce the observed pattern — while jets are suppressed, the v 2 is 
anomalously large compared to the quenching. We argue that the data can be 
reproduced by a focusing of the partonic jets caused by processes associated 
with a deconfined quark-gluon plasma. 

The second puzzle is the surprisingly similar suppression of light mesons 



and nonphotonic electrons, which precludes perturbative predictions predi- 
cated on gluon bremsstrahlung radiation as the dominant energy loss channel. 
Near qualitative agreement results from including collisional energy loss and 
integrating over the fluctuating jet pathlengths. 

Another conjecture for heavy quark energy loss comes via explicit con- 
struction using the AdS/CFT correspondence; the momentum loss of a hang- 
ing dragging string moving through the deconfined plasma leads to qualitative 
agreement with heavy quark decay data. We propose a robust test to experi- 
mentally differentiate these two competing ideas: the ratio of charm to bottom 
suppression rapidly approaches 1 for pQCD but is independent of momentum 
and well below 1 for AdS/CFT. 

Finally as a warmup problem to calculating the photon bremsstrahlung 
associated with jet energy loss we quantify improvements to the perturbative 
estimates of the Ter-Mikayelian effect. By not neglecting the interference from 
away-side jets we find agreement with our results and the classical limit, regu- 
late divergences at low momenta, and note the importance of terms neglected 
in previous in-medium radiative energy loss derivations. 
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Chapter 1 
Introduction 

1.1 Philosophy of Physics 

It's hard to overemphasize the importance of the rational-reductionist tradition 
begun when Thales claimed that "All things are made of water" [1]. However 
the ancient Greek philosophers did not believe in the usefulness of experimen- 
tation; their investigations into Nature came from reasoning alone. Of course 
this emphasis on logic over the physical caused certain difficulties, culminating 
most famously with Zeno and his arguments against the possibility of motion 
[2]- 

Modern philosophy of science places primacy in experiment: the usefulness 
of a scientific theory is measured by its falsifiability [3]. And while many in 
the string theory community have reverted to measuring progress based on 
aesthetics [4], I am a devout Popperite. 

This thesis is broadly organized as follows. The experimental and theo- 
retical advances leading to QCD as the unique theory of the strong nuclear 
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force are reviewed. The existence of a transition from normal nuclear matter 
to a novel state is motivated, and some of the traditional theoretical tools and 
signatures of such are described. This places the work of the first half of the 
thesis in context. Then the recent application of conjectured strongly-coupled 
methods of AdS/CFT to heavy ion physics is detailed, which the latter half 
of this thesis proposes to test experimentally. 

1.2 Experimental Measurements Leading 
to QCD 

Any history of science is necessarily revisionist. The path to our current un- 
derstanding was not straight and certainly not chronological; in fact, there 
were plenty of dead ends. This is of course not what we as physicists would 
like to think, and in our attempt to make sense of the past the work of many is 
elided while a lucky few are picked out as the representatives for the discovery 
of now-obvious results. As this is not a thesis on the history of physics I will, 
with some regret, continue in this tradition; additionally I will focus on the 
developments leading to QCD as the fundamental theory of the strong force. 
For a more detailed narrative of physics from the late 19th century and beyond 
see, e.g., [5-7] and references therein. 

Reaching all the way back to the Greeks again, Leucippus and his pupil 
Democritus were the first to postulate atomic theory: matter (as opposed to 
the void) is made up of indestructible, individual particles [8]. It took until 
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the early 19 century for Dalton, inspired by his own experiments and the 
experimentally derived laws of conservation of mass [9] and definite proportions 
[10] in chemical reactions, to propose the precursor to the modern, scientific 
atomic theory [11, 12]. 

About a hundred years later, atoms began falling apart. In 1897 J. J. 
Thomson definitively demonstrated that cathode rays are made up of neg- 
atively charged particles, electrons [13]. Soon afterward, Rutherford's [14] 
observation of large angle scattering from gold foil showed that the majority 
of atomic mass is found in a minuscule, positively charged nucleus. It became 
clear from the transmutation experiments begun by Rutherford that the hy- 
drogen nucleus was one of the fundamental building blocks of all other nuclei; 
as such he named it the proton [15]. Chadwick's discovery of the neutron, 
simultaneously explaining the charge-mass asymmetry and providing a means 
of keeping the positive charge within the nucleus, completed the discovery of 
the constituents of atoms [16]. 

Yukawa [17] developed a theory for the force binding nuclei together, posit- 
ing the existence of an as-yet unmeasured particle whose mass in natural units 
is of the order of the nuclear radius, 1 fm _1 . After a period of confusion, Powell 
conducted the definitive experiments [18, 19] that disentangled the pion from 
the muon. Then things got weird. In December of that same year a kaon was 
first seen in a cosmic ray cloud chamber photograph [20]. Strange particles 
proliferated: the i], 0, and u mesons were found as was the A baryon. 

During the time that the muon-pion puzzle was being sorted out, Stuck- 
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elberg proposed the conservation of baryon number to explain the stability of 
protons. Experimentally, strange particles are produced on short timescales 
but decay relatively slowly; Pais suggested [21] that their production and decay 
mechanisms were different. Gell-Mann [22, 23] and Nishijima [24] expanded on 
this and Stiickelberg's idea by positing a new conserved quantity: strangeness. 

Gell-Mann, and independently Ne'eman (see [25] and references therein), 
began to understand the proliferation of hadrons by organizing them in mul- 
tiplets that are representations of the group SU(3), famously predicting the 
existence of the Q~ baryon (discovered in 1964 [26]). Even before this particle 
was observed, several authors (see [27, 28] and references therein) improved 
upon the phenomenological Eightfold Way by boldly postulating the existence 
of subnuclear structure. Gell-Mann named these smaller, fundamental build- 
ing blocks of hadrons quarks. 

The quark model was initially enormously successful. By taking these con- 
stituents to have spin-1/2 and fractional charge 2/3 for the u and -1/3 for the 
d and s, the spin, charge, and strangeness of all known hadrons was under- 
stood. Moreover the multiplets and their mass-ordering were well described. 
Using current algebra techniques scaling laws were proposed [29, 30] and found 
[31, 32]. These two SLAC papers also found large-angle scattering from pro- 
tons, which, like Rutherford's earlier experiments with atoms [14], showed 
decisively that nucleons have substructure. 

There were two problems with the quark model at this time: quarks had 
not been observed and some baryons apparently violated the Pauli principle. 
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Although the terminology came later in [33] , Nambu and Han [34] proposed a 
solution to the latter problem: quarks come not only in different flavors but 
also different colors. Phenomenologically one could posit that nature requires 
color neutrality to explain the former. In the same paper Nambu and Han in- 
troduced gauge vector fields associated with the quark color charge. Feynman 
rederived one of the previously mentioned scaling relations after developing 
the parton theory of nucleons: at high q 2 nucleons look like they are made up 
of point-like objects [35-38]. To agree with data nucleons have to be made up 
of not just the three valence quarks, but also sea quarks [39] and gluons [40]. 
While no experiment has yet to directly observe a bare quark or gluon (nor do 
we expect them to) there is strong indirect evidence for quarks from two jet 
events [41] and gluons from three [42] and four [43, 44] jet events; see Fig. 1.1. 

To describe the strong force binding quarks in nuclei, then, one wants a 
quantum field theory that becomes weaker at short distances. The idea of a 
non-Abelian gauge theory was introduced by Yang and Mills [45] ; quantization 
was achieved by Faddeev and Popov [46] , and 't Hooft proved their renormaliz- 
ability [47] . After asymptotic freedom was first shown for non-Abelian theories 
[48, 49], it was proved that in four dimensions this property is unique to non- 
Abelian theories [50]. The infrared divergences of non-Abelian theories led to 
a natural explanation of confinement [51, 52], the inability to directly observe 
quarks or gluons. 

Experimental evidence for QCD as the theory of strong interactions is 
legion. At high momenta, for which perturbative methods are applicable, log 
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Figure 1.1: Momentum representation of two (a) two-jet events and (b) three-jet 
events. Figure adapted from [42]. 

violations of Bjorken scaling in deep inelastic scattering were predicted [52- 
54] and observed [55, 56]. Next-to-leading order (NLO) calculations reproduce 
the world data of prompt photon production [57]. Heavy quark jet production 
rates are calculable and agree with experiment (see [58] and references therein). 
For a review comparing QCD theory and data at e + e~ colliders see [59]. 

At low momenta lattice calculations [60] show additional evidence for QCD. 
qq pairs experience a linear potential, numerically suggesting a mechanism 
for confinement; see Fig. 1.2 (a). Hadron masses found from lattice QCD 
are converging to those seen in the lab. Fig. 1.2 (b) compares experimental 
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measurements of nonperturbative quantities computed on the lattice and finds 
agreement to data to within statistical and systematic errors of 3% or less. 




0.5 1 1.5 2 2.5 3 LQCD/Exp't, (n f = 3) 

(a) (b) 

Figure 1.2: (a) Heavy quark potential at T = from LQCD calculations in [61]. 
The lattice results are well approximated by the Cornell potential, V(r) = —a/r + r. 
Figure adapted from [62]. (b) Lattice QCD results divided by experimental data 
for a range of nonperturbative quantities; the values agree to within statistical and 
systematic errors of 3% or less. Figure adapted from [63]. 

To resolve the Pauli problem mentioned earlier one needs at least three 
colors; several measurements demonstrate that Nature in fact uses no more 
than three (for a review see [64] and [65]). The ratio of the cross sections, 
cr(e + e~ — > hadrons)/cr(e + e~ — » = N C ^2Q 2 ^ tests both the number 

of colors (N c ) and active flavors {Nf) as a function of center of mass energy, 
shows that at energies above the bottom mass but below the top that N c = 3 
and Nf = 5; see Fig. 1.3. The decay rate of the Z° to hadrons compared to 
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Figure 1.3: R = a(e + e~ — > hadrons) / o(e + e~ — > / u + plotted against electron 
energy in GeV. Above 5 GeV, R is well approximated by N C ^2,Q 2 with N c = 3. 
Figure adapted from [7]. 

e + e~ depends on N c ; again, experimentally N c = 3. Similarly, the inclusive 
semi-hadronic decay rate of the r lepton to its semi-leptonic one is a direct 
measure of N c with N c = 3. And the decay of neutral pions to two photons is a 
direct measure of the square of the number of colors; the theoretical prediction 

64^) ~P = 7 - 7 eV ' (L1) 

where /„- = 92.4 MeV is the pion decay constant controlling the ir~ — > fiu, 
is in remarkable agreement with data (7.7 ± 0.6 eV). Four jet measurements 
identified directly the triple gluon vertex and found that the gauge group is 
SU(3) instead of SO(3) or perhaps U(l)3 [43, 44]; see Fig. 1.4. Similarly the 
running of the coupling agrees well for N c = 3 [66]. 
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Figure 1.4: Experimental verification of SU(3) as the gauge group of QCD from 
four jet events measured at LEP. N c , Cf, and Tp are properties of the group struc- 
ture functions f and generators T defined by f abc f* abd = 5 cd N c , (T a T^ a ) ij = 6ijC F , 
and Tr{T a T^ b ) = 5 ab TF, where repeated indices are summed over. Groups with 
N c ^ 3 are excluded by decay ratios while S0(3) and U(l)s do not have the correct 
ratios of group theoretic factors. Figure adapted from [44]. 



1.3 QCD Phase Diagram 



The possibility of a new state of nuclear matter accessible through the collisions 
of heavy nuclei has a long history (see [67] and references therein). From 
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the low energy side, Lee and Wick showed that scalar field theories could 
support an 'abnormal' nuclear state with properties far different from those in 
the usual vacuum [68]. In more modern language there is a nonzero vacuum 
expectation value for a quark condensate which one expects to melt at higher 
temperatures, thus restoring chiral symmetry. Both the statistical model of 
Hagedorn [69, 70] and the hadronic mean field approach of Walecka [71-73] 
predict a phase transition; see Fig. 1.5. 

From the QCD side, and before the advent of asymptotic freedom, Itoh was 
the first to suggest the possibility of deconfined quark matter [74] , followed by 
Carruthers (as cited in [67]). Collins and Perry [75] were the first to recognize 
the importance of asymptotic freedom, leading at large energies to a 'quark 
soup.' Shuryak coined the name 'quark-gluon plasma' (QGP) [76] as noted in 
[77]. 

The characteristic length for the polarization of the medium caused by 
color charges is the inverse of the Debye mass, \xd ~ gT, which is related to 
the temperature of the plasma. LQCD calculations of the static qq potential 
as a function of distance and temperature show that for T > T c the potential 
is screened; see Fig. 1.6 (a). Lattice calculations show a sharp rise of entropy 
density as a function of temperature, with the density afterward given by an 
ideal gas with quark and gluon degrees of freedom to within 20%; recent results 
using almost physical quark masses [78] are shown in Fig. 1.6 (b). There is 
disagreement within the lattice community as to whether the chiral transition 
and the deconfinement transition occur simultaneously. Nevertheless all these 
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very different lines of reasoning point to something very interesting occuring 
in nuclear matter at high temperatures and densities. 




Figure 1.5: (a) Phase transition from confined baryonic/mesonic matter to 
quark/gluon matter from the mean field quantum hadrodynamics (QHD) approach 
of Walecka [71]; figure adapted from [79]. (b) A sharp rise in energy reminiscent 
of a phase transition is seen at T ~ 160 MeV in Hagedorn's statistical model [69]; 
figure adapted from [70]. 

Interestingly these disparate descriptions of nuclear matter T high energy 
density all obtain similar values for the transition temperature T c . Simple di- 
mensional analysis leads one to expect the temperature of such a QGP to be 
around Aqcd ~ 200 MeV, as this is a necessary momentum scale to resolve dis- 
tances of order the nucleon and at the same time is the scale of normal nuclear 
energy densities of 1 GeV/fm 3 . Walecka yields ~ 150 MeV [79]. Frautschi 
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Figure 1.6: (a) The qq potential from LQCD for heavy quarks shows evidence of 
plasma formation at high temperatures. For T > T c the potential becomes indepen- 
dent of distance for separations longer than a temperature-dependent characteristic 
length; figure adapted from [80-83]. (b) Lattice calculations show a sharp rise in 
entropy density as a function of temperature, indicating a significant change in the 
number of degrees of freedom at T c ~ 180 MeV; figure adapted from [78]. 



found that Hagedorn's bootstrap gives T c rj 160 MeV [70]. Current data 
for T c from the lattice are in qualitative agreement with the previous results, 
although they are inconsistent quantitatively. The Wuppertal group found 
[84] chiral restoration at 151 (3) (3) MeV and a crossover phase transition at 
176(3) (4) MeV while the BNL/Bielefeld group found [78] T c = 192(7) (4) with 
the chiral limit smaller by 3%. 

In order to explore the new physics of nuclear matter in extremis an am- 
bitious program of heavy ion collisions began with Bevalac and continued 
through Acs, Aps, Rhic, and soon will commence at Lhc. 
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Chapter 2 
Motivation 

Since there are so many theoretical indications for a phase transition in 
QCD at low baryon chemical potential and at T > T c ~ 160 MeV, and since 
that unexplored phase would be a truly novel and interesting state of matter in 
which the ordinarily confined quarks and gluons — and not protons, neutrons, 
pions, etc. — are the pertinent degrees of freedom, one naturally asks how one 
might observe the creation of such conditions. Due to the techniques used to 
calculate them, the experimentally measured quantities associated with heavy 
ion collisions naturally separate themselves into low momentum, or bulk, and 
high momentum jet observables. 

Although heavy ion experiments have a long pedigree this discussion will 
focus mainly on Rhic data (see the white papers from the four experimen- 
tal collaborations at Rhic for a review [77, 85-87]); this thesis focuses on jet 
observables, for which cross sections were simply too small at previous exper- 
iments, and the data from Rhic, as discussed below, is qualitatively different 
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from these previous experiments in exciting new ways. 

2.1 Bulk Observables 

The most basic bulk observable is the energy deposited by all particles. Sur- 
prisingly, this elementary quantity provides a qualitative estimate of the energy 
density created in heavy ion collisions [77]. More differential measurements 
have the potential to provide much more information on the medium. Sim- 
ply taking ratios of particle species gives insight into the thermal properties 
of their creation [88-95]. Single particle spectra and their distribution over 
the reaction plane, through the use of hydrodynamics, hold the hope of deter- 
mining the bulk evolution and its equation of state. Two particle correlations, 
often quoted as Hanbury Brown- Twiss (HBT) radii [96, 97], measure the phys- 
ical dimensions of the fireball at freezeout and provide a consistency check for 
hydrodynamics, a test which it has consistently failed (for a review see [98]). 

2.1.1 Bulk Evolution 

A thermalized medium with a small mean free path compared to the system 
size can be understood using (relativistic) hydrodynamics (see [99-102] for 
reviews). One then might hope to learn about the physics of a heavy ion 
collision by measuring bulk observables and comparing them to results from 
hydrodynamic evolution. Schematically, hydrodynamics takes a given set of 
initial conditions and evolves them according to the known conservation laws of 
the system and a set of externally specified equations, for instance the equation 
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of state for ideal hydrodynamics and the stress fields for second order viscous 
Israel-Stewart hydrodynamics; see Fig. 2.1. Strong evidence for a deconfined 
QGP could come from a robust result whereby hydrodynamics with a QGP 
equation of state (EOS) reproduces experimental data while hydro with a 
hadronic EOS does not. It is interesting to note that the original conception 
of the QGP was of a weakly interacting plasma of deconfined quarks and 
gluons. For hydrodynamics to be useful, though, thermalization after the 
violent nuclear collisions must be both rapid and sustained. This of course 
requires a strongly coupled quark gluon plasma, or sQGP, as named in [103, 
104]; here 'strong' refers not to the strong force, but rather to the ratio of 
potential to kinetic energy, T, being greater than 1. However, this is not a 
universally accepted naming convention. [86] takes these properties as assumed 
and refers only to a QGP; for a more detailed discussion see [77]. 

Since hydrodynamics predicts the evolution of the entirety of the bulk there 
are a number of observables that can be compared to data. The simplest is 
single the particle spectra and their angular distribution with respect to the 
reaction plane. It is useful to Fourier expand the detected distribution, 



where the normalization (1N(pt) and t>2, or azimuthal anisotropy, are the two 
most important moments for heavy ion collisions. Early results from RHIC 
and ideal hydrodynamics were quite promising: while the v% of particle yields 
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Coordinate space: Collective interaction Momentum space: 

initial asymmetry pressure final asymmetry 

Figure 2.1: In non-central heavy ion collisions there is initial spatial anisotropy 
that the resultant collective pressure gradients evolve into anisotropy in momentum 
space. Hydrodynamics aims to quantitatively model this process to gain information 
on the medium and its properties. Figure adapted from [105]. 

generated by hydro was too large when compared to previous heavy ion exper- 
iments, those at Rhic matched quite well [106]; see Fig. 2.2. Unfortunately 
the HBT radii were not so well reproduced; see Fig. 2.3. 

However the picture can never be this crystal clear. Hydrodynamics is a 
set of evolution equations; the initial conditions are input, and, not surpris- 
ingly, what comes out of hydro is highly dependent on what goes in; see Fig. 
2.4. There have been some suggestions for testing the initial state [114, 115]; 
without good theoretical or, better, experimental control over the initial con- 
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Figure 2.2: (a) V2 scaled to the initial spacial eccentricity, e, as a function of the 
charge particle density per unit transverse area, low-px flow at previous colliders 
fell far short of the hydrodnamic limit, which is reached for the Brst time at Rhic. 
Figure adapted from [107]. (b) Early azimuthal asymmetry, V2{pt), predictions from 
[106] compared to Rhic data [108] as adapted from [77]. 



ditions statements drawn from hydrodynamical modeling will be inconclusive 
at best. 

There are two additional major complications in interpreting hydrodynam- 
ics results. The first comes from evolution beyond a thermalized medium into 
a viscous gas of hadronic particles. As the system expands with time it neces- 
sarily becomes too dilute for equilibrium to be maintained and hydrodynamics 
can no longer be a valid description. Often taken to occur at the same time 
is the breakup of the quark gluon plasma into ordinary hadronic matter. But 
the transition from quark and gluon degrees of freedom to baryons and mesons 
is not well understood. Unlike high-pr fragmentation — discussed later — this 
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Figure 2.3: HBT parameters [109, 110] compared to various hydrodynamics cal- 
culations [111-113]; none describe the data well. Figure adapted from [77]. 

hadronization is in the low-p^ sector, and there is no reason to believe that 
experimental measurements in, say, p + p collisions can provide independent 
information on the process. Similarly the low momenta involved makes it non- 
perturbative; there is little theoretical insight, either, although recombination 
[117, 118] and quark coalescence [119, 120] are two proposed mechanisms for 
hadronization. More sophisticated 3 + ID hydrodynamics treatments that in- 
clude hadronic rescattering (a so-called hadronic 'afterburner') have had good 
success in reproducing both single particle spectra and azimuthal anisotropy 
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Figure 2.4: Recent calculations including hadronic rescattering [116] that show 
the influence of the initial conditions (IC) on the output. For diffuse, Glauber-like 
IC ideal hydrodynamics slightly underpredicts data whereas sharper CGC-like IC 
require viscous effects to follow the experimental trend. 

[116]. Unfortunately even these highly computationally intensive calculations 
that follow evolution all the way until free streaming still do not describe the 
observed HBT radii. 

The second open issue is the role of viscous corrections. Past 3 + ID rel- 
ativistic hydrodynamic treatments of Rhic were all ideal; the viscous terms 
were explicitly set to zero. There have been some attempts to quantify the 
magnitude of these corrections [121]. However the rigorous treatment of vis- 
cosity in relativistic hydrodynamics is still an open problem. The inclusion of 
only the first order correction terms leads to unstable solutions with acausal 
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modes [122-124]. The second order Israel- Stewart formalism was shown to 
be stable and causal [125-128], but the need to specify 10 independent stress 
field component initial conditions remains a formidable problem. Much work 
is underway on modeling, with 2 + ID viscous treatments recently published 
[129, 130]. And there is some hope that with the inclusion of these viscous 
terms hydrodynamics will finally reproduce the HBT correlation radii as well 
as the single particle distributions [121]. However fully 3 + ID relativistic 
viscous hydrodynamics models are still unavailable, and will be so for some 
time. 

Nonetheless a truly quantitative understanding of the size of viscous effects 
allowed by data in heavy ion collisions is crucial for progress in the field. 
Teaney [121] estimated that in Rhic collisions the viscosity to entropy density 
ratio, rj/s, often shortened to viscosity to entropy ratio, could be no more 
than about 0.1. Utilizing elastic channels with perturbative cross sections 
only, weakly coupled pQCD gives rj/s ~ 1 [131-134]. The most famous and 
promising result of AdS/CFT is the calculation of rj/s = l/4n in a strongly 
coupled plasma; for more on the AdS/CFT conjecture see Section 2.3. 

Currently the best bridge between 3 + ID ideal and viscous hydrodynamics 
comes from parton transport theory [135-137]. For parton cascades including 
only 2—7-2 processes, the ideal hydrodynamics limit is only approached when 
using very large elastic cross sections [138]. Nevertheless first comparisons 
between viscous hydrodynamics calculations and those from the cascade have 
been made [139]. Recent work by Xu and Greiner [137] has incorporated 
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2 <H- 3 processes; utilizing this promising new code they have found, using 
perturbative ideas alone, rj/s ~ .1 [140] and also hope to include both bulk and 
jet dynamics within a single theoretical framework [141]. There are plenty of 
difficulties, though, as including gluon splitting and merging makes numerical 
acausal artifacts more pronounced (their work has yet to be confirmed by 
other groups) and the gluon radiation from jets does not model coherence 
effects well. As noted previously it is hard to overemphasize the importance 
of initial conditions (IC) when considering hydrodynamics output. On the 
crucial issue of viscosity the two are necessarily coupled: consistency with 
data with more diffuse IC requires a more ideal fluid; sharper edges a more 
viscous one. Hopefully by independently experimentally testing the IC — see 
[115] for a proposed discerning observable — and by theoretically motivating the 
fluid viscosity we can arrive at a self-consistent picture of the bulk dynamics 
in heavy ion collisions. 



2.2 High-pr Observables 

The exciting first evidence of jets came from Isr in 1972 [142-146]. Orders 
of magnitude more high-p-^ pions were observed than were expected from a 
low-pr extrapolation; the production spectrum had turned over from exponen- 
tial to power law. Soon afterward, two jet events were explicitly seen at e + e~ 
colliders [41, 42]. While Bjorken was the first to suggest using jet suppression 
to learn about a QCD medium [147], the precision pQCD predictions of pro- 
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duction rates [148-151] held out the possibility for jet tomography: much like 
in medical applications such as a Pet scan, a careful measurement of the jet 
quenching pattern would reveal information on the medium through which the 
probe traveled. As high transverse momentum parton jets are produced early 
in the collision (by Heisenberg's Uncertainty Principle) and preferentially deep 
within the fireball — more on this later — they are potentially excellent probes 
of the medium. 

Naturally then before investing the tremendous time and resources nec- 
essary to investigate high-p^ particles theoretically and experimentally, one 
should ruminate on the epistemology of hard probes. There were claims in 
[77] that jet suppression can merely give information on the density of scat- 
tering centers. This is an important measurement, as noted previously a high 
density is a prerequisite for QGP formation, but jets have the potential to in 
fact reveal much more than simply a mean density. The possibility of testing, 
e.g. deconfinement, depends heavily on the energy loss model and the mapping 
made between its input parameters and the physical medium. This mapping 
is a critical, although often overlooked, component of any model attempting 
tomography. As an example, the WHDG model [152] (see Chapter 4) explic- 
itly assumes a connection between the medium density, its temperature, and 
its Debye mass: light quarks and gluons are taken to have masses proportional 
to no, and the energy loss results are surprisingly sensitive to its changes; the 
elastic energy loss is derived from classical considerations that assume plasma 
polarization. Additionally Chapter 3 argues from a phenomenological ap- 



Chapter 2: Motivation 



23 



proach to jet suppression that the large magnitude of the observed anisotropy 
at intermediate-pT is a result of deconfinement physics. 

2.2.1 Factorization in p + p and A + A 

Factorization provides the theoretical framework within which pQCD calcula- 
tions are made and states that reactions with large momentum transfer can 
be factorized separately into long distance and short distance pieces [153-160]. 
There are two crucial components to factorization: (1) the (presently) incal- 
culable nonperturbative low momentum contributions are universal; i.e. they 
are the same for any QCD process, and (2) asymptotic freedom guarantees 
that the high momentum contribution can be reliably found using pertubative 
methods in QCD. A factorization scale, /ip, is introduced separating the short 
and long distance contributions; it is of order the hard scale of the problem, 
but is not specified further from within the theory. Following [150] the cross 
section for producing a hadron of type h in a p + p collision is then 



where the fs are parton distribution functions (PDFs), D% is a fragmenta- 
tion function (FF), and da c ah is the perturbative partonic cross section. Using 
standard notation the partonic variables are in lower case and the hadronic 
ones are in upper case. In this way a PDF gives the probability of finding 



do' 



h 




(2.2) 
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parton a in a proton with momentum fraction between x a and x a + dx, where 
Xa = Pa/Pp arid something of an abuse of notation occurred by which p a stands 
for the momentum of parton a and P p stands for the momentum of the proton. 
The FF gives the probability of parton c hadronizing into h with fractional 
momentum between z c and z c + dz, where Ph = z c p c . All allowed combinations 
of a, b, and c are summed over where a + b — » c + X such that c may hadronize 
into h. We note that we have slightly extended the factorization theorem to 
include the fragmentation process, which introduces the new scale fi' F . Addi- 
tionally a renormalization scale associated with the running coupling a s is 
included in the perturbative cross section. 

Parameterizations of PDFs may be found from, e.g., the CTEQ collabora- 
tion [161], GRV [162], MRST [163], or Alekhin [164]. Hadronic fragmentation 
functions have similarly been parameterized by, e.g., DSS [165], AKK [166], 
HKNS [167], and KKP [168]; [169] provides a photon FF. We note that the 
posted code for KKP has a bug in its leading order (LO) kaon output that 
has still not been corrected; ported Mathematica codes for the DSS and KKP 
fragmentation functions, as well as a corrected Fortran code for KKP can be 
found online [170]. Although PDFs and FFs at a specific scale must be found 
from experiment, their form is constrained by QCD and general sum rule con- 
siderations and their evolution is governed by pQCD and DGLAP dynamics 
[171-173]. 

Derivations of perturbative cross sections vary depending on the specific 
reaction considered. For those involving gluons and light guarks NLO calcu- 
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lations exist [150, 151]; [174, 175] give derivations to NLO for prompt photon 
production; fixed order next to leading log (FONLL) production spectra of 
heavy quarks can be found in [176, 177]. 

In heavy ion collisions Eq. (2.2) is altered in three significant ways. First, 
the geometry of collisions of large nuclei of mass number A is very different 
from p+p; for a review of heavy ion collision geometry, see Appendix A. Large 
nuclei suffer a number of hard p + p-like reactions as they overlap, which is a 
function of the nuclei involved and the impact parameter b. Second, the PDFs 
used must be exchanged from nucleon to nuclear. To lowest order these new 
distributions should account for the isospin asymmetry of nuclei. In discussing 
these in detail beyond this simple correction it is useful to define the quantity 



F? /A (x,Q) 
F 2 a/p (x,Q) 



(2-3) 



where we have used the structure functions defined as the sum over all fla- 
vors of the PDF's weighted by x and the quark charge squared, F 2 (x,Q) = 
J2 f xQjff(x,Q). In the shadowing region x < .1, < 1 [178-181]. .1 < 
x < .3 is the anti-shadowing region in which > 1 [182]. In the EMC 
•3 ^ x < .7 region R£ < 1 ([183] and references therein), and in the Fermi 
motion x > .7 region > 1 [184-186]. [187-190] give some parameteriza- 
tions of nuclear structure functions. Finally, and most importantly for this 
thesis, as the high-pr parton travels through the fireball it interacts with col- 
ored scattering centers in the medium. This is encapsulated in a probability 
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of energy loss -P(e), where p T = {1 — e)p % T . Although it is called a probability 
of energy loss there are processes — e.g. 3 — > 2 detailed balance and rare elastic 
events — in which momentum is gained and e < (e is certainly restricted from 
above by 1; more on this later). Given this Eq. (2.2) becomes 



where P(e) in general depends on the position of hard parton production x , 
its direction of propagation <fi, and the properties {p}i of the medium; we will 
often refer to P(e; xq, <f), {p}i) in shorthand as P{e). See Fig. 2.5 for a schematic 
picture of Eq. (2.4). 

2.2.2 Raa and Final State Suppression 

A common means of reexpressing the spectrum of hadrons produced in a nu- 
clear collision A + B — ?■ /i + X is to divide out by the trivially binary scaled 
proton-proton spectrum: 



For symmetric collisions (such as Au + Au or Pb+Pb) we will denote the above 
quantity as Raa- Often Raa is given as a function of transverse momenta pr 




a,b,c 





(2.4) 




(2.5) 
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Figure 2.5: An illustration of A + B — > h + X in a heavy ion collision. We assume 
factorization still holds such that the main alteration to p+p collisions, Eq. (2.2), is 
the addition of a theoretically calculable -P(e) encapsulating the in-medium energy 
loss. 

and the angle with respect to the reaction plane 0. Older, lower statistics data 
is sometimes reported as RAA(N pBX t), although [191] demonstrated its limited 
utility. 

Raa is a nice quantity to work with as it displays the effect of nuclear col- 
lisions at a glance: Raa = 1 means no modification from trivially scaled p + p 
collisions, Raa > 1 means enhancement, and Raa < 1 means suppression. 
For this reason, Raa is known as the nuclear modification factor. When the 
partonic production spectra are power laws of approximately constant power, 
which is the case at Rhic and Lhc, Raa is approximately given by an espe- 
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daily simple relation; see Appendix B. 
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Figure 2.6: Suppression ofn and charged hadron jets in central Au+Au collisions 
compared to binary scaled p + p reactions at \fs = 130 GeV at Rhic from 2001. 
Figure adapted from [192]. 



Early data from Rhic showed a stunning suppression of hadronic jets; see 
Fig. 2.6. The immediate suspicion was that this was due to final state energy 
loss. However the large quenching of jets could result from an alteration in 
other piece of Eq. (2.4): binary scaling, initial state, hard cross section, or 
fragmentation. Taking these in reverse order, uncertainty arguments make 
changes to the fragmentation function unlikely for high-p^ pions: in its rest 
frame Heisenberg gives that a ~ 140 MeV pion should take ~ 1.4 fm/c to form; 
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a measured 5 GeV pion has a boost factor of 7 ~ 36 which means that in the lab 
frame it is formed ~ 50 fm away from where it was originally produced. This is 
well outside the medium and at a time much later than any QGP is reasonably 
expected to exist so it seems the use of vacuum fragmentation functions is 
quite safe for light mesons. We note that the same analysis for a 5 GeV 
detected proton yields a distance of just ~ 1 fm, most likely well within the 
medium. As noted previously, recombination [117, 118] and quark coalescence 
[119, 120] are two proposed models of medium-modified hadronization; there 
are even indications that heavy quark mesonic states, whose large mass also 
precludes a necessary vacuum fragmentation, might also survive in the QGP 
phase [193, 194]. Changes to the hard cross section are also unlikely as for 
large momentum transfers asymptotic freedom must make the coupling small. 

While it seemed unlikely that binary scaling did not hold in Rhic colli- 
sions many believed the initial state PDFs were significantly altered due to, 
for instance, the color glass condensate (CGC) [195-203, 203-216]. The ef- 
fects of large nuclei on all the aforementioned possibilities were tested by the 
deuteron-gold, d + Au, so-called 'control run.' By smashing these together at 
the same energy per nucleon as the A + A run changes to binary scaling, hard 
cross sections, and nuclear distributions were minimized while removing the 
medium, and thus final state energy loss or any modifications of the fragmen- 
tation functions. Fig. 2.7 (a) shows that R^a ~ 1 for high-pr hadrons (with 
a rather large Cronin enhancement at lower momenta [217]), thus falsifying 
the claims of CGC-type suppression and showing that the hard production in 
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heavy ion collisions is under theoretical control. 




(a) (b) 

Figure 2.7: (a) d + A results from STAR [218] show that high-p? suppression is 
not due to initial state effects; figure adapted from [86]. (b) R^ A (N part ) ~ 1 extends 
the results of (a) to show that jet suppression is not due to initial state effects in 
full A + A collisions. NLO prompt photon predictions [219-222] compare well to the 
PHENIX data [223]; figure adapted from [223]. 

The high-pT conclusions based on the d + Au control-run were indepen- 
dently verified and extended by the measurement of direct photons [223]. 
As electromagnetic probes are only weakly coupled (oe&m ~ 1/137) to the 
medium they do not suffer appreciable final state energy loss, even in A + A 
collisions. Fig. 2.7 (b) compares trivially scaled pQCD predictions for direct 
photons from p + p reactions [219-222] to those measured in Au + Au. The 
striking consistency shows that even in A + A collisions the initial state pro- 
duction is well understood; the variation in nuclear PDFs from p + p to d + A 
to A + A do not affect high-py observables much. 

Since the observed jet suppression cannot come from other sources, it must 
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be due to final state energy loss. Now IF energy loss is under theoretical con- 
trol and IF the combined theoretical and experimental system is not fragile, 
then jet tomography is possible and the observed suppression pattern can be 
inverted to gain knowledge about the medium. 

2.2.3 Elastic Energy Lost History 

Final state suppression calculations began with Bjorken's estimate of the 
energy lost by a high-p-r parton through elastic, 2 — > 2, processes [147]. 
He considered the simple t-channel diagram associated with a parton trav- 
eling through a thermalized, deconfined QGP of temperature T. In this case 
\M\ 2 ~ l/t 2 ~ l/(q 2 + M 2 ) 2 , where M ~ /xp ~ gT is the infrared cutoff given 
by the Debye mass. With some simplifying assumptions, weighed with the 
appropriate flux and kinematic factors this gives an austere analytic formula 
for the energy lost per unit distance, 



where Cr is the jet Casimir, 4/3 for a quark or 3 for a gluon, Nf is the number 
of active flavors in the medium, and (k) ~ 3T is the average momentum of 
the plasma particles. 

This estimate has been improved upon by a number of subsequent papers. 
Thoma and Gyulassy [224] calculated the linear response to a jet in classical 
Abelianized QCD, incorporating the hard thermal loop (HTL) work of Braaten 




(2.6) 
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and Pisarksi [225-227] (for a review of HTL, see [228]); the energy loss was 
found by deriving the work done on the source current by the induced fields 
in the dielectric medium. Braaten and Thoma [229, 230], expanding upon the 
work of Svetitsky [231], separately evaluated the high momentum part of the 
dynamics using the vacuum matrix element and the low momentum piece with 
linear response and connected the two at a scale q* which, to leading order, 
drops out of the problem [232]. 

These leading log results were compared to asymptotic radiative loss ap- 
proximations and found to be small [233, 234], not surprising as this is a well 
known result from E&M [235]. Further work on elastic mechanisms slowed 
until [236] showed that with realistic kinematic limits at Rhic energies ra- 
diative and elastic losses are of the same order. This and the experimental 
evidence of surprisingly strong heavy quark quenching [237-240] motivated 
work to incorporate both radiative and collisional processes in jet quenching 
models (see Chapter 4) and has generated a flurry of interest in improving 
elastic calculations. 

Some more recent developments in elastic loss include considerations of run- 
ning coupling, finite creation time effects, and the importance of the higher 
moments of the distribution. Peshier [241] saw a large energy loss enhancement 
when he allowed the coupling to run; [242] found a qualitatively similar result. 
Wicks [243] , however, showed that for Rhic and Lhc kinematics this is actu- 
ally a small effect. All previously mentioned calculations used asymptotic jets 
created in the infinite past. Peigne, Gossiaux, and Gousset [244] claimed that 
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finite formation time effects were large and persisted far beyond the expected 
single Debye length. Consideration of only the elastic pole contributions by 
Adil et al. [245] reproduced the intuitive result; subsequent work by Gossiaux 
et al. [246] confirmed this. Djordjevic tackled the problem from a quantum 
mechanical standpoint and came to the same conclusion [247]. Nevertheless, 
the effects of finite time and off mass-shell creation are still under investigation 
[248]. 

Svetitsky [231] was the first to include fluctuations about the mean in 
elastic energy loss; he employed the Fokker-Planck equation. Moore and 
Teaney [249] thoroughly investigated the relativistic Langevin and Fokker- 
Planck equations in heavy ion collisions. This work was applied to heavy quark 
energy loss with additional nonperturbative mesonic resonances by Rapp and 
van Hees [250] . As shown by Wicks, however, for the pathlengths and densities 
at Rhic and Lhc the number of 2 — > 2 collisions is of order a few, and the 
Gaussians resulting from these applications of the Central Limit Theorem are 
not a good approximation to the elastic energy loss distribution [243]. 

2.2.4 Background Radiation 

Rigorously deriving radiative energy loss is a tough business. Even in Abelian 
QED tremendous effort and theoretical contortions ('reinterpretation') were 
required to satisfactorily deal with the infrared and ultraviolet divergences 
(c.f., e.g., [251]). Undaunted, a number of nuclear theorists have devoted their 
lives to conquering nonabelian QCD radiative bremsstrahlung. The formalisms 
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created to tackle this problem can be roughly categorized into four groups: (1) 
BDMPS-Z-ASW (often simply BDMPS) [252-270], (2) DGLV (GLV) [271- 
279], (3) WWOGZ (Higher-Twist) [280-285], and (4) AMY [133, 134, 286, 
287]. 

Gunion and Bertsch [288] began the translation of QED into QCD by 
deriving the strong force field theory diagrams associated with the nuclear 
analog of incoherent Bethe-Heitler radiation, 

dN g = C 2 c a s qj 
dr]d 2 q ± n q 2 ± {k± - q ± ) 2 ' 

where q = (q°, q) = (q°, q\\, qj_), and C c are the color algebra matrix elements. 
Multiple coherent scattering over lengths shorter than the radiation forma- 
tion time leads to interference that the suppresses the emission of radiation; 
this is the well-known LPM effect in electromagnetism, named after Landau 
and Pomeranchuk [289] and Migdal [290]. Brodsky and Hoyer [291] began the 
work of including these effects in QCD calculations, and it was continued by 
Gyulassy and Wang [292]. This paper also introduced the notion of a static 
color scattering center with a Yukawa-like screened potential and noted the 
importance of the unique nonabelian extension of the LPM effect in QCD, 
by which the radiation reinteracts with the medium via the 3-gluon vertex. 
Baier, Dokshitzer, Mueller, Peigne, and Schiff [253-255] were the first to in- 
clude this effect in an energy loss calculation. Unlike the GW paper that 
examined the thin plasma limit, like Landau and Pomeranchuk they built up 
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the soft gluon radiation from single hard scatterings, BDMPS examined the 
multiple soft scattering limit, similar to Migdal and Moliere scattering [293- 
295]. Contemporaneous to the BDMPS papers, Zakharov developed his own 
formalism employing path integration [261-264], later shown to be equivalent 
to the BDMPS approach [256]. 

The thin plasma limit of GW was extended in the opacity expansion work 
of Gyulassy, Levai, and Vitev [271-274]. The reaction operator approach they 
developed allowed the derivation of a closed form solution of the resummed 
single inclusive gluon radiation distribution dN g /dxd 2 k to all orders in opacity, 
X — L/X. At a similar time, Wiedemann examined the dipole path integral, 
opacity expansion, and the Zakharov and BDMPS limits [266-268]. He and 
Salgado [269] numerically investigated the BDMPS and GLV results, pub- 
lishing a popular public code for calculating BDMPS energy loss 'quenching 
weights.' 

In QED the infrared divergences from single scattering bremsstrahlung 
exponentiate into a Poisson distribution of multiphoton fluctuations about 
the semiclassical expectation value [251]. Multigluon fluctuations were in- 
corporated in [275] by assuming this Poisson form: P{e\pr) = T n -^( e |Pr) , 
P\{Avt) = exp(-(N g (p T )))dN g /dx(x;p T ), and 



where the final momentum is expressed in terms of the initial momentum as 




(2.8) 
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p T = (1 — e)p T . The importance of multigluon correlation effects is currently 
unknown and is an important open theoretical problem. 

Also associated with the Poisson convolution is probability leakage, in 
which P(e > 1;Pt) unphysically has nonzero weight. For large regions of 
parameter space at Rhic and Lhc this leakage is quite large due to kinematic 
constraints neglected in the Poisson approximation. Two common approaches 
to dealing with this leakage, P ex cess = P{e)de, are to reweigh or truncate 
the distribution. When reweighing the excess is redistributed evenly by renor- 
malizing the probability: -P r ( e ) — -P(^)/-Pexcess- Truncation takes the excess as 
total jet absorption: -Pt(e) = P(e)9(l — e) + -P xccss<5(e — 1). Since the Salgado- 
Wiedemann (SW) quenching weights [269] use infinite jet energy, calculations 
for realistic jets with finite energy suffer the same problem. Dainese, Loizides, 
and Paic [296] took two methods for redistributing this excess probability and 
compared the resulting RaaS] these were taken as an error band, which turned 
out to be quite large. 

Neither of these proposed schemes is particularly satisfying. While P(e) 
with weight for e > 1 implies that the estimated number of emitted gluons 
is high, there is no reason to expect the shape of the true distribution to be 
the same but with a different normalization. And obviously the true shape is 
different from the truncation method. Relative branching is another Poisson 
approximation that is an improvement on truncation and renormalization. In 
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this case 



1 



/ 



dz 



Pn(z\ PT ) 



dN g e-z 



\Pt), 



(2.9) 



71 + 1 



1 - 



Z 



dx 1 — z 



where p n = (1 — z)p l and pf = (1 — e)/(l — z)p n = (1 — e)p l . In [243] rela- 
tive branching results closely resembled those from truncation; redistributing 
probability evenly gives excess weight to low e probabilities, which — due to 
the steeply falling spectrum — disproportionately affects the results. As can be 
seen from the above equations the Poisson convolution must have all interme- 
diate energy loss steps evaluated at the same pr, or else it will not be properly 
normalized. For large fractional energy loss, the apparent case at Rhic, this 
is not a good approximation; it is not clear how large an effect this has on 
physical observables. 

Around the same time Wang began the development of the 'Higher-Twist' 
formalism [280-284]. Similarly to GLV, the derivation builds up the energy 
loss from single hard scatterings but differs by making some alternative as- 
sumptions in their evaluation. Most important, the arbitrary potential in 
GLV — usually taken as Yukawa — is replaced by an arbitrary gluon distribu- 
tion function. This obscures the relation between jet suppression patterns and 
physical medium quantities such as density and temperature. For a review 
comparing the two results see [297]. 

Motivated by the 'dead cone' work of Dokshitzer and Kharzeev [298] — just 
as in QED, in QCD a massive charge also radiates less — and hoping for a 
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consistent theoretical description of gluon, light quark, and heavy flavor sup- 
pression these groups' work was extended. Zhang, Wang, and Wang [285] and 
Armesto, Salgado, and Wiedemann [270] included heavy quark mass effects 
in the WWGO-Z and BDMPS-Z-ASW formalisms. Djordjevic and Gyulassy 
included both the effects of a heavy quark jet and a gluon mass term in D-GLV 
[278]. 

Since then Arnold, Moore, and Yaffe [133, 134, 286] developed a thermal 
field theoretic derivation of energy loss in which dominant and subdominant 
contributions are carefully tracked; [299] found the result to be similar to 
BDMPS. A large advantage of this formalism is its simultaneous treatment of 
both gluon and photon bremsstrahlung, providing an added experimental con- 
sistency check. Unfortunately their use of asymptotic states neglects the large 
interference effects from the initial production radiation, making comparison 
to data problematic. 

Besides in-medium inelastic energy loss two other radiation effects have 
been studied: transition radiation and Ter-Mikayelian radiation reduction. 
Transition radiation occurs in E&M when a relativistic charged particle propa- 
gates through an inhomogeneous medium, in particular the boundary between 
two spaces with different electrical properties [300, 301]. In heavy ion colli- 
sions just such a boundary forms between the deconfined QGP medium and 
the vacuum; [302] quantified the extra energy loss caused in this transition and 
detailed its regulation of infrared divergences ordinarily absorbed into DGLAP 
evolution. 
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The Ter-Mikayelian (TM) effect [303, 304] is a direct result of radiative 
quanta gaining mass in a plasma. As in beta decay, production of high mo- 
mentum charged particles also has radiation associated with the process. For 
QCD this infrared divergent vacuum radiation is absorbed into fragmenta- 
tion functions, but in-medium the finite gluon mass regulates and suppresses 
this radiation. Djordjevic and Gyulassy [305] calculated the QCD analog of 
the TM effect for single quark pairs. Chapter 6 extends this work to include 
back-to-back jet production. We find the away-side jet, necessary for charge 
conservation, fills in the dead cone for heavy quarks, further suppresses the 
production radiation, and naturally regulates the momentum loss as the jet 
momentum approaches zero. 

2.2.5 Fragility 

Previously we argued that by comparing theoretical model results to data jets 
have the potential to provide information on not just the density of scattering 
centers but also of the nature of the matter in the medium. The process would 
be as follows: models qualitatively inconsistent with experimental observations 
are falsified; models in quantitative agreement with data set limits on medium 
properties by their input parameter range allowed by data. Clearly the more 
precise the experimental observable, the more tightly constrained the medium 
property. Similarly the more sensitive the theoretical results to changes in 
medium property the stronger the statement about the medium. Of great 
concern but poorly investigated is the influence of theoretical imprecision — 
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more on this later. 

[306] was the first to coin the term 'fragility,' which we will take to mean 
the limit on information about the medium that can be learned by inverting 
the experimental data with theoretical modeling. The paper raised the issues 
of surface emission, sensitivity, and fragility, and used these terms rather in- 
terchangeably. In this thesis each term has a separate meaning and will be 
addressed individually. 

As noted before the production of high-pr partons is biased toward the 
center of a heavy ion collision because the process scales like the number of 
binary collisions. Due to the large background of low-p^ particles full jet 
reconstruction algorithms are not yet available for Rhic, although both Star 
and Phenix are actively researching them [307-309]. Limited to measuring 
the leading hadron of a jet, experimental cuts combined with the in-medium 
energy loss and steeply-falling production spectra bias the production point 
of observed high-py particles toward the surface of the medium; see Fig. 2.8. 
Naively one expects that there will always be a corona of emission, or surface 
emission: near the very edge of the fireball the medium is so thin high-p^ 
partons escape with little or no energy loss. This naturally leads to the idea 
of a minimum nuclear modification value, -R™™, below which jets cannot be 
further suppressed. In fact this is not a correct conclusion for systems with 
realistic nuclear densities because at the outer edge of the medium no hard 
scatterings occur: the thinness of the nuclear overlap results in no binary 
collisions. In fact Section 5.2 explicitly shows Raa for large energy loss. 
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Even worse the phrase 'surface emission' leads to a Boolean mindset: Is there 
surface emission or not? This is not a useful question for a qualitative, let 
alone quantitative study of data; a better question would be 'How biased are 
the jets?' Surface emission as a term will not be used further in this thesis. 



Hard Cylinder Hard Sphere Woods-Saxon 




x(1m) 



Figure 2.8: Examples of surface bias in different nuclear geometries. Simplified 
nuclear geometries such as hard spheres, and especially hard cylinders, create greater 
surface bias. 

Theoretical sensitivity quantifies the change in a model-predicted observ- 
able given a change in the input parameter(s). Sensitivity is a slippery slope 
that may easily slide into the idea of an insensitive model. Just this type of 
Boolean thinking led to ironic conclusions from [306] who claimed that their 
calculation of Raa became insensitive to increases in q. But by quantifying 
the decrease in Raa as a function of q in a very similar model [310] found a 
power law relationship in the 'insensitive' results, Raa ~ q~ 1 ^ 2 . In this way a 
fixed fractional change in q leads to a fixed fractional decrease in Raa , which 
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is in agreement with the naive interpretation of a sensitive theory. The error 
in interpretation becomes clearer upon looking at plots of Raa vs. q] on a 
linear-linear scale Raa does appear to saturate (Fig. 2.9 (a)), but is clearly a 
power law on a log-log plot (Fig. 2.9 (b)). 




PQM Model < q > (GeV 2 /fm) pqm Model ( q > (GeV 2 /fm) 

(a) (b) 

Figure 2. 9: (a ) The energy loss model on a linear-linear scale appears to be losing 
sensitivity; increasing q appears to be less and less effective at decreasing Raa ■ (b ) 
The same plot but on a log-log scale. It is clear from the fit that Raa(q) ~ V*? 1 ^ 2 
and is equally sensitive at all values of q explored here. Figure adapted from [310]. 

Fragility is another Boolean concept. By accepting the use of that word 
one must conclude either one way or another than the generic convolution of 
theoretical sensitivity£g>theoretical precision®experimental precision is fragile 
or not. A more quantitative approach would be to determine how much the 
model©data system constrains physical attributes of the medium. The latter 
is emphasized to underscore the important, but often neglected, physics of 
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mapping theoretical input parameters back to usable information on plasma 
properties. Recent work by Phenix has attempted to rigorously statistically 
quantify the knowledge of the controlling model input parameters to be gained 
from data [310, 311]. For WHDG the result is dN g /dy = 14001$!? and ±|28 
at the 68% and the 95% confidence levels, respectively. While the 1-sigma 
constraint is within ~ 20%, the parameter range consistent with data at the 
2-sigma level varies by a factor of 2; this turns out to be roughly true for all 
the energy loss models. 

Regrettably this calculation was not able to take into account any theo- 
retical error, which has not been extensively researched by any group. Recent 
estimates of the error solely due to the running coupling are not promising; 
see Fig. 2.10. And there is the large systematic error of vastly different ap- 
proximations going into the modeling. This suggests three courses of research 
for theoretical high-p^ physics: (1) search for observables that are more theo- 
retically sensitive, (2) more rigorous estimates of major sources of theoretical 
error and the possibility of minimizing them, or (3) search for qualitative tests 
of theoretical formalisms and approximations. We pursue the latter option in 
this thesis. 

2.3 Introduction to AdS/CFT 

The AdS/CFT conjecture is possibly the most important new theoretical tool 
developed in the past two decades. It claims a correspondence, or duality, 
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Figure 2.10: A simple estimate of the large systematic theoretical error due to the 
running of the coupling. Convolved collisional + radiative energy loss calculations 
[243] of ir° Raa{pt) at Rhic for hxed a s = .3, .4, and .5. Figure adapted from 
[243]. 

between SU(N) field theories defined in d dimensions with Type IIB string 
theory on d + 1 dimensional anti de-Sitter (AdS) space and a 9 — d dimen- 
sional compact manifold [312]. The promise of the conjecture comes from the 
two different analytically tractable limits of the duals: field theories for small 
A = ffYM-N 1; classical supergravity for A ^> 1 and N ^> 1. In this way 
strongly coupled string theory problems are approachable via weakly coupled 
perturbative field theoretic techniques. More important for heavy ion physics, 
strongly coupled field theory calculations, previously only analyzable through 
equilibrium Euclidean numerics on the lattice, are now analytically solvable 
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using classical supergravity. But what makes the conjecture so tantalizing 
also makes it so enigmatic: no proof yet exists, and the different regimes of 
tractable applicability mean one may never be found. 

2.3.1 History of String Theory 

There's a certain irony that string theory is again purported to describe the 
strong nuclear force. Created in the 1960's as a phenomenological model to 
explain the mass spectrum of hadrons, string theory originally had some suc- 
cesses. A number of consequences unseen in Nature falsified string theory as 
a theory of hadrons. Undiscouraged, the practitioners audaciously proposed 
string theory as the theory of everything. Having possibly lost its usefulness 
as a theory in the landscape, the theory is back as a tool to analytically cal- 
culate, appropriately, properties of strongly coupled QCD. Surprisingly, string 
theory has rather humble roots (for a very nice historical overview of string 
theory developments until 1985 from one of its founders, see Schwarz [313] and 
references therein). In the early 1960's there was a successful quantum field 
theory of electrodynamics, but Geoffrey Chew argued that a weak-coupling 
field theory approach was inappropriate for the strong force. Even worse, no 
hadron appeared more fundamental than another: there was a paralyzing nu- 
clear democracy (while Gell-Mann had already published his Eightfold Way 
quarks were considered merely a mathematical construct; the pivotal Ruther- 
fordian Slac measurements had not yet been made). 

The order of the day was to focus on physical quantities, especially the S 
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Matrix of on-shell asymptotic particles, as opposed to off-shell physics. General 
physical arguments of causality and nonnegative probabilities implied unitarity 
and maximal analyticity of the S matrix. Chew and Frautschi argued for the 
further requirement of maximal analyticity in the angular momentum. Then 
the partial wave amplitudes ai(s), become analytic functions a(l, s) of angular 
momentum and center of mass energy with isolated poles, Regge poles. Regge 
trajectories, I = a(s), then give the position of these poles. 

Flush with hadrons from the Bevatron, Acs, and Ps, plots of angular 
momentum vs. center-of-mass energy showed a stunning linear dependence 
with a common slope 

a(s) = a(0) + a's, a' ~ 1.0 (GeV) -2 ; (2.10) 

see Fig. 2.11 for one example. It was argued from the crossing symmetry of 
analytically continued scattering amplitudes that the t-channel exchange of 
Regge poles dominated the high-energy, fixed momentum transfer, asymptotic 
behavior of physical amplitudes. In this way 

A(s, t) ~ f3(t)(s/s ) a{s) ; s^oo,t<0. (2.11) 

The partial wave expansion of this amplitude revealed a tower of Breit-Wigner 
resonances at physical values of I, integer for mesons and half-integer for 
baryons. 

In order to fully characterize the S Matrix theory Chew proposed the boot- 
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Figure 2.11: Chew-Frautschi plot of the Regge trajectories a(M 2 ) for the A and 
£ hadrons. Poles at half-integer values of a correspond to the masses of the narrow 
resonances. Figure adapted from [314]. 



strap model, in which the hadrons themselves are exchanged to create the bind- 
ing strong force; in a self-consistent way the hadrons create the force which 
allows their existence. This model implied a duality: the poles in the s and t 
channels would be identical. Unlike a Feynman-like calculation, including all 
channels would in this case result in double counting. 

Then in 1968 Veneziano [315] proposed an actual formula for scattering 
amplitudes with all the required analytic and duality properties that also had 
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linear Regge trajectories: 



A(s,t) 



T 



A(s,t) + A(s,u) + A(t,u) 
T(-u(s))T(-a(t)) 
r (-a(s) - a(t)) 



(2.12) 



(2.13) 



Soon afterward Virasoro [316] published an alternative form with full s, t, and 
u symmetry: 



There were iV-particle generalizations of these formulae, but most surprisingly 
they were found to be factorizable into a spectrum of single particle states of 
an infinite number of harmonic oscillators {c^}, m — 1,2,.... This suggested 
a tantalizing interpretation: these were not just phenomenological fitting func- 
tions but rather a tree approximation of a full-fledged quantum field theory. 

Unfortunately the Lorentz transformation properties of these operators im- 
plied the existence of negative norm states. By imposing algebraic constraints, 
now named after him, on the oscillator operators Virasoro removed these states 
[317]. Solving one problem led to another, though: this algebra implied that 
the leading trajectory associated with Eq. (2.13) had an intercept at a(0) = 1, 
which meant in turn that in addition to a massless vector the theory admit- 
ted a tachyonic ground state as well. Additionally, unitarity required not the 
usual d = 4 dimensions of hadronic physics but rather d = 26. Finally the 



T 



T(-\a{s))T(-\a{t))T(-\a{u)) 



(2.14) 



T(-\a{s) 



\a{t))T(-\a{s) - \a(u))T{-\a{t) 



\a(u))' 
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intercepts for the trajectories of the tt and p mesons gave them unphysically 
large masses. 

Although there were attempts to resolve these issues and advances were 
made, with the advent of QCD and the Standard Model — and the convinc- 
ing experimental evidence for them — string theory rapidly went out of favor. 
Undaunted the remaining string practitioners abandoned it as a theory of 
hadrons and shot for the moon in 1974 by proposing it as a theory of every- 
thing. By doing so the serious flaws of string theory as a theory of hadrons 
became assets: the massless spin-2 particle is the graviton and an inevitable 
consequence of strings; the theory is free of the usual UV divergences associ- 
ated with point-particle quantum gravity; the massles spin-1 particles of open 
strings are Yang-Mills gauge bosons; and the extra dimension dynamics are 
set by the geometry of gravity. It is interesting to note that this promotion 
of string theory changed the string tension by 38 orders of magnitude, from 
a' ~ 1 GeV~ 2 to a' ~ 1/M 2 , where the Planck mass in AD is M p ~ 1.22 x 10 18 
GeV. 

The string theory story until the 1990's is of less interest to the nuclear 
physicist: SUSY and SUGRA emerged in the late 1970's; superstrings, anomaly 
cancellation in the early 1980's; and the discovery of heterotic strings and the 
compactification on Calabi-Yau manifolds of the mid-80 's, the latter of which 
produced the first set of multiplets reminiscent of Standard Model particle 
classification. This second superstring revolution caused a frenzy of activity 
that continues today. While appreciably closer to unifying the known forces, 
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string theory still faces the daunting issues of the landscape and hierarchy 
problems. The more recent work of the 90's is of greater pertinence for heavy 
ion collisions. It is during this time that string theory ideas were applied to 
black holes. 

2.3.2 AdS/CFT Background 

Holography 

The most important idea to emerge from attempting to consistently treat 
quantum mechanics and black holes is that of holography (for a review of 
holography, and especially its connection with AdS/CFT, see [318]). A holo- 
gram in the vernacular is a 2D object that projects a 3D image of the recorded 
subject; in this way the full information of the 3D subject is encoded on a 2D 
surface. Holography in string theory has the same connotation: given some 
volume the physics of the bulk is determined by the physics of the boundary. 
A standard example of this principle is the derivation of the maximum entropy 
of a volume of space, V. Imagine V filled with mass and just enough energy 
that, should the system collapse into a black hole, its event horizon would be 
exactly the boundary of V, dV. We know that the black hole entropy asso- 
ciated with this horizon is Aj AG, where A is the surface area of dV. By the 
second law of thermodynamics this entropy must be greater than or equal to 
the entropy of the system before the collapse; therefore it gives a maximum 
entropy bound for V. (There is some subtlety skimmed over in this argument. 
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As it will evaporate, the black hole is not really in equilibrium. Nonetheless 
the entropy of the final state of radiation is (9(1) times the entropy of the 
black hole). Clearly AdS/CFT is another example of the holographic princi- 
ple; the correspondence claims the boundary d-dimensional gauge field theory 
physics is the same as the bulk d + 1 dimensional physics of supergravity on 
an AdS 5 x S* 5 background. 

n-dimensional GR 

Another important notion is that of a p brane. These are just black hole 
solutions to classical supergravity, found by extremizing the action. Recall the 
physics of ordinary gravity from general relativity. In n dimensions the action 
is 



where Sh is the Einstein-Hilbert action, Cm is the Lagrange density for the 
matter, and g = det (in this AdS/CFT introduction a mostly plus met- 
ric is used), and R is the Ricci scalar. In these 'unnatural' units [G^] = 
L n ~ 1 M~ 1 T~ 2 \ which leads to a length scale of a quantum theory of gravity, 
the Planck length, 




(2.15) 




(2.16) 
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and a Planck mass, 

/ fcn-3 \ ~> 

M ^W^j • (2 - 17) 

Unlike particle physics for which all quantities are in units of energy, or equiv- 
alently length, quantum gravity is dimensionless. For completeness we note 
that a particle of mass M has a Compton wavelength of A c = h/Mc and an 
n-dimensional Schwarzschild radius 

where 2k 2 = 16irG^/c 3 . Then for M ~ M p , X c ~ R s ~ l p and quantum 
effects must be taken into account of a graviational treatment. 

Variation of the action Eq. (2.15) with respect to the metric [319] yields 
the usual Einstein equations (modulo a subtlety due to varying the derivative 
of the metric on the boundary; see [320]) 



R^-\R9^u = K 1 T llu , (2.19) 



where the symmetric, gauge-invariant energy-momentum tensor is 

T " = - 2 (2 ' 20) 

With these (clearly poor) conventions the gravitational potential for a point 
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particle is [320] 

k 2 Mc 3 1 

(2.21 



2(n-2V (n _ 2) |rt- 3 ' 
where Ud is the volume of a <i- dimensional sphere of radius 1 {oj% = 2tt, bj 2 = 47r, 
etc.). We see that n = 4 gives the usual 1/r potential. Eq. (2.21) yields the 
n-dimensional gravitational force on a test mass m, 

1 ; Vl ; (n - 2V (n _ 2) Irl"- 1 1 ; 

Taking Cm = and a mass M at the origin leads to the usual Schwarzschild 
solution in AD (with G = G^\ and H — c—1 for the rest of the thesis), 

l2 / 2GM\ l2 / 2GMV 1 , 2 2lfl , nno x 
<is 2 = - ( 1 \dt 2 +il J dr 2 + r 2 dVL 2 . (2.23) 



There is a nonsingular event horizon at r = 2GM corresponding to g„ switch- 
ing sign and a true singularity at r = (which can be seen from the invariant 
scalar BP vpa = 48(GM) 2 /r 6 ; other scalars that can be examined for sin- 
gularities are the Ricci scalar R = R^^ and the Ricci tensor squared, R^R^). 

Now consider an n = 4 gravity action including the E&M Lagrangian, 
which is simply a two form field strength squared with a proportionality con- 
stant, £ e &m = -(1/4)F^F^: 
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For a point particle of mass M and charge Q located at the origin variation of 
this action leads to the Reissner-Nordstrom metric, 



ds 2 = -Adt 2 + A^dr 2 + r 2 dn 2 2j (2.25) 
2GM GQ 2 

A = 1 — + (2-26) 



There are now two special radial coordinates at which g rr = 0, 



r± = GM± ^G 2 M 2 - GQ 2 ; (2.27) 

there can be zero, one, and two solutions to this equation. For GM 2 < Q 2 
there are no solutions, which leads to a naked singularity at the origin (under 
very general conditions, the cosmic censorship conjecture prevents these from 
forming from gravitational collapse). For GM 2 > Q 2 there are two event 
horizons, both of which are removable coordinate singularities. The r + = r_ 
single solution case defines an extremal black hole. 

p branes 

Following [321, 322] consider the bosonic part of a supergravity action in 10 
dimensions, a slight generalization of Eq. (2.24): 



Chapter 2: Motivation 



55 



where is the dilaton field, F p+2 is a p+2 form field strength, and l s is the string 
length; the string length l s is related to the string tension (2na')~ l by l 2 s = a'. 
This form is suggestive of the result = 87r 6 c/^Zf [320], where g s is given 

by the vacuum expectation value of the dilaton field, g s = exp(0). Just as a 
dimensional point charge sources a 1 form potential A which creates a 2 form 
field strength F = dA in E&M, this p + 2 form is sourced by a p dimensional 
object, a p brane. Similarly to E&M, setting a spherically symmetric charge 
at the origin the integral of the flux of the p + 2 form yields the charge: 



Extremizing the above action for a p brane with the p parallel dimensions 
denoted by x\ the radial distance away from the p brane denoted by p, and 
the angular dimensions encoded in dfl, yields the metric 




N. 



(2.29) 




+ P 2 f4p)^dnl P , 



(2.30) 



9s 



(2.31) 
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r± characterize the mass and charge of the p brane by 

M = T, vhr, 2/8 ((8-P)^ P -^ P ) (2-33) 

N = -^(r + r_)V (2.34) 

where d v is a numerical factor given by 

d P = 2 5 "%^r (1^) . (2.35) 

The interpretation for p branes is exactly the same as for charged black holes: 
r_ > r + gives a naked singularity at r = r_ ; r + > r_ a horizon and removable 
coordinate singularity at r = r + . 

For an extremal p brane r + = r_ and Eq. (2.30) becomes 

ds 2 = y/f+(p) i^-dt 2 + dx*dx*j + f+(p)%~^dp 2 

+p 2 f+(p)^dni_ p . (2.36) 

Changing coordinates to a distance beyond the horizon 

r 7- P = p 7-p _ r J-P (2.37) 
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we find the metric takes the form 



ds 2 

H(r) 

J-p 



-dt 2 + dx i dx l 



+ ^H^)(dr 2 + r 2 dn 2 8 _ p ); 



9- 



;H{r) 



3-P 



1 /?"+\ 7 ~ P 

g s d p Nl 7 -t>, 



(2.38) 
(2.39) 

(2.40) 
(2.41) 



which we will find useful later on. 



D branes and Chan-Paton Factors 





(a) 



(b) 



Figure 2.12: (a) Open strings ending on a D brane. (b) D branes as the source of 
closed strings. Figure adapted from [322]. 



Besides strings, string theory necessarily involves nonperturbative objects 
upon which open strings may end; these are D branes, D for Dirichlet, or 
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sometimes Dp branes (for a review of D branes see [323]). It is thought that 
p branes and Dp branes are different descriptions of the same object, with 
the p brane description of the backreaction on the geometry valid for small 
curvature in Planck units (r_/Z p = R/l p ^> 1). The endpoints of open strings 
can be charged with a non-dynamical degree of freedom. Each end is indexed 
by i or j running from 1 to N, where N is the number of D branes (only 
Type I string theory has open strings not ending on D branes, but we are not 
interested in this string theory); see Fig. 2.13 (a). The NxN matrices A^ form 
a basis for this part of the string wavefunction, and in a scattering amplitude 
each vertex carries one of these matrices. Consider a scattering amplitude 
in which the N D branes are coincident. In this case one must sum over 
all possible Chan-Paton indices and, in addition, since the right endpoint of a 
string is the same as the left endpoint of another string — see Fig. 2.13 (b) — the 
total scattering amplitude results in a trace over these matrices Tr (A ai ■ ■ ■ A ™). 
Therefore the amplitude is invariant under U(N) transformations of these 
matrices, A — > U\U~ l . One may then elevate this global symmetry to a local 
one in spacetime by defining the vertex operator = X^dtX 1 * exp(ik ■ X) 
that transforms under the adjoint of U(N). Since U(N) = U(l) x SU(N)— 
U G U(N) is an N x N unitary matrix. Unitarity implies that its determinant 
is a pure phase: UU* — 1 =>• \det(U)\ = 1. Therefore any U(N) matrix may 
be written as U — exp(i<ft)SU, where SU G SU(N) — we now have a p + 1 
dimensional SU(N) field theory defined on the stack of Dp branes. 

Specifically, the low energy effective action for a D3 brane in AdS space can 
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(a) 

Figure 2.13: (a) Open string with i and j Chan-Paton degrees of freedom, (b) 
Open string scattering amplitude; the Chan-Paton index for the right end of each 
string becomes that of the left end of the next. In the hadronic string theory this 
corresponds to 3 meson to 2 meson scattering. Figures adapted from [323, 324]. 

be reliably approximated by the Born-Infeld action [322, 325] for the massless 
modes of the string theory. For Type IIB the boson part of the action on the 
D brane worldvolume reads: 



S 



1 



D3 



(27r) 3 g s a 



r2 



'4„ t-1 



d*xf 



f 



det( Va p + fd a rdfsr + r 2 f gij d Q 6^ + 2vr«' v / 7^) - 1 
Ang s a' 2 N 



(2.42) 
(2.43) 



where the 9 l are the angular coordinates of the 5-sphere. The r and 9 l are 
then six scalar fields that are functions of the brane coordinates. Notice that 
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by changing coordinates to U — r/a' the a' dependence drops out of Sr> 3 . The 
controlling parameter for Eq. (2.42) then turns out to be g s N(dU) 2 /U 4 [322]. 
The lowest order term from the action yields M = 4 SU(N) super- Yang-Mills 
from the F a p term. The quadratic term has no quantum corrections and the 
quartic term has only a one-loop correction, consistent with M = 4 SYM 
[326, 327]. This loop correction has been calculated from the gauge theory 
and string theory, and the two agree [328]. Moreover, it can be argued that 
all higher order terms in Eq. (2.42) are determined from the fourth-order term 
[312]. In order for this series to make sense g s N(dU) 2 /U 4 <C 1. In particular 
in the supergravity regime of the AdS/CFT duality g s N 3> 1, and the higher 
order terms beyond SU(N) may become important. 

2.3.3 Motivation for the AdS/CFT Correspondence 

In this section we will follow [322] closely for a schematic motivation of the 
AdS/CFT conjecture (see the same for a very good review of the correspon- 
dence as well as an extensive bibliography). Consider from a field theoretic 
standpoint Type IIB string theory in a 10 dimensional spacetime with a stack 
of TV coincident extremal (r + = r_ = R) D3 branes at the origin. At ener- 
gies small compared to the string scale l/l s = 1/ \fa' only the massless string 
modes can be excited. The full low energy effective action is then 



S — Sbulk + SWanes + S int . (2.44) 
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<Sbuik is the 10D supergravity action with higher derivative corrections. Shanes 
is the 3 + ID SYM action with higher derivative corrections. Finally S- m t 
describes the interaction between the brane and bulk modes. It turns out 
that all of these actions are controlled by a': the bulk action is controlled 
by Newton's constant k, the brane action by a' 2 , and S- m t also by k. Then 
taking the low energy limit k ~ g s <y' 2 — > (specifically a' — > with U = r/a' 
fixed), the higher derivative and interaction terms all drop out. We are left 
with noninteracting supergravity in the bulk and M = 4 SU (N) SYM on the 
branes, which are decoupled from each other. 

On the other hand consider the geometrical interpretation of the low energy 
limit of Type IIB string theory in 10 dimensions. From Eq. (2.38) we have 
that 



ds 2 = f- 1/2 (-dt 2 + dxi + dx 2 2 + dx 2 3 ) + f 1/2 (dr 2 + r 2 dVLl), (2.45) 



There are two regions of importance: the near horizon limit, r — > 0, and 
the asymptotically flat limit, r — > oo. In the flat limit the metric is the 
usual Minkowski one corresponding to, because of a' — > 0, free supergravity. 
In the low energy limit these massless excitations decouple from the brane 
physics as their cross section goes ~ oj 3 R 8 [329] (note that cross sections 
in n spatial dimensions for problems with p translational symmetries have 
dimension L n ~ 1 ~ p ). On the other hand the nontrivial gu component of Eq. 



/ 



1 + 




R 4 = 4irg s l 4 s N. 



(2.46) 
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(2.45) means that excitations in the near horizon throat region appear red 
shifted to an observer at infinity, E- lv f = f~ 1 ^ 4 E r . Since it is the energy as 
measured by an observer at infinity that is important, in the limit r — > 
the full Type IIB string theory must be kept. Nevertheless the higher energy 
modes from the string theory cannot escape the throat region without being 
redshifted away. We are thus left with supergravity in flat asymptopia and 
IIB string theory compactified on the near-horizon geometry, and the two are 
decoupled. 

Stepping back we find that we have two pictures of the same low energy 
limit of one theory, Type IIB string theory. In the field theory picture we 
have two decoupled theories: supergravity in the far region and M = 4 SU(N) 
SYM on the D branes. In the geometry picture we also have two decoupled 
theories: supergravity in asymptotically flat space and Type IIB string theory 
in the throat region. Noticing that both pictures have identical asymptotic 
supergravity in them we boldly propose, in the spirit of Maldacena, that the 
other decoupled theories are also identical: Type IIB string theory compact- 
ified on the near horizon background of Eq. (2.45) is dual to 3 + ID M = 4 
SU(N) SYM. 

The regions for which analytic tools exist for these two different pictures 
turn out to be completely incompatible. Comparing the Born-Infeld action of 
iV coincident D p branes and that for an SU(N) two form field strength yields 

g 2 SYM = 2g s (2Tr) p -H 3 - p g\ YM = Ang s for p = 3. (2.47) 
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Perturbative field theory is a consistent approach for calculations when 

9 2 symN = 4ntg 8 N = ^ « 1, (2.48) 

'■s 

where the final relation came from Eq. (2.46). 

Contrariwise string theory is tractable in the classical supergravity limit, in 
which case the characteristic length scale of the problem, R, is large compared 
to the two length scales of quantum string theory: the string length, l s , and the 
Planck length, l p ; this ensures no string corrections and no loop corrections, 
respectively. Demanding this implies 

* = A7rg s N = g 2 SYM N » 1; ^ ~ N » 1. (2.49) 

There is an extra subtlety due to the SL(2, Z) self-duality symmetry of 
both the gauge and string theory. For the coupling r = 47u/g| yM — i/g s the 
theories are invariant under r — > (ar + b)/ (cr + d) where ad — be = 1. We 
can see that one can take g s — > l/g s with a = d = and b = 1, c = —1; a 
strong string coupling can be swapped for a weak one and vice-versa, but the 
number of colors must always remain large. This is the origin of the A large 
and fixed while iV — > oo limit; this guarantees an S transformation leaves R 
large compared to both l s and l p . 

Now we will show that the near horizon limit that is the background for 
the dual string theory is AdS$ x S 5 . To see this, taking r — > in Eq. (2.45) 
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yields 



r 2 R 2 
ds 2 = -^(-dt 2 + dxl + dx\ + dx\) + — dr 2 + R 2 dQ 2 5 . (2.50) 



Compare this to the metric of a p + 2 dimensional hyperboloid 

p+i 

X ° 2 + X p+22 -J2 X ' 2 = R 2 (2-51) 

8=1 

embedded in a flat space with SO(2,p+l) isometries 

p+i 

ds 2 = -dX° 2 - dX p+22 + J2 dxi2 ■ ( 2 - 52 ) 

i=l 

If we reparameterize the hypersurface with coordinates y a = (t,r,x), r > 
such that 



x° = 


1( 

2r \ 




rx % 
~R' 


X p+l = 


if 

2r \ 


X P+2 = 


rt 



(i = l...p) (2.54) 
i? 2 -^( J R 2 -^ 2 + t 2 )), (2.55) 



(2.56) 



then the induced metric is obtained from the pullback [319] 



{y * 9U = %*% 9 ^ (2 ' 57) 

2 P+ 1 p2 

= 4(-^ 2 + "E^) + ^ r2 - (2-58) 
R z ^-^ r 1 

i=l 
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Therefore Eq. (2.50) is exactly AdS 5 x S 5 . 

While there is no proof yet of the AdS/CFT conjecture a number of tests 
have found quantities calculated from the field theory side and the string theory 
side agree. These are nontrivial in the sense that they must be independent 
of coupling. These tests include checking the equivalence of symmetries in 
both theories; some correlation functions, usually associated with anomalies, 
that are independent of quantum corrections and A = g§y M iV; the spectrum 
of chiral operators match for those that are currently calculable; see [322] and 
references therein for a more complete list. 

Having used the above picture to motivate the correspondence we now 
throw away that geometry and simply take our whole space to be AdS§ x S 5 . 
It is useful to know that the boundary at asymptotic r — > oo is the usual 
3 + ID Minkowski. To see this first change coordinates to y — R/r. Then Eq. 
(2.50) becomes 



ds 2 



= —{-dt 2 + dx 2 ) + —dy 2 + R 2 dQl 



(2.59) 




(2.60) 



where for the second step we Weyl rescaled the metric by an overall factor of 
y 2 /R 2 . Now r — )• oo corresponds to y — > 0, and we see that at the boundary 
the radius of the sphere shrinks to and usual 3 + ID Minkowski is recovered. 
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2.3.4 Applications of AdS/CFT to Heavy Ion Collisions 

To make contact with thermal QCD physics we extend the conjecture to nonex- 
tremal D3 branes, in which case conformality is broken. Starting with the 
nonextremal 3-brane metric Eq. (2.30), take R = r_, r = rl— ri, and change 
coordinates to r 4 = p 4 — i? 4 ; it then becomes 



ds 2 = 1 \-h(r)dt 2 + dx 2 ] + ^/W{ 
y/H(r) 



2 



h(r) 



(2-61) 



H(r) = 1 + ^; h(r) = l-£. (2.62) 



This yields in the near-horizon, or decoupling, limit of r <C R 



ds 2 = ^ \-h(r)dt 2 + dx 2 ] + -^r-Ar 2 + R 2 dQl (2.63) 
R l L J r 2 n{r) 



We find the temperature associated with the horizon by requiring that the 
t = it coordinate after Wick rotation has the correct periodicity for a path 
integral partition function. We could integrate r in a circle around r = r$ [330] 
but will rather follow [331]: for a metric of the form ds 2 = Vdr 2 + V~ 1 dr 2 + . . . 
the temperature is T = V'\ r=ro /4:7i, where r is the radial position of the 



horizon. Both methods yield 



T=^. (2.64) 



Using this equation for T one may follow [330, 332] to compare the entropies 
of weakly coupled Af = 4 U(N) and strongly coupled gravity at the same 
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temperature. The first comes from the usual counting of states with 6N 2 
massless scalars and 4N 2 massless Weyl fermions. Taking the spatial volume 
of the 3-brane to be V3 we have 



27T 3 

S = —N 2 V 3 T 3 . (2.65) 



The Bekenstein-Hawking relation, Sbh = A/4Gn, gives the entropy of the 
black p brane in the supergravity picture. The 8-dimensional area of the 
horizon may be read off of Eq. (2.63) as 

A = g) 3 V 3 R 5 Q 5 = n 6 V 3 R 8 T 3 , (2.66) 

where the second equality comes from Eq. (2.64). By Eq. (2.46) GSJ 0) = 
8Tt 6 g 2 l 8 s = tt 4 R s /(2N 2 ), so 

Sbh=^ = yiW 3 T 3 . (2.67) 

The factor of 4/3 difference between the weakly-coupled and black hole entropy 
limits is a prediction for the change in the number of degrees of freedom from 
weakly-coupled SU(N) to strongly-coupled SU(N). This provides a natural 
explanation for the discrepancy of about 20% between lattice results at a 
few times T c , still at strong coupling, and the expected weak-coupling Stefan- 
Boltzmann limit as seen in Fig. 1.6 (b). 

One may next find the strong coupling limit for the entropy to viscosity 
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ratio, r]/s. By comparing the correlation function for the low energy scatter- 
ing cross section used in the motivation for the correspondence to the Kubo 
formula for the viscosity of a strongly-coupled SYM plasma [333] one finds 
that 

^abs(O) A 

71 = l^G~ N = To^GV (2 - 68) 
where [334] found that the cross section for extended p-branes has a universal 
low energy limit of the area of the horizon. Using this with Eq. (2.67), one quite 
trivially arrives at the famous strong coupling limit of viscosity to entropy 



S 47T 



The result led to another conjecture [335], that this is in fact a universal lower 
bound. As noted in Section 2.1.1, qualitative estimates of the viscosity at 
Rhic [121] appear to require just such values of rj/s ~ .1 ~ l/4ir. This gives 
great impetus to look for additional observables calculable in AdS/CFT. 

In fact we can learn a great deal by extending the AdS/CFT conjecture 
further by breaking some of the field theory supersymmetry encoded in the 
S 5 part of the geometry. Specifically the introduction of M <C N D7 branes 
that fill space from the asymptotic 3 + ID Minkowski boundary to some finite 
distance above r = and wrapping an equatorial S 3 C S 5 corresponds to 
introducing M fundamental flavors in the gauge theory [336]. The M ^ N 
requirement allows the D3 branes to backreact and warp spacetime in the 
large N limit while the D7 branes act only as probes. An open string ending 
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on the r < oo boundary of one of these D7 branes corresponds to a heavy 
fundamental quark. At zero temperature the geometry supports mesons, with 
an open string having both its endpoints on D7 branes. At finite tempera- 
ture T a string may hang from the D7 brane and terminate on one of the 
nonextremal D3 branes. This is identified with an open heavy quark in a 
deconfined SYM plasma held at T. One may calculate the Debye screening 
length for this process, see Fig. 2.14 (a), [337] (and references therein) and also 
find a small enhancement to qq dissociation for moving heavy mesons such as 
J/ty [338, 339]. Given some initial velocity a hanging string will both perturb 
the Minkowski boundary metric, leaving a wake behind its path in the SYM 
plasma complete with Mach cones [340, 341] (see Fig. 2.14 (b)), and also lose 
momentum from its D7 open heavy flavor endpoint into the D3 brane horizon 
[342-344]; see Fig. 2.15. The big payoff from introducing this new complexity 
is analytic formulae for these intuitively familiar effects. 

Other applications of the conjecture include calculating the light-like Wil- 
son loop [338, 345] conjectured to be connected with the jet quenching pa- 
rameter q [265] via AdS/CFT, although there is controversy over its ultimate 
magnitude [346-350]. Also the conjecture was used to investigate the diffusion 
coefficient in a Langevin formalism [351, 352]; this is only valid in the v — > 
limit [347], and it is not clear what the corrections for nonzero velocity are. 
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// 




(a) 



(b) 



Figure 2.14: (a) Heavy qq potential for N c = 3 and A = 10 (blue, dotted), 10 2 
(green, short dash), 10 3 (red, long dash), and 10 4 (black, solid). As A — ^ oo the 
potential develops a kink at |x| = .24/T; for larger distances the potential continues 
to rise, but very slowly. Figure adapted from [337]. (b) The AdS/CFT scaled energy 
density of the SYM plasma with Mach cone for a heavy quark with v = 3/4. Figure 
adapted from [341]. 

2.3.5 Testing the Conjecture 

The most well known result of using this complicated machinery is the famous 
viscosity to entropy bound [332, 335] 



As noted earlier, qualitative estimates of the viscosity at Rhic [121] appear to 
require such small values of rj/s. This encourages continued investigation of 
applications of the AdS/CFT conjecture, but the big question is: how much 
can we trust the AdS/CFT results? The duality of M = 4 SYM and Type IIB 
string theory on a pure AdS^ x S 5 background has strong motivation and passed 




(2.70) 
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Figure 2.15: Visualization of a heavy quark propagating through a thermal 
plasma. The fundamental quark is represented by a string hanging from a D7 brane 
(in this case at r = oo) down through the horizon onto the stack of coincidentD3 
branes. Its motion causes a metric disturbance on the Mikowski boundary dual to 
a wake in the SYM plasma and a loss of momentum down the length of the string. 
Figure adapted from [340]. 

a number of theoretical tests, although it remains unproven. But maximally 
supersymmetric Yang-Mills theory is very far from QCD. As symmetries are 
broken and the field theory dual to string theory becomes more and more 
like QCD, the further the conjecture must be extended and the less evidence 
exists for its holding. And while QCD is far from a quantitative description of 
heavy ion collisions there is a well known procedure for estimating the errors 
involved. On the other hand, calculating quantities in a theory that is similar 
but not the same as QCD is an uncontrolled approximation, and it is far from 
clear how much those results will be modified by convergence to QCD. 

We are thus motivated to find qualitative features of heavy ion collisions for 
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which AdS/CFT predictions are not and can not be reproduced by weakly- 
coupled QCD. In this sense the entropy to viscosity ratio is the worst ob- 
servable to use: the Frankfurt parton cascade group finds pQCD also gives 
rj/s ~ .1; and the value of the ratio, due to the uncontrolled heavy ion ini- 
tial conditions that so strongly affect hydrodynamics results and the unknown 
importance of viscous effects in hydrodynamics predictions, is very poorly de- 
termined. Nevertheless AdS/CFT can be applied to heavy flavor jets and 
jet-bulk interactions, and we will find that these are falsifiable predictions 
qualitatively different from pQCD. 

Specifically AdS / CFT can be applied to the energy loss of heavy quarks 
propagating through a thermalized plasma [342, 343]. It turns out that the 
drag experienced by these quarks is independent of their momentum and in- 
versely proportional to their mass; this is entirely different from the logarithmic 
dependence on both mass and momentum found in pQCD. A major result of 
this thesis is to propose an experimental observable, the ratio of charm to bot- 
tom quark Raa{pt), as a robust observable for differentiating between these 
two theories; see Chapter 5. 

Another success of the AdS/CFT paradigm came from the consistency 
check of Mach cones produced in the bulk matter by heavy quarks propagating 
faster than the speed of sound, v = It turns out though that connecting 

this phenomenon to the observed conical flow in dijet correlations at Rhic is 
highly nontrivial [353]. In fact one finds that hadronization removes all traces 
of a signal from any Mach wake in the plasma; the temperature variations are 
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simply too small. Recent work has surprisingly shown that the observable con- 
sequences from the so-called 'neck' region, the bulk plasma trailing behind the 
quark jet, are qualitatively different for thermodynamic sources from pQCD 
and from AdS/CFT [354]; see Fig. 2.16. 
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Figure 2.16: The away side correlation functions of particles produced by Cooper- 
Frye hadronization associated with a high-p? jet propagating through a thermalized 
bulk plasma. The hydrodynamic source on the left comes from pQCD calculations 
[355]; the right from AdS/CFT [340]. Figure adapted from [354]. 



2.4 Outline of Thesis 

In the next chapter a simultaneous qualitative description of intermediate-pT 
Raa(pt) and V2(pt) is developed predicated on a nonperturbative momentum 
punch at the boundary of the QGP and the vacuum. Another interesting 
source of this focusing of jets could come from the reduction of the Debye 
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mass rri£) ~ gT as the parton propagates through a cooling plasma. This 
unexpected and heretofore inexplicably large vi given the overall suppression 
is argued as an experimental signal for deconfinement. 

In Chapter 4 we reduce the qualitative disagreement between pQCD calcu- 
lations of heavy quark energy loss and observations of light and heavy meson 
suppression from Rhic. The inclusion of the elastic energy loss channel, whose 
neglect disproportionately affects the heavy quarks, increases the quenching 
of all parton flavors. For single fixed pathlength calculations with the usual 
L ~ 5 fm, pion predictions are falsified by data. The realistic treatment of ge- 
ometry is crucial for a simultaneous description of both pions and nonphotonic 
electrons: just as fluctuations in the number of emitted gluons softened the 
energy loss so does allowing fluctuations in the in-medium distances traveled 
by jets. 

The proposal of taking the ratio of charm to bottom quarks as a robust 
test of pQCD and AdS/CFT ideas is the focus of Chapter 5. Most clearly seen 
at Lhc, scheduled to run heavy ions soon, the ratio goes to 1 in all pQCD 
predictions as the mass scale of the quark becomes irrelevant compared to the 
momentum scale of the particle. On the other hand the ratio is flat and much 
smaller than one for AdS/CFT drag calculations. The critical issue comes 
from the momentum regimes of applicability: the momentum above which the 
drag approximations fail is not clearly known, nor is the momenta below which 
perturbative predictions break down entirely understood, either. Upgrades at 
Rhic will allow for discrimination between charm and bottom jets, possibly 
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ending the controversy over the contributions to heavy quark energy loss; the 
advantage there will come from lower multiplicities so that the drag regime 
will extend to higher momenta. 

Finally we discuss heavy quark production radiation in Chapter 6. By 
including the influence of the away-side jet we find corrections to previously 
published Ter-Mikayelian effect results. This is also the first step in computing 
the photon radiation associated not only with 2 — » 2 processes but also from 
momentum changes due to gluon bremsstrahlung. 



76 



Chapter 3 

V2 vs. Raa as a Constraint on 
sQGP Dynamics 

3.1 Introduction 

A complete theoretical description of Rhic mid- to high-p T particles must 
reproduce Raa(<P), the ratio of observed high momentum leading hadrons in 
A + A collisions to binary scaled p + p collisions as a function of the angle 
made with respect to the reaction plane, 0. A useful means of representing 
this quantity is via Fourier expansion: 

Raa{4>, Pt, y) = Raa(pt, y) (1 + 2v 2 (p T , y) cos(20) + ...)• (3.1) 

The goal of theorists, then, is to correctly replicate these moments, the most 
important being the first two: the normalization and the azimuthal anisotropy. 

This is actually quite a difficult problem due to the anticorrelated nature 
of Raa and v 2 . In order to reproduce the anisotropy of the jet data, the 
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partons must be made sensitive to the medium. Then the anisotropy of the 
background created by noncentral RHIC collisions translates into anisotropy of 
the jets. However, this sensitivity to the medium also leads to jet attenuation 
and a decrease in Raa- As we will see, the very high t> 2 values measured 
at Rhic mean that, at best, previous models either oversuppressed Raa or 
underpredicted v 2 . 

Hydrodynamics has been applied to the \ow-p T particles with great success. 
Unfortunately, the lack of thermal equilibrium precludes the use of hydrody- 
namics in the momentum regime we are interested in. Moreover, a naive 
application would highly oversuppress Raa due to the Boltzmann factors. 

Parton transport theory attempts to extend hydrodynamics' range of ap- 
plicability to higher transverse momenta [135, 138]. The Molnar parton cas- 
cade (MPC) succeeded in describing the low- and intermediate-pr f 2 results of 
Rhic only by taking the parton elastic cross sections to be extreme, a t ~ 45 
mb [138]. Even more damning, we will show that the MPC approach reaches 
the reported t>2 values at the expense of reproducing the experimental Raa 
data. 

pQCD becomes valid for moderate and higher pt partons, and models 
based on pQCD calculations of radiative energy loss have had success in re- 
producing the experimental Raa{.Pt) data [356]. These models use a single, 
representative pathlength; as such, they give t>2 = 0. 

Compelling but not convincing, past geometric approximations of energy 
loss sought to demonstrate the inability of pure jet quenching to match the 
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Rhic v 2 findings [191, 357, 358]. In [357], the author found that v 2 reached a 
maximum value that was smaller than the data at the time. Our formulation 
achieves large enough values of v 2 , but over-quenches Raa to do so. In [191], 
the authors found that when they fit their model to the most central Raa data, 
they could not replicate the large v 2 observations. Instead of modeling Raa or 
v 2 separately, we ask a different question: when considering the experimental 
Raa and v 2 on equal footing, is there any way a pure geometric energy loss 
model can be consistent with the data within the new, 2D error ellipses? We 
find the answer is no. 

[191, 359] demonstrated the powerful influence of the nuclear density ge- 
ometry used in calculations. Sharp-edged approximations such as the hard 
sphere (HS) and especially the hard cylinder (HC) greatly enhance model v 2 
results. Thus past works that used these geometries and could reproduce both 
Raa and v 2 , such as energy loss with hydrodynamics [359], energy loss with 
thermal absorption [360], and energy loss with low-pr flow effects [361] will 
no longer produce large enough v 2 values in realistic, diffuse medium densi- 
ties. We will investigate the case of energy loss with thermal absorption more 
thoroughly later in this paper. 
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3.2 Plot of vi vs. Raa an d Failures of Previous 
Models 

In Fig. 3.1, we combine Star charged hadron Raa(pt) [362] and v 2 {px) (4 
particle correlations only) [363], Phenix charged hadron Raa(pt) [364] and v 2 
[108] vs. centrality, and Phenix it (p T > 4 GeV) Raa(<P) [365] centrality data. 
We naively averaged the Star and Phenix Raa{pt) results to approximately 
match the pr bins of their corresponding v 2 measurements. We report the 
Raa and v 2 modes for the Phenix tt° Raa{4>) data. The error estimates are 
unfortunately schematic only. 

The charged hadron data from both Star and Phenix appear to fall on a 
single line extending from low Raa and v 2 to high Raa and v 2 . The v 2 values 
are larger in general than one would expect a priori, and the absence of a 
decline in v 2 as Raa — >• 1 is especially surprising. For the tc° data, the possibly 
flat in Raa trend of v 2 is also quite strange. Taken as a conglomerate with 
the charged hadrons, however, the n° Raa(4>) data appears consistent within 
error. 

The MPC has a single free parameter: the opacity, X = J dza t p g [138]. 
Allowing the MPC to run at an approximate 30% centrality and varying the 
opacity gives the dashed curve in Fig. 3.1. We see that this curve does not pass 
near the area in the v 2 vs. Raa diagram corresponding to the experimental 
results for the same centrality. While in [138] a single x predicts v 2 (pt) well, 
the MPC fails at simultaneously reproducing the first two moments of Raa{4>)- 
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In order to make a rough calculation of both Raa and resulting from 
partonic energy loss to the medium, we approximate the first order in radiative 
energy loss [268, 274, 281] with a purely geometric model that neglects number 
fluctuations. The asymptotic approximation of this loss in a static medium 
for a parton of energy E traversing a path length L in the limit of EL ^> 1 
and E ^> fi is 



(3.2) 



IE 

In — - 0.048 



E AXE 

where \l is the Debye screening mass, and A is the mean-free-path of the 
parton [283]. Taking the bracketed term to be constant (this is not such a 
good approximation: it changes by a factor of 2 for the range of parameters 
considered here) and p part oc dN g /dy (note that p 2 /A oc dN g /dy), we use an 
energy loss scheme similar to [191]: 

e = = kI. (3.3) 

E 

k is a free parameter encapsulating the E dependence, etc. of Eq. (3.2) and the 
proportionality constant between dN g /dy and p part . I represents the integral 
through the Bjorken expanding medium, taken to be 

I= dl 1 TT~TPp^ x + ( l + l ^ n ^ y + ( l + l o) n y), ( 3 - 4 ) 
Jo 1 + '0 

where we use / = .2 fm as the formation time of the medium. We consider 
only ID expansion here because [359] showed that including the transverse 
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expansion of the medium had a negligible effect on the total energy loss. 

Since the partonic production spectrum observed at Rhic can be approx- 
imated by a power law, we use the momentum Jacobian (p£ = (1 — e)p l T ) as 
the survival probability of hard partons; see Appendix B. Following the fac- 
torization theorem and Glauber approach [366], we distribute partons in the 
overlap region according to p co n = Taa an d isotropically in azimuth; hence 



where 4 < n < 5. The difference between n = 4 and n = 5 is less than 10%, 
and in this paper we will always use the former. 

We evaluate Raa{<P) at 24 values of from 0-2tt and then find the Fourier 
modes Raa and v 2 of this distribution. We label the results of this model in 
the plots as GREL for geometric radiative energy loss. A different method 
for finding v 2 , not used here, assumes the final parton distribution is given 
exactly by Raa and v 2 , and then determines v 2 from the ratio of Raa(0) and 
Raa{k/2); this systematically enhances v 2 , especially at large centralities. 

Choosing 40 — 50% centrality as a good representative, we calculate the 
line in (Raa,Pt) space corresponding to < k < oo, Fig. 3.1. One sees in the 
figure that the curve does not cross the error ellipse for the relevant Phenix 
7t° data. This is despite the simplifying use of hard sphere geometry, where 
we took Rhs = 6.78 fm to ensure that (r\ ws ) = (rj_ HS ), that amplifies 
the t>2. Since the experimental errors are the most basic estimate — it would 



Raa(4>; b) 



JJ dxdy T AA (x, y; b) (1 - e(x, y, 0; 6)) 
N co n(b) 



ii 



(3.5) 
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be very interesting to see the results of a careful analysis — and the theoretical 
models do not currently have any error estimates associated with them prevent 
a statistical quantification of the disagreement, one is strongly persuaded that 
neither of these models encompasses all the relevant physics. 




R AA 



PHENIX p T =4 QeVrP I STAR Pt=5 GeV h 1 

PHENIX p T =4 GeV h 1 ■ STAR p T =? GeV h* 

■ STAR p T =4GeVh ± 

Figure 3.1: Star /i± data for 0-5%, 10-20%, 20-30%, 30-40%, and 40-60%, Phenix 
/i ± data for 0-20%, 20-40%, and 40-60%, and Phenix vr° data for 10-20%, 20-30%, 
. . . , 50-60% centralities. MFC results for increasing x that should go between the 
20-30% and 30-40% Phenix vr° results. GREL results for increasing k that misses 
the 40-50% Phenix tt° data. 
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3.3 Exclusion of Detailed Balance 

In [283] , Wang and Wang derived the first order in opacity formula for stimu- 
lated emission and thermal absorption associated with the multiple scattering 
of a propagating parton in a hot, dense QCD medium. They found that this 
effects a net reduction of the energy loss. Since their energy gain has a linear 
length dependence in a static medium, as opposed to the quadratic dependence 
for radiative energy loss, one might hope to reproduce the experimental Raa 
and V2 data by increasing the opacity enough that radiative energy loss creates 
the t>2 and the energy gain properly reduces the oversuppression of Raa- 

From [283] , we have that in the limit of EL 3> 1 and E ^> the asymptotic 
approximation of the effect in a static medium to the first order in opacity is: 



AE2 _ , 7ra s C R LT* 



E 3\E 2 



/i 2 L i 6C'(2) 



In -=- - 1 + 1e 

T 7T 



2 



(3.6) 



Since /i 2 = 4ira s T 2 and the bracketed term is actually approximately constant 
for our range of inputs, we use 



AE rad AE abs 

— — = kI- kl 2 , (3.7) 



where kI is the same as in Eq. (3.3), k is a free parameter encapsulating the 
proportionality constants in Eq. (3.6), and I 2 represents an integral through 
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the ID expanding medium: 

poo ^ 

h= dl 7 -^ r p part (a;+ (I + l )n x , y + (l + l )n y ). (3.8) 

Note that I2 has one less power of / in the integrand; this allows a unique 
determination of the two free parameters (k, k) when fitting the model to a 
particular (Raa, v%) point. 

Using the same hard sphere geometry as before, we are able to fit the 
20-30% centrality Phenix tt° Raa{4>) data point with k — .5 and k = .25. 
Keeping k and k fixed, Fig. 3.2 shows the results when the impact parameter 
is varied (long-dashed, blue curve); one can see that this matches the Rhic 
Raa and vi trends quite well. 

Taking Eq. (3.7) seriously, we invert the relationships and solve for dN g /dy. 
Recalling that p = gT, A -1 ~ 97ra^(p/2) 1 / 3 /4, and the density for an expand- 
ing medium is p = -^-j (A±(^)), we find that, as a function of parton energy 
E, path length L, formation time Iq, and a s , 

dN ra d AE kL 

dN abs 4E 2 lL 

~ k o^ , r ^Part, (3.10) 



dy 3nCRalv2 lo + L 

where v\ and V2 correspond to the bracketed terms in the energy loss and 
energy gain approximations, Eqs. (3.2) and (3.6), respectively. 

For our fitted values of k and k, the choice of E = 6 GeV, L = 5 fm, 
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and a s = .4 gives dN r g ad /dy ~ 1000 from Eq. (3.9) and dNf s /dy ~ 3000 
from Eq. (3.10) in 0-5% most central collisions. Not only are the two values 
inconsistent, but the dN g /dy value needed to create the necessary absorption 
is rather large. Keeping the other quantities fixed but changing E to 10 GeV, 
one finds dN g /dy ~ 1000 from Eq. (3.9) and dNf s /dy ~ 9000 from Eq. 
(3.10). The huge increase of dN g /dy to values inconsistent with the Rhic 
entropy data reflects the E 2 dependence of the Detailed Balance absorption. 
It seems the only way to have a large enough energy gain while maintaining 
a dNg/dy ~ 1000 is to increase a s above 1. Moreover, these calculations 
were performed using a hard sphere nuclear geometry profile, which naturally 
enhances v%. For the more realistic Woods-Saxon profile, k would have to 
increase significantly. Correspondingly, k would rise to counter the quenching 
of Raa, aggravating the problem of unseemingly large dNg bs /dy. 

3.4 Success of the Punch 

Building on the success of radiative energy loss in reproducing Raa{pt), and 
supposing that latent heat, the bag constant, the screening mass, or other 
deconfinement effects might provide a small (~ 1 GeV) momentum boost to 
partons in the direction normal to the surface of emission, we created a new 
model based on the GREL model that includes a momentum "punch," Ap?- 
After propagating to the edge of the medium with GREL, the parton's final, 
"punched-up" momentum and angle of emission are recomputed, giving a new 
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PHENIX p T =4 QeVrP ■ STAR p T =5 GeV h 1 

PHENIX p T =4 GeV h 1 ■ STAR p T =? GeV h* 

STAR p T =4 GeV h ± 



Figure 3.2: Addition of thermal absorption (dashed blue curve) or momentum 
punch (solid red curve) to GREL; models' free parameters fit to Phenix 20-30% 
centrality ir° Raa{4>) ■ 

probability of escape. Fitting to a single {Raa-,^) point provides a unique 
specification of k and punch magnitude. The results are astonishing: one sees 
from Fig. 3.2 that a tiny, .5 GeV, punch on a 10 GeV parton reproduces the 
data quite well over all centralities. Fitting the Phenix 20-30% tt° data sets 
k = .18 and the aforementioned Apt = .5 GeV. The size of the representative 
parton's initial momentum is on the high side for the displayed RHIC data; 
however, the important quantity is the ratio Ap?/E. Moreover, although 
the geometry used artificially enlarges the t>2, we feel confident that when 
this model is implemented for a Woods-Saxon geometry (a nontrivial task 
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because the "edge" of the medium is then poorly defined), the necessarily 

larger final punch magnitude will still be of order 1 GeV. The magnitude 

of this deconfinement-caused momentum boost must be independent of the 

parton's momentum; hence ^(pt) will decrease like 1/pt- Moreover, since e is 

larger out of plane than in, a fixed Apt enhances Raa{^/^) more than Raa{0)- 

These are precisely the trends seen by Phenix [367]. Keeping the same values 

for k, k, Apt, etc. as for Au + Au, we show in Fig. 3.3 the centrality-binned 

Raa and t>2 results for Cu + Cu in the three geometric energy loss models. 
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Figure 3.3: Cu + Cu results for GREL (black, dotted), GREL with the punch 
(red, solid), and GREL with absorption (blue, dashed). 



We can in fact go further in thinking about the origin of this deconfinement 
mechanism. Specifically we will find the change in the Debye mass of the jet 
fj,i(x,t) = gT(xi,ti) as it propagates through the medium provides the right 
order of magnitude momentum punch. In moving from position i to i + 1 one 
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calculates the deflection of the jet following a slight generalization the usual 
means. At i the parton has momentum /cf = (u, k\\, k±^), where the parallel 
direction is defined by the local line of constant temperature in 2D and the 
perpendicular direction is its normal; the normal is in the direction of the 
2D gradient of the medium temperature, equivalently density. By symmetry 
crossing the plane cannot affect the parallel component of the momentum so 
at i + 1 we have kf +1 = (uu, k\\, k± ji+ i). Then 

fci )i+ i = *ii+^-^ + i- (3-n) 

As the jet moves to the less dense edges of the medium and the whole fireball 
cools as it Bjorken expands /ij+i < /ij, and the jet receives a continuous punch 
in the local normal direction. This will tend to bend the jets, focusing them 
into higher t> 2 and simultaneously increase the Raa- For a quark (gluon) 
produced in the center of a 20-40% centrality collision at t = .2 fm and 
followed until T = T C = 160 MeV the change in Debye mass is 0.38 (0.55) GeV 
which is slightly smaller (larger) than the nominal 0.5 GeV investigated earlier. 
It will be interesting to see if a more quantitative calculation in a realistic 
diffuse medium will meet the expectation that the data can be understood 
from this mechanism. 
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3.5 Conclusions 

By failing to simultaneously match the Raa and v-i values seen at Rhic we 
discounted the MPC and pure GREL models. We showed that while includ- 
ing medium- induced absorption reproduces the Raa{4>) phenomena, it does 
so at the expense of inconsistent and huge dN g /dy. Recent parton cascade 
results including 2 <-> 3 processes were claimed [141] to simultaneously fit 
both Raa(pt) an d v 2 (pt)- Careful examination of the paper reveals that their 
large v 2 comes from using huge a s = .6 and gluon jets only, but the claims of 
consistency with Raa use only a s = .3. In this chapter we showed that the 
addition of a mere 5% punch created a RHic-following trend. The reduction 
in Debye mass from the hot, dense center of the medium to the cool, diffuse 
edge provides just such an order of magnitude momentum change. 



90 



Chapter 4 

Elastic, Inelastic, and Path 
Length Fluctuations in Jet 
Tomography at RHIC and LHC 



4.1 Introduction 

Light quark and gluon jet quenching observed via ir, rj suppression [368-370] 
in Cu+Cu and Au+Au collisions at y/s = 62 — 200 AGeV at the Relativis- 
tic Heavy Ion Collider (Rhic) has been remarkably consistent thus far with 
predictions [274, 275, 277, 297, 371-373]. However, recent non-photonic single 
electron data [237-240] (which present an indirect probe of heavy quark en- 
ergy loss) have significantly challenged the underlying assumptions of the jet 
tomography theory (see [374]). A much larger suppression of electrons than 
predicted was observed in the pt ~ 4 — 8 GeV region (see Fig. 4.1). These data 
falsify the assumption that heavy quark quenching is dominated by radiative 
energy loss when the bulk QCD matter parton density is constrained by the 
observed dN/dy ~ 1000 rapidity density of produced hadrons. 
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Figure 4.1: The suppression factor, Raa(p±), for non-photonic electrons from the 
decays of quenched heavy quark (c+b) jets is compared to Phenix [237, 238] and 
preliminary STAR data [239, 240] data in central Au+Au reactions at 200 AGeV. 
Shaded bars indicate systematic errors, while thin error lines indicate statistical ones. 
All calculations assume initial dN g /dy = 1000. The upper yellow band from [374] 
takes into account radiative energy loss only, using a fixed L = 6 fm; the lower 
yellow band is our new prediction, including both elastic and inelastic losses as well 
as jet path length fluctuations. The bands provide an estimate of current theoretical 
uncertainties. The dashed curve shows the electron suppression using inelastic and 
TG elastic loss with bottom quark jets neglected. 

The observed "perfect fluidity" [77, 85-87, 375] of the sQGP at long wave- 
lengths (pt < 2 GeV/c) provides direct evidence for highly nonperturbative 
dynamics leading to rapid equilibration and PdV work consistent with the 
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QCD equation of state [104, 376]. Above p T > 5 - 7 GeV [77, 85-87, 377], 
pQCD appears to provide increasingly reliable predictions, at least for the nu- 
clear modification of light parton jets [274, 275, 277, 297, 371-373]. So the 
question raised by the electron data is: to how short wavelengths, 1/pt, is 
the novel nonperturbative physics dominant in the strongly interacting Quark 
Gluon Plasma (sQGP) [104]? This question is important because hard jets can 
be utilized as effective "external" tomographic probes of the bulk sQGP matter 
only if their dynamics can be predicted reliably. Otherwise, jet quenching can 
only be an additional signal of new physics instead of providing a calibrated 
diagnostic probe of that physics. 

The upper band of Fig. 4.1 shows that the predictions from [374] consid- 
erably underestimated the electron nuclear modification even out to px ~ 8 
GeV. This discrepancy points to either (1) missing perturbative QCD physics, 
(2) incomplete understanding of the initial heavy quark production and/or (3) 
novel non-perturbative mechanisms affecting partonic physics out to pt > 10 
GeV. We note that pt ~ 8 GeV (single non-photonic) electrons originate in 
our calculations from the fragmentation and decay of both charm and bottom 
quarks with transverse momenta px ~ 12 ± 4 GeV (see Fig. 4.5 in [374]). 

Possibility (3) is of course the most radical and would imply the persistence 
of non-perturbative physics in the sQGP down to extremely short wavelengths, 
i.e. 1/pt ~ 1/50 fm. These processes can be postulated to improve the fit to 
the data[378], but at the price of losing theoretical control of the tomographic 
information from jet quenching data. DGVW [374] showed that by arbitrarily 
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increasing the initial sQGP densities to unphysical dN/dy^ 4000, the non- 
photonic electrons from heavy quarks can be artificially suppressed to Raa ~ 
0.5. Thus conventional radiative energy loss requires either violation of bulk 
entropy bounds or nonperturbatively large a s extrapolations of the theory. 
Even when ignoring the bottom contribution, Ref. [379], found that a similarly 
excessive transport coefficient [380], q e ff ~ 14 GeV 2 /fm, was necessary to 
approach the measured suppression. 

The main theoretical problem is that bottom quark jets are too weakly 
quenched by radiative energy loss. Their contribution significantly reduces 
the single electron suppression [374] compared to that of the charm jets alone. 
While the current data cannot rule out the possibility (2) that the bottom 
production in this kinematic range is overestimated in NLO theory [381], in this 
paper we pursue the more conservative approach of examining the inclusion of 
a previously neglected pQCD component of the physics. 

The recent work by Mustafa [236, 382] and others [383] motivated us to 
revisit the assumption that pQCD elastic energy loss [147] is negligible com- 
pared to radiative. In earlier studies, the elastic energy loss [147, 224, 229, 
230, 234, 384, 385] was found to be dE el /dx ~ 0.3 - 0.5 GeV/fm, which was 
erroneously considered to be small compared to the several GeV/fm expected 
from radiative energy loss. The apparent weakness of conventional pQCD col- 
lisional energy loss mechanisms was also supported by parton transport theory 
results [138, 249], which showed that the typical thermal pQCD elastic cross 
section, a e \ ~ 3mb, is too small to explain the differential elliptic flow at high 
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Pt > 2 GeV and also underestimates the high p T quenching of pions. 
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Figure 4.2: AE/E for u,c,b quarks as a function of E. A Bjorken expanding 
QGP with path length L = 5 fm and initial density fixed by dN g /dy = 1000 is 
assumed. The curves are computed with the coupling as = 0.3 held fixed. For 
Dehye mass {id oc (dN g / dy)^ 1 /^, the gluon mass is hd/\/2, the light quark mass 
is /Ad/2, the charm mass is 1.2 GeV, and the bottom mass is 4.75 GeV. Radiative 
DGLV first order energy loss is compared to elastic parton energy loss (in TG or 
BT approximations). The yellow bands provide an indication of current theoretical 
uncertainties in the elastic energy loss for bottom quarks. 

In contrast, Mustafa found that radiative and elastic energy losses for heavy 
quarks were in fact comparable over a very wide kinematic range accessible at 
Rhic. In Fig. 4.2, we confirm Mustafa's finding [236, 382] and extend it to 
the light quark sector as well; the fractional energy loss, AE/E, from DGLV 
radiative for u, c, b quarks (solid curves) is compared to Thoma-Gyulassy (TG) 
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[224] and Braaten-Thoma (BT) [229, 230] estimates of elastic (dashed curves). 
For light quarks, the elastic energy loss decreases more rapidly with energy 
than radiative energy loss, but even at 20 GeV the elastic is only 50% smaller 
than the radiative. 

From Fig. 4.2, we see that above E > 10 GeV, the light and charm quarks 
have similar elastic and inelastic energy losses. But due to the large mass 
effect [285, 298, 305, 379, 386, 387], both radiative and elastic energy losses 
remain significantly smaller for bottom quarks than for light and charm quarks. 
We present both TG and BT as a measure of the theoretical uncertainties of 
the Coulomb log. These are largest for the heaviest b quark, for which the 
leading log approximation [224, 229, 230] breaks down in the kinematic range 
accessible at Rhic as they are not ultrarelativistic. With advanced numerical 
covariant transport techniques [138] the theoretical errors on the elastic energy 
loss effects can be reduced considerably in the future. 



4.2 Theoretical Framework 

The quenched spectra of partons, hadrons, and leptons are calculated as in 
[374] from the generic pQCD convolution 

Ed 3 a(e) E^oiQ) , . . , , N 

= ^ 3 ® E f) ® D (Q ^ H 0) ® /( Q ^ e )- (4- 1 ) 

where Q denotes quarks and gluons. For charm and bottom, the initial quark 
spectrum, Ed 3 a(Q)/dp 3 , is computed at next-to-leading order using the code 
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from [381, 388]; for gluons and light quarks, the initial distributions are com- 
puted at leading order as in [277]. P(Ei — > Ef) is the energy loss probability, 
D(Q — > Hq) is the fragmentation function of quark Q to hadron Hq, and 
f{Hq e) is the decay function of hadron Hq into the observed single elec- 
tron. We use the same mass and factorization scales as in [389] and employ 
the CTEQ5M parton densities [161, 390] with no intrinsic k?- As in [389] we 
neglect shadowing of the nuclear parton distribution in this application. 

We assume that the final quenched energy, Ef, is large enough that the 
Eikonal approximation can be employed. We also assume that in Au+Au 
collisions, the jet fragmentation function into hadrons is the same as in e + e~ 
collisions. This is expected to be valid in the deconfined pQCD medium case, 
where hadronization of Q — > Hq cannot occur until the quark emerges from 
the deconfined sQGP medium. 

The main difference between our previous calculation [374] and the present 
one is the inclusion of two new physics components in the energy loss proba- 
bility P(Ei — > Ef). First, P(Ei — > Ef) is generalized to include for the first 
time both elastic and inelastic energy loss and their fluctuations. We note that 
Vitev [391-393] was the first to generalize GLV to include initial state elastic 
energy loss effects in D+Au. In this work, Eq. (4.2) extends the formalism to 
include final state elastic loss effects in A + A. 

The second major change relative to our previous applications is that we 
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now take into account geometric path length fluctuations as follows: 



where L is the locally determined effective path length of the jet given its 
initial production point x± and its initial azimuthal direction relative to 
the impact parameter plane (x, y). The geometric path averaging used here is 
similar to that used in [359] and by Eskola et al. [306]. However, the inclusion 
of elastic losses together with path fluctuations in more realistic geometries 
was not considered before. 

We consider a diffuse Woods-Saxon nuclear density profile [394], which 
creates a participant transverse density, p p (x±), computed using the Glauber 
profiles, Ta{x), with inelastic cross section gnn — 

42 mb. The bulk sQGP 

transverse density is assumed to be proportional to this participant density, 
and its form is shown (for the y = slice) in Fig. 4.4 by the curve labeled 
Pqgp- The distribution of initial hard jet production points, x±, is assumed 
on the other hand to be proportional to the binary collision density, Taa = 
Ta{x + b/2)TA{x — 6/2). This is illustrated in Fig. 4.4 by the narrower curve 
labeled p Jct . 

The combination of fluctuating DGLV radiative [278] with the new elas- 
tic energy losses and fluctuating path lengths (via the extra d 2 x±d(() integra- 




®P d (E^E f ;L(x ± , (/)))] 



(4.2) 




(4.3) 
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Figure 4.3: Distribution of path lengths (given by Eq. (4.3)) traversed by 
hard scatterers; the lengths, L(x±,(p), are weighted by the probability of pro- 
duction and averaged over azimuth. An equivalent formulation of Eq. (4.4) is 
R I Q = f dLl/N bin dN Un /dL f de(l - e) n P^e; L) . Since this is a purely geometric 
quantity, it is the same for all jet varieties. Also displayed are the single, represen- 
tative pathlengths, Lq, used as input in approach II. Note the hierarchy of scales 
with glue requiring the shortest, then charm, light quarks, and bottom the longest 
effective pathlength. 

tions) adds a high computational cost to the extended theory specified by 
Eqs. (4.1,4.2). In this first study with the extended theory, we explore the 
relative order of magnitude of the competing effects by combining two simpler 
approaches. 

In approach I, we parameterize the heavy quark pQCD spectra by a simpler 
power law, Ed 3 aQ/d 3 k oc l/p™ +2 , with a slowly varying logarithmic index 
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n = n(px)- F° r the pure power law case, the partonic modification factor, 
Rq = d<TQ Ual I da % Q tial , (prior to fragmentation) is greatly simplified. This 
enables us to perform the path length fluctuations numerically via 

Rl Q = J^J 4^Jb) TAA{ *^ ] j d < l -'Y lpI Q^ L ^^^ ( 4 - 4 ) 
where 

P^e-L) = J dxP Q ]L)P Qt ei(e - x;L). (4.5) 

Both the mean and width of those fractional energy losses depend on the local 
path length. We emphasize that no externally specified a priori path length, 
L, appears in Eq. (4.4); the path lengths explore the whole geometry. Fig. 4.3 
shows the broad distribution of lengths traversed by hard partons in approach 
I. 

In the second approach, we determine effective path lengths, Lq, for each 
parton flavor, Q, by varying fixed L predictions Rg(p±, L) and comparing to 
approach I; see Fig. 4.3. In approach II, Rg(p±, L) is calculated directly from 
Eq. (1) with P(Ei — > Ef) in Eq. (2) replaced by the fixed L approximation 

P{E t -> E f , L) w P rad (Ei E;L) <S> P d (E -> E f - L). (4.6) 

Here, jet quenching is performed via two independent branching processes in 
contrast to the additive convolution in Eq. (4.4). For small energy losses the 
two approaches are similar. They differ however in the case of long path lengths 
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Figure 4.4: Transverse coordinate (x, 0) distribution of surviving px = 15 GeV, 
Q = g,u,c,b jets moving in direction <p = as indicated by the arrows. Units are 
arbitrary for illustration. The transverse (binary collision) distribution of initial jet 
production points, pj e t(x, 0), is shown at midrapidity for Au+Au collisions at b = 2.1 
fm. The ratio pq / pj e t (see Eq. (4.8) ) gives the local quenching factor including elastic 
and inelastic energy loss though the bulk QGP matter distributed as pqgp(x, 0). 

when the probability of complete stopping approaches unity. In the convolu- 
tion method, the probability of e > 1 is interpreted as complete stopping, 
whereas in the branching algorithm the long path length case is just highly 
suppressed. In both cases we take into account the small finite probability 
that the fractional energy loss e < due to fluctuations. 

To illustrate the difference in approach II, consider the case of power law 
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initial Q distributions as in Eq. (4.4). In this case 

J2£(pjL. Lq) = ((1 - e r Q (L Q )r(l - e e Q (L Q )r) AE . (4.7) 

The branching implementation is seen via the product of two distinct factors 
in contrast to the one quenching factor in Eq. (4.4). For small (eg 6 ) both 
approaches obviously give rise to the same Rq = 1 — n(eQ). 

In both approaches, fluctuations of the radiative energy loss due to gluon 
number fluctuations are computed as discussed in detail in Ref. [374, 387]. 
This involves using the DGLV generalization [278] of the GLV opacity ex- 
pansion [274] to heavy quarks. Bjorken longitudinal expansion is taken into 
account by evaluating the bulk density at an average time t = L/2 [374, 387]. 
For elastic energy loss fluctuations, the full Fokker-Planck solution applica- 
ble for small fractional energy loss is approximated by a Gaussian centered 
at the average loss with variance a 2 el = 2T(AE el (pT, L)) [249]. In approach I 
the correct, numerically intensive integration through the Bjorken expanding 
medium provides AE (pt, L). In approach II the r = L/2 approximation is 
again used; numerical comparisons show that for L ~ 2 — 7 fm this reproduces 
the full calculation well. Finally, we note that we use the additional numerical 
simplification of keeping the strong coupling constant as fixed at 0.3. 
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4.3 Numerical Results: Parton Level 

In Fig. 4.5, we show the quenching pattern of Q from the second approach 
for a "typical" path length scale L = 5 fm, similar to that used in previous 
calculations [374]. The curves show Rq(pt), prior to hadronization, for Q = 
g, u, c, b. The dashed curves show the quenching arising from only the DGLV 
radiative energy loss. The solid curves show the full results after including 
TG elastic as well as DGLV radiative energy loss. Adding elastic loss is seen 
to increase the quenching of all flavors for fixed path length. Note especially 
the strong increase of the gluon suppression and the factor ~ 2 increase of the 
bottom suppression. The curious switch of the u and the c quenching reflects 
the extra valence (smaller index n u ) contribution to light quarks. 

Fig. 4.5 emphasizes the unavoidable result of using a fixed, "typical" path 
length scale, L, in jet tomography: the pion and single electron quenching can 
never be similar. If pions were composed only of light quarks and electrons 
only of charm, then we would expect comparable quenching for both. However, 
contributions from highly quenched gluons decrease the pion Raa while weakly 
quenched bottom quarks increase the electron Raa- Therefore, in the fixed 
length scenario, we expect a noticeable difference between pion and single 
electron suppression patterns. 

In Fig. 4.4 we use fluctuating path lengths. The solid curves labeled by 
the parton flavor Q = g,u,c,b show the relative transverse coordinate density 
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Figure 4.5: Partonic nuclear modification, Rq(pt) via Eq.(4.7), for g,u,c,b as 
a function of pt for fixed L=5 fm path length and dN g /dy = 1000. Dashed curves 
include only radiative energy loss, while solid curves include elastic energy loss as 
well. 

of surviving jets defined by 

p Q {x, 0) = p Jct (f) J de(l - e) n P^{e; L{x, 0)). (4.8) 

Pq is given by the initial transverse x production distribution times the quench- 
ing factor from that position in direction with final momentum p?. The case 
shown is for a pr = 15 GeV jet produced initially at (x, 0) and moving in the 
direction = along the positive x axis. The quenching is determined by 
the participant bulk matter along its path pqgp(x + vt,0,t), and varies with 
x because the local path length L = L(x, 0, 0) changes according to Eq. (4.3). 
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What is most striking in Fig. 4.4 is the hierarchy of Q-dependent length 
scales. No single, representative path length can account for the distribution 
of all flavors. In general heavier flavors are less biased toward the surface 
(in direction <fi) than lighter flavors since the energy loss decreases with the 
parton mass. Gluons are more surface biased than light quarks due to their 
color Casimir enhanced energy loss. In addition, note the surprising reversal of 
the u and c suppressions, also seen in Fig. 4.5. Fig. 4.2 shows that the energy 
loss of c is somewhat less than for u; however, the higher p? power index, n, 
of c relative to u - as predicted by pQCD and due to the valence component 
of u - compensates by amplifying its quenching. 

However, it is clear that none of the distributions can be usefully cate- 
gorized as surface emission. The characteristic widths of these distributions 
range from Ax ~ 3 — 6 fm. Such a large dynamic range of path length fluc- 
tuations allows for the consistency of simultaneously reproducing both the 
electron and pion data. 

We turn next to Figs. 4.6 and 4.7 to show the interplay between dynamical 
geometry seen in Fig. 4.4 and the elastic-enhanced quenching of partons. In 
Figs. 4.6 and 4.7 the solid green curves labeled "DGLV+TG/BT: Full Geom- 
etry" are the results using approach I based on Eq. (4.4). The curves labeled 
TG and BT are from approach II based on Eq. (4.7). The effective fixed Lq 
in II were taken to match approximately the green curves in which full path 
length fluctuations are taken into account. This procedure is not exact be- 
cause of the different numerical approximations involved, but the trends are 
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well reproduced. The Lq are determined only to ~ 0.5 fm accuracy, as this 
suffices for our purposes here. We show the comparison between approaches 
I and II for heavy quarks in Fig. 4.6 using L c = 4.5 and Lb = 6.5 fm, and 
for gluons and light quarks in Fig. 4.7 using L g = 4.0 and L u = 5.0 fm; see 
Fig. 4.3 for a visual comparison of the input length distributions used. This 
hierarchy of Q-dependent length scales is in accord with that expected from 
Fig. 4.4. 
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Figure 4.6: Heavy quark jet quenching before fragmentation into mesons for 
dNg/dy = 1000. Solid green curves show the results of approach I based on Eq. (4.4) 
including full geometric path length fluctuations and DGLV radiative and TG elastic 
energy loss for c and b quarks. Upper and lower yellow bands show predictions using 
approach II via Eq. (4.7) with effective path lengths taken as L c = 4.5 and = 6.5 
fm. As previously noted in Fig. 4.2, the difference between TG and BT curves 
indicates the magnitude of theoretical uncertainties in the elastic energy loss. 
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Note that geometric fluctuations reduce charm quark quenching in Fig. 
4.6 relative to the fixed L = 5 fm case in Fig. 4.5. This is seen even more 
dramatically for gluons in Fig. 4.7: the overquenching of gluons in Fig. 4.5 is 
reduced by a factor ~ 2 when geometric fluctuations are taken into account. 
Nevertheless, even with path fluctuations the gluons are still quenched by a 
factor of 10 when elastic energy loss is included in addition to radiative. 
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Figure 4.7: As in Fig. 4.6 but for light u, d quarks and gluons. The yellow bands 
are computed in this case with effective g, u path lengths L g = 4.0 and L u = 5.0 fm 
based on Eq. (4.7). Note that charm and light quark quenching are similar in this 
Pt range. 

The amplified role of elastic energy loss is due to its smaller width for 
fluctuations relative to radiative fluctuations. Even in moderately opaque me- 
dia with L/X ~ 10, inelastic energy loss fluctuations are large because only 
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a few, 2-3, extra gluons are radiated [275]. Thus, gluon number fluctuations, 
0(l/y/N g ) lead to substantial reduction in the effect of radiative energy loss. 
On the other hand, elastic energy loss fluctuations are controlled by collision 
number fluctuations, 0(a/A/L), which are small in comparison. Therefore, 
fluctuations of the elastic energy loss do not dilute the suppression of the nu- 
clear modification factor as much as N g fluctuations. The increase in the sensi- 
tivity of the final quenching level to the opacity is a novel and useful byproduct 
of including the elastic channel. The inclusion of elastic energy loss significantly 
reduces the fragility of pure radiative quenching [306] and therefore increases 
the sensitivity of jet quenching to the opacity of the bulk medium. 

4.4 Numerical Results: Pions and Electrons 

We now return to Fig. 4. 1 to discuss the consequence of including elastic energy 
loss of c and b on the electron spectrum. The largest source of uncertainty, 
represented by the lower yellow band, is the modest but poorly determined 
elastic energy loss, AE/E « 0.0 — 0.1, of bottom quarks (see Fig. 4.2). Nev- 
ertheless, it is remarkable that the lower yellow band can reach Raa ~ 0.5 
in spite of keeping dN g /dy = 1000, consistent with measured multiplicity, 
and using a conservative as = 0.3. While the preliminary data suggest that 
R e AA < 0.4, the Rad+Elastic band is not inconsistent with the preliminary 
data above pt~ 7 GeV within present experimental and theoretical errors. 
Our tentative conclusion, therefore, is that the combined elastic and inelastic 
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pQCD approach may solve a substantial part of the heavy quark quenching 
puzzle. 
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Figure 4.8: The consistency of the extended jet quenching theory is tested by 
comparing its prediction to the nuclear modification of the ir° spectra observed by 
Phenix [368-370]. 

However, as emphasized in [374] , any proposed solution of this puzzle must 
also be checked for consistency with the extensive pion quenching data [368- 
370], for which preliminary data now extend out to pr ~ 20 GeV. That this is 
a challenge is clearly seen in Fig. 4.5, where for fixed L = 5 fm, the addition of 
elastic loss would significantly overpredict the quenching of pions. However, 
the simultaneous inclusion of path fluctuations leads to a decrease of the mean 
g and u, d path lengths that partially offsets the increased energy loss. There- 
fore, the combined three effects considered here makes it possible to satisfy 
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R AA < 0.5 without violating the bulk dN g /dy = 1000 entropy constraint and 
without violating the pion quenching constraint R A ° A ~ 0.2 ±0.1 now observed 
out to 20 GeV; see Fig. 4.8. We note that the slow rise of R AA with px in 
the present calculation is due in part to the neglect of initial kj* smearing that 
raises the low px region and the EMC effect that lowers the high px region 
(see [277]). 

4.5 Sensitivity and Imprecision in WHDG 

The ultimate goal of hard probe physics is jet tomography, the use of the 
high-pr hadronic attenuation pattern to learn about the medium density. One 
needs then a theoretical description of the processes involved that describes 
the data, has a mapping from its parameter(s) to the matter density, and is 
under calibrated control. 

For these tools to be truly tomographic, however, there is the additional re- 
quirement that the resulting medium measurement be reasonably precise. The 
overall precision of density determination is set by a combination of experimen- 
tal precision and theoretical sensitivity: if the experimental data are infinitely 
precise, any model with the previous paragraph's characteristics could asses 
the medium density with similar precision; but experimental error exists and 
must be taken seriously, limiting the precision of medium determination by 
the innate sensitivity of the theoretical model to changes in its input parame- 
ter. Additionally there will also be theoretical error. Sources for this include, 
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but are not limited to: neglect of higher order terms; treatment of the run- 
ning coupling, especially at low q 2 ; and the translation of the realistic medium 
density to the brick problem for which the analytic theoretical derivation was 
applied. 

[306] claimed that BDMPS-based radiative energy loss models are a fragile 
probe of the medium; i.e. large changes in q, the input parameter for the theory, 
are not reflected in Raa predictions incompatible with data. Surface emission, 
the process in which only the partons produced very close to the medium edge 
escape for observation, was posited as the reason for their observed insensitivity 
to increased medium opacity. This is an appealing argument, and [296, 395] 
provided some evidence for this idea. The very pessimistic conclusion, then, 
was that jet tomography is not possible at Rhic and will not be possible at 
the Lhc. 

We counter the dire claims in [306] by first noting that several approxi- 
mations were made in their description of the medium produced in heavy ion 
collisions. In [306], the authors used a hard cylinder initial nuclear density and 
did not include Bjorken expansion. [296, 395] used the realistic Woods-Saxon 
nuclear density distribution but again neglected Bjorken expansion. These 
certainly aided in alleviating some of the high numerical cost involved in the 
calculations. However we feel these geometrical simplifications strongly biased 
their results to the surface of the medium, resulting in a significant loss in 
model sensitivity. Moreover, experimental error bars have shrunk consider- 
ably in the time since [306] were published. These issues raise serious doubts 
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about the validity of their conclusions. 
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Figure 4.9: Model sensitivity to varying dN g /dy, which is much greater than 
observed in [306]. The susceptibility to changes in the fixed ot s suggests the need to 
investigate the theoretical error stemming from the coupling. 

One can readily see in Fig. 4.9 that the use of fluctuating path lengths 
in a realistic and expanding background medium geometry results in a pre- 
dicted Raa{pt) inconsistent with data when the medium density is artificially 
increased by a factor of 2. Sensitivity to medium density changes is enhanced 
over pure radiative loss models by the inclusion of elastic energy loss, due 
to its smaller width for fluctuations relative to radiative Fig. 4.10. We note, 
however, that models using GLV-type radiative loss or higher-twist loss only 
are similarly not fragile [356]. More recent work in which jets were propa- 
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gated through a medium evolved through hydrodynamics found that energy 
loss using the ASW quenching weights in conjunction with up-to-date exper- 
imental results are in fact not fragile, either [396]. The pion Raa{pt) then 
has the potential to be a quite good experimental signal for jet tomographic 
studies at Rhic. However the large systematic theoretical error seen in Fig. 
4.9 associated only with changes in the coupling raise serious doubts about the 
theoretical imprecision of calculated results. Before any firm conclusions may 
be made regarding the possibility of tomography these issues must be more 
rigorously investigated. 

4.6 LHC Predictions 

We show in Fig. 4.12 the large qualitative difference in dependence for 
predicted Lhc pion Raa from the WHDG model and from two different im- 
plementations of the ASW model; one easily sees the dramatic rise in Raa 
with increasing from the first as opposed to the flat in p? results of the 
latter two. The consistency of Vitev's curve [397] with ours over a range of 
dNg/dy and the consistency of the two ASW calculations [296, 306] suggest 
that this will be a robust test of the energy loss mechanisms. The origin of 
the difference can be easily understood given our discussion of Rhic results. 
For the case of WHDG, the pion Raa at Rhic does not require much over- 
absorption. The modest increase in medium density, ~ 2 — 3 based on either 
an extrapolation of Phobos results [398, 399] or predictions from the CGC 
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Figure 4.10: The ratio of rms width, <j(L), to the mean fractional energy loss (e) 
for radiative, elastic and convoluted energy loss distributions is shown as a function 
of the path length, L , for a Bjorken expanding plasma with dN g jdy = 1 000 . The case 
of an up quark jet with E = 15 GeV is shown. Notice that the elastic distribution 
is signihcantly narrower than the radiative one. This amplifies the effect of elastic 
energy loss on Raa relative to radiative. 

[400, 401], for the Lhc leads to small energy losses at high momenta that can 
be well approximated by the pocket asymptotic energy loss formulae 

e rad = AE rad /E ~ a 3 s log{E/fi 2 L)/E (4.9) 
e el = AE el /E ~ a? t log(VET/m g )/E; (4.10) 

see Fig. 4.11. As pt increases the \og(E)/E reduction in energy loss is not 
compensated by the slow (in comparison to Rhic) increase in the power law 




Chapter 4: Fluctuations in Jet Tomography 



114 



partonic production spectrum, dN/dpx oc l/p^ PT ^ (see Appendix B); thus 
Raa increases with p?. On the other hand, the ASW models mimic the small 
normalization of the Rhic data by highly suppressing their jets; as discussed 
in the subsequent paragraphs, the significant quenching leads to a loss of in- 
formation on the details of the energy loss process, flattening the results. The 
two ASW models represented in Fig. 4.12 used EKRT-type medium density 
scaling [402] making the Lhc medium ~ 7 times more dense than at Rhic. 
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Figure 4.11: Radiative DGLV and TG and BT elastic partonic fractional energy 
loss as implemented in WHDG [152] at Lhc momenta for all jet species at fixed 
L = 5 fm and (a) dN g /dy = 1750 and (b) dN g /dy = 2900; the former density comes 
from the Phobos extrapolation [87, 403], the latter from the KLN model of the 
color glass condensate (CGC) [400]. AE/E exhibits asymptotic behavior for both 
energy loss channels as given by the analytic pocket formulae Eq. (4.9) and Eq. 
(4.10). 
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The important quantity in any of these calculations is P(e), the probability 
that the final momentum is some fraction of the initial momentum, = (1 — 
e)p l T . For calculations that include only radiative processes, P(e) = P ra d(^), 
where 




P°5(x) + P rad (x) 

e- Ng E ^PradM = e- N ^S(x) + e- N °N g P rad:1 (x) + .... 



(4.11) 

The possibility that no gluons are emitted (subsequent to the original gluon 
radiation created by the initial hard scatter) is encapsulated in the coefficient, 
Pq , of the S(x) term. 

Due to the approximations used in the radiative calculations overabsorp- 
tion, P{e > 1), has a large support for highly suppressed jets. Usually, one 
of two prescriptions is applied to remove this unphysical artifact. Either the 
integrated excess probability weighs an explicit delta function, 

/oo 
dxP old (x)S(l-e), (4.12) 

or reweighs (rw) the original distribution, 

P rw {e) = f0O , \ , , Poid(e)0(l - e). (4.13) 
J 1 dxP i d {x) 

Clearly the latter approach leads to larger Raa values for the two. For large 
overabsorption energy loss details are lost to the removal of the unphysical 



Chapter 4: Fluctuations in Jet Tomography 



116 



e > 1 region. 
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Figure 4.12: Lhc predictions of pion Raa(pt) for several energy loss models. 
Curves correspond to WHDG [152], GLV [397], PQM [296], and ASW [306]. For 
the latter two, rw indicates the use of Eq. (4.13) for overabsorption; otherwise Eq. 
(4.12) was employed. Also for the latter two, q = 100 and q = 68 were the values of 
the ASW input parameter, respectively. Notice the sharp rise in pt for the WHDG 
and GLV curves as opposed to the flatness of the ASW and PQM ASW results. For 
heavy quark Raa(pt) predictions from WHDG see Chapter 5. 

4.7 Conclusions 

While the results presented in this paper are encouraging, further improve- 
ments of the jet quenching theory will be required before stronger conclusions 
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can be drawn. More work is needed to sort out coherence and correlation 
effects between elastic and inelastic processes that occur in a finite time and 
with multiple collisions. Classical electrodynamics calculations presented in 
[244] suggested that radiative and elastic processes could destructively inter- 
fere over lengths far longer than previously thought. However, the formalism 
used does not disentangle known radiative Ter-Mikayelin effects [305] from 
expected finite size collisional effects. More important, the authors inconsis- 
tently treat their current as conserved (using j° = kji/u) while simultaneously 
dropping the backwards-moving partner parton term from the current (thus 
violating current conservation by producing, at a finite time, a charge without 
a corresponding opposite charge); this results in a spurious subtraction of the 
(negative) binding energy of the pair and a correspondingly drastic reduction 
in the elastic energy loss. As seen in [245], a proper accounting of the current 
shows finite size effects persist out only to the expected lengths of order the 
screening scale, 1///d~ 1 f m - Recent work on the quantum mechanical treat- 
ment of elastic loss in a finite-sized medium [247, 404] also concluded that 
these effects are small. At the same time, the authors of [247, 404] disagree on 
the effects of interference between the radiative and elastic contributions; a so- 
lution of this problem may involve the proper inclusion of transition radiation 
[302]. 

In addition, the radiative and elastic energy losses depend sensitively on 
the coupling, AE rad oc a| and AE el cx a|. Future calculations will have 
to relax the fixed as approximation and allow it to run. The energy loss 
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involves integrals that probe momentum scales that are certainly nonpertur- 
bative. Therefore it will be important to study the irreducible uncertainty as- 
sociated with the different maximum 0:5 cutoff prescriptions commonly used. 
Further, improved numerical techniques will have to be devised to bypass 
the approximations employed here to incorporate essential path length fluc- 
tuations. Nevertheless qualitatively different high momentum behaviors are 
predicted for 7r° suppression at the Lhc from the BDMPS-based ASW models 
and the GLV radiative and DGLV radiative convolved with elastic loss models: 
for the former Raa{pt) is approximately flat in p T ; for the latter it increases 
significantly with p?. From an experimental perspective, direct measurement 
of D spectra would be valuable to separate the different bottom and charm 
jet quark dynamics that is at present convoluted in the electron spectra. 
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Chapter 5 

pQCD vs. AdS/CFT Tested by 
Heavy Quark Energy Loss 

5.1 Introduction 

Recent discoveries at Rhic [77, 104, 375, 405] have led to suggestions [335, 406- 
414] that the properties of strongly coupled quark gluon plasmas (sQGP) 
produced in ultra-relativistic nuclear collisions may be better approximated 
by string theoretic models inspired by the AdS/CFT gravity-gauge theory 
correspondence [312, 415, 416] than conventional Standard Model perturba- 
tive QCD (pQCD). Four main classes of observables have attacted the most 
attention: (1) Entropy production as probed by multiplicity distributions 
[398], (2) "Perfect" Fluidity [405] as probed by collective elliptic flow ob- 
servables [106, 121, 138, 417-421], (3) Jet Quenching and Tomography as 
probed by high-p-r hadrons [192, 277, 297, 373, 422] and nonphotonic lep- 
tons [237, 423, 424], and (4) Dijet-Bulk Correlations as probed by two and 
three particle correlations [425-431]. 
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Figure 5.1: Nonphotonic electron predictions from the AdS/CFT drag and pQCD 
WHDG [152] models compared to Rhic data from Phenix [237, 238] and Star 
[239, 240] in central 200 GeV Au + Au collisions. For a density ofdN g /dy = 1000 
the perturbative estimate is in qualitative disagreement; the drag calculation is in 
agreement in its range of reasonable choices of input parameters. See text for details 
of the drag calculation; the electron fragmentation was done as in [152] 

Qualitative successes of recent AdS/CFT applications [335, 406-414] to 
nuclear collision phenomenology include the analytic account for (1) the sur- 
prisingly small (~ 3/4) drop [332, 432] of the entropy density in lattice QCD 
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calculations relative to Stefan-Boltzmann, (2) the order of magnitude reduc- 
tion of the viscosity to entropy ratio r]/s predicted relative to pQCD needed to 
explain the seemingly near perfect fluid flow of the sQGP observed at Rhic, 
(3) the unexpectedly large stopping power of high transverse momenta heavy 
quarks as inferred from heavy quark jet quenching, see Fig. 5.1, and elliptic 
flow, and (4) the possible occurrence of conical "Mach" wave-like correlations 
of hadrons associated with tagged jets. 

While quantitative and systematic comparisons of AdS/CFT gravity dual 
models with nuclear collision data are still incomplete and while the conjec- 
tured double Type IIB string theoryOconformal Supersymmetric Yang-Mills 
(SYM) gauge theory<B-non-conformal, non-supersymmetric QCD correspon- 
dence remains under debate (see, e.g. [433]), the current successes provide 
strong motivation to seek more sensitive experimental tests that could help 
guide the theoretical development of novel theoretical approaches that may be 
needed to explain recent Rhic and soon Lhc heavy ion data. 

The aim of this chapter is twofold. First is to propose a robust observable 
that can more readily reveal the kinematic boundaries and reaction conditions 
where specific Standard Model weak coupling pQCD approximations may fail 
and where specific strong coupling AdS/CFT approximations may fail (or if 
they are applicable at all). Asymptotic freedom and the factorization theorems 
of QCD ensure that pQCD should apply above some hard scale Q(A, y/s) that 
depends in general on both atomic number, A, and center of mass energy, 
y/s. For jet production this scale is the greater of the gluon saturation scale 
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Ql ~ Ql \og(Ay/s / pt) , that defines a scale below which strong initial state 
nuclear modification of the parton distributions must be taken into account, 
and qL = J dzdap(z)q 2 oc Q 2 S log n (Q s ), below which strong final state nuclear 
modifications must be taken into account. For finite A < 238 and finite 
y/s < 10 ATeV the numerical values of these scales remains uncertain. Second 
is to suggest this observable will discriminate between current pQCD and 
AdS/CFT dynamical models. 

We start with the predicted nuclear modification, Rj[ a (pt), of the trans- 
verse momentum distribution of identified heavy quark jets produced in central 
Pb+Pb reactions at 5.5 ATeV at Lhc. Specifically, we propose that the double 
ratio of identified charm and bottom jet nuclear modification factors R%a/Raa 
is a remarkably robust observable that can distinguish between a rather wide 
class of pQCD energy loss mechanisms and a recently proposed class of grav- 
ity "drag" models [336, 342, 343, 347, 434] of heavy quark dynamics. Similar 
tests can be performed at Rhic when identified jet flavor detector upgrades 
are completed. The main advantage of Lhc in comparison to Rhic is of course 
the much higher p? kinematical range that will be accessible. The main advan- 
tage of Rhic is a better control of the initial state saturation physics because 
of extensive d + A and p + p control data at the same y/s. 

The current failure of pQCD based energy loss models [152, 298, 374, 435] to 
account quantitatively for the recent Rhic data from STAR [423] and Phenix 
[237, 424] on the nonphotonic electron spectrum provides additional moti- 
vation to focus on heavy quark jet observables. Unlike for light quark and 
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gluon jet observables, where pQCD predictions were found to be remark- 
ably quantitative [277, 297, 373], heavy quark jet quenching, especially as 
inferred indirectly for bottom quarks, appears to be significantly underpre- 
dicted [152, 298, 374, 435]. Current issues that cloud the pQCD based energy 
loss predictions are (1) the uncertainty in the initial state nuclear production of 
bottom to charm quarks, (2) the current controversy over the relative magni- 
tude of elastic versus radiative loss channels [152, 435], and (3) the possibility 
that short formation time nonperturbative hadronization effects may have to 
be taken into account [250, 436]. 

The AdS/CFT correspondence has so far been applied to heavy ion jet 
physics in three ways. The first involves the calculation of the QCD Wilson line 
correlator that corresponds to the radiative transport coefficient q [270, 345]. 
The second concentrates on estimating the heavy quark diffusion coefficient 
D [351] that is an input to a Langevin model of drag [249]. The third is 
a prediction of the heavy quark drag coefficient based on the gravity dual 
dynamics of a classical string in an AdS black brane background [342, 343, 347]. 
All three approaches remain under active debate (see, e.g. [347-349]). 

We focus in this chapter on the third proposed AdS / CFT application that 
involves the most direct string theoretic inspired gravity "realization" of heavy 
quark dynamics [336, 342, 343]. A heavy quark in the fundamental representa- 
tion is a bent Nambu-Goto string with one end attached to a probe brane and 
that trails back above the horizon of a D3 black brane representing the uniform 
strongly coupled SYM plasma heat bath. This geometry maps the drag force 
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problem into a modern string theoretic version of the 1696 Brachistochrone 
problem and yields a remarkable, simple analytic solution for the string shape 
and momentum loss per unit time. 



5.2 AdS/CFT compared to pQCD 

Exploiting the AdS/CFT correspondence, the drag coefficient for a massive 
quark moving through a strongly-coupled SYM plasma in the A = g 2 SY M^c 3> 
1, N c > 1, M Q > T* limit is given in [342, 343, 347] as 

dp T nV\(T*y 

where T* is the temperature of the SYM plasma as fixed by the Hawking 
temperature of the dual D3 black brane. This is parametrically quite similar 
to the Bethe-Heitler limit of incoherent radiative energy loss, in which case 



dp T T 3 

~dT ~ (5 ' 2) 



differing from Eq. (5.1) only by a factor of T/M, and quite dissimilar to the 
usual LPM-dominated asymptotic pQCD result, 

dPT rrr3 



—LT \og(p T /M). (5.3) 



dt 

There exists a maximum momentum, or 7 C w p^ lt /Mq beyond which Eq. 
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(5.1) cannot be applied. Self consistency within the classical string picture 
requires a time-like boundary for the string worldsheet [347]. For constant 
velocity this limits heavy quark "speeds" to 7 < j c , where 



The work of [343] relaxed the assumptions of infinite quark mass and con- 
stant velocity. The analytic form of Eq. (5.1) was found to well reproduce the 
full numerical results (most importantly \xq remained independent of pr) but 
with Mq no longer the bare quark rest mass. All the calculations in this paper 
were based upon the infinite bare mass approximation of [342, 347] with Mq 
replaced with realistic quark masses. 

We note that one may arrive at this same value through a different line 
of reasoning. Assume that the quark's constant velocity is maintained by an 
electric field. The largest electric field sustainable from the Born-Infeld action 
limits the magnitude of the momentum loss; the critical speed after which the 
field cannot be strong enough to keep the quark velocity constant is given by 



While in the infinitely strongly coupled plasma dual the quasiparticle pic- 
ture is not applicable, a similar "speed limit" arises for ordinary incoher- 
ent Bethe-Heitler (BH) radiative energy loss (See Appendix D.l). Landau- 
Pomeranchuck coherence effects invalidate the linear in E m p T rise of the 
energy loss when the formation time, te ~ E/M 2 , exceeds the mean free 




(5.4) 



Eq. (5.4) [352]. 
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path, A = I /per ~ l/cx s T (in a Debye-screened ultrarelativistic plasma). Just 
as in the form of loss itself the requirement that te < A therefore limits the 
applicability of BH to 7 < = Mq/o: s T, similar to the AdS/CFT speed 
limit, Eq. (5.4), but with one less power of T* /Mq. 

The speed limit Eq. (5.4) was found assuming a constant plasma temper- 
ature; this is not the case in experiment. Nevertheless, to get a sense of the 
Pt scale where the AdS/CFT approximation could break down, we will plot 
the momentum cutoffs from Eq. (5.4) for different SYM input parameters with 
two different assumptions for the mapping of QCD to AdS/CFT parameters, 
described in detail below. There is still additional ambiguity in the 7 C due 
to the time evolution of the QGP temperature. The smallest 7 C corresponds 
to the largest temperature; we take as the generous lower bound the extreme 
T(x = ; r = r ), shown as a '(' in the figures. For the largest 7 C we use 
T c , show as a "]." To further emphasize the possibility of corrections in the 
drag model we gradually fade the curves from the "(" to the "]." Surprisingly 
the introduction of a thermal plasma in AdS / CFT results in an effective mass 
smaller than the bare mass; this will result in a reduction of the momentum 
reach of the drag formalism. 

Applying Eq. (5.1) to heavy ion collisions requires an additional assump- 
tion about how to map the QCD temperature and coupling to the gedanken 
SYM world and its SUGRA dual. The "obvious" first prescription [437] takes 
9 sym = 9s constant, T* = T QCD , and N c = 3. However it was suggested in 
[437] that a more physical "alternative" might be to equate energy densities, 
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Figure 5.2: R c AA {px) and R aa (pt) predicted for central Pb+Pb at Lhc comparing 
AdS/CFT Eq. (5.1) and pQCD using the WHDG model [152] convolving elastic 
and inelastic parton energy loss. Possible initial gluon rapidity densities at Lhc are 
given by dN g /dy = 1750, from a Phobos [87, 403] extrapolation, or dN g /dy = 2900, 
from the KLN model of the color glass condensate (CGC) [400]. The top two curves 
from pQCD increase with px while the bottom two curves from AdS/CFT slowly 
decrease with p?- The AdS/CFT parameters here were found using the "obvious" 
prescription with q>sym = -05, tq = 1 fm/c, giving D = 3/2itT (abbreviated to D = 
3 in the figure). Similar trends were seen for the other input parameter possibilities 
discussed in the text. The "(" and "]" denote momenta after which possible string 
theoretic corrections may need to be considered; the curves' increasing transparency 
from "(" to "]" is meant to additionally emphasize this, see text. 

giving T* ~ T QCD /3 1 / 4 , and to fit the coupling A = 9 2 S ym^c ~ 5.5 in order to 
reproduce the static quark-antiquark forces calculated via lattice QCD. 

The string theoretic result for the diffusion coefficient used in the Langevin 



Chapter 5: pQCD vs. AdS/CFT 



128 



model is D = 2/\f\-nT* [351]. This illustrates well the problem of connecting 
the T* and A of SYM to "our" QCD world. Using the "obvious" prescription 
with a s = .3, N c = 3, one finds D ~ 1.2/2ttT. However, D = 3/2ttT was 
claimed in [351, 424] to fit Phenix data somewhat better. Note that D = 
3/27rT requires an unnaturally small a s ~ 0.05 that is very far from the 
assumed A ^> 1 't Hooft limit. 

We proceed by computing the charm and bottom nuclear modification 
factors, neglecting initial state shadowing or saturation effects. In order to 
correctly deconvolute such effects from the final state effects that we compute 
below, it will be necessary to measure nuclear modification factors in p + A as 
a function of (u,Pt) at Lhc just as d + A was the critical control experiment 
[77] at Rhic [104]. 

Final state suppression of high-p^ jets due to a fractional energy loss e, 
Pj, = (1 — e)p l T , can be computed knowing the Q-flavor dependent spectral 
indices nq+1 = — dl0 g PT log f ^^r j from pQCD or directly fromp+p — >■ Q+X 
data. The nuclear modification factor is then R® a (pt) =< (1 — e) nQ >, where 
the average is over the distribution P(e; Mq,p t , £) that depends in general on 
the quark mass, pr, and the path length I of the jet through the sQGP. As 
in [152] we find ng from FONLL production cross sections [381, 388], average 
over jets produced according to the binary distribution geometry, and compute 
I through a participant transverse density distribution taking into account the 
nuclear diffuseness. Given dN/dy of midrapidity partons, the temperature is 
computed assuming isentropic Bjorken ID Hubble flow (see Appendix C) with 
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Nf = 3. As emphasized in [152], detailed geometric path length averaging 
plays a crucial role in allowing consistency between 7r°, 77 and heavy quark 
quenching in pQCD. 

For AdS/CFT drag, Eq. (5.1) gives the average fractional energy loss as 
e = 1 — exp(— (XqI). Energy loss is assumed to start at thermalization, tq ~ 
0.6 — 1.0 fm/c, and stops when the confinement temperature, T c ~ 160 MeV, 
is reached. The exponentiated T 2 dependence in \iq leads to a significant 
sensitivity to the opacity of the medium, as well as to r and T c . 

To understand the generic qualitative features of our numerical results it 
is instructive to consider the simplest case of a geometric path average over a 
static, finite, uniform plasma of thickness L; then 

Raa(Pt) = r- ~ F , (5.5) 

where the pt dependence is carried entirely by the spectral index Uq(pt)- 

Two implementations of pQCD energy loss are used in this paper. The 
first is the full WHDG model convolving fluctuating elastic and inelastic loss 
with fluctuating path geometry [152]. The second restricts WHDG to include 
only radiative loss in order to facilitate comparison to [270]. Note that when 
realistic nuclear geometries with Bjorken expansion are used, the "fragility" of 
Raa for large q reported in [306] is absent in both implementations of WHDG. 

Unlike the AdS/CFT dynamics, pQCD predicts [152, 298, 374] that the 
average energy loss fraction in a static uniform plasma is approximately e ~ 
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Figure 5.3: The double ratio of R c aa (pt) to R aa (pt) predictions for Lhc using 
Eq. (5.1) for AdS/CFT and WHDG [152] for pQCD with a wide range of input 
parameters. The generic difference between the pQCD results tending to unity con- 
trasted to the much smaller and nearly px -independent results from AdS/CFT can 
be easily distinguished at Lhc. The "(" and "]" denote momenta after which pos- 
sible string theoretic corrections may need to be considered; the curves' increasing 
transparency from "(" to "]" is meant to additionally emphasize this, see text. 

nL 2 q\og{j>T I 'Mq) /px, with k a proportionality constant and q = fi 2 D /X g . The 
most important feature in pQCD relative to AdS/CFT is that € p qcd — > 
asymptotically at high-p^ while €Ads remains constant. uq(pt) is a slowly 
increasing function of momentum, Fig. 5.4; thus R P A ^ A CD increases with pt 
whereas R^f decreases. This generic difference can be observed in Fig. 5.2, 
which shows representative predictions from the full numerical calculations of 
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charm and bottom Raa{pt) at Lhc. 




Figure 5.4: The power law production index uq(pt) for Rhic and Lhc. While it 
is quite Bat-but slowly increasing-over a large momentum range at Lhc, at Rhic 
uq{pt) hardens appreciably as momentum increases. 



5.3 Double Ratio of charm to bottom R^ A 

A disadvantage of the Raa{pt) observable alone is that its normalization and 
slow pt dependence can be fit with different model assumptions compensated 
by using very different medium parameters. In particular, high value ex- 
trapolations of the q parameter proposed in [435] could simulate the flat pt- 
independent prediction from AdS/CFT. 

We propose to use the double ratio of charm to bottom Raa to amplify the 
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observable difference between the mass and pt dependencies of the AdS / CFT 
drag and pQCD-inspired energy loss models. One can see in Fig. 5.3 that 
not only are most overall normalization differences canceled, but also that the 
curves remarkably bunch to either AdS/CFT-like or pQCD-like generic results 
regardless of the input parameters used. 

The numerical value of R cb shown in Fig. 5.3 for AdS/CFT can be roughly 
understood analytically from Eq. (5.5) as, 

, M c n b (p T ) M c . , 

^w^rwr" 26 ' (5 - 6) 

where in this approximation all A, T*, L, and n c {p T ) ribipr) dependences 
drop out. 

The pQCD trend in Fig. 5.3 can be understood qualitatively from the 
expected behavior of € p qcd noted above giving (with n c « rib = n) 



RpQCD ~ 1 > (5-7) 

Pt 



where p c & = kti(pt)L 2 log(M&/M c )g sets the relevant momentum scale. Thus 
R cb — > 1 more slowly for higher opacity. One can see this behavior reflected 
in the full numerical results shown in Fig. 5.3 for moderate suppression, but 
that the extreme opacity q = 100 case deviates from Eq. (5.7). 

Supposing that Lhc data are similar to the pQCD predictions one might 
be able to distinguish between convolved elastic and inelastic loss, for which 
R cb monotonically increases, and purely radiative energy loss, for which R cb 
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dips to a minimum near p T ~ 10 GeV; one must be cautious, though, as 
many of the assumptions in the elastic energy loss derivations break down for 

p T > Mq. 

Future Rhic detector upgrades will allow for individual charm and bottom 
quark detection. Predictions for Raa(Pt) an d R aa (pt) at Rhic from Eq. (5.1) 
and pQCD are shown in Fig. 5.5. As seen in Fig. 5.4 the power law production 
index grows quickly at RHIC; we now expect the rapid increase in nq to over- 
come the (relatively in comparison) slow decrease in € p qcd so that, unlike at 
Lhc, dR AA /dpx ^ 0. In fact Fig. 5.5 shows that the full numerical results for 
R ( AA from pQCD and AdS/CFT drag both decrease with px- Nonetheless one 
may still examine the double ratio R cb , Fig. 5.6. While the larger index makes 
the grouping less dramatic at Rhic one may still differentiate between pQCD 
and AdS/CFT drag. Due to its smaller multiplicities, the temperature of the 
medium at Rhic is smaller than will be seen at Lhc; hence the AdS/CFT 
drag "speed limit", Eq. (5.4), is higher at Rhic than Lhc. 

5.4 Conclusions 

Possible strong coupling alternatives to pQCD in nuclear collisions were stud- 
ied based on a recent AdS/CFT model of charm and bottom energy loss. The 
predicted nuclear modification factors, R AA , were found to be decreasing as 
a function of pr, as compared to increasing as predicted from pQCD. We 
showed that the momentum dependence differences in the individual R AA can 
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Figure 5.5: R c aa (pt) and R aa (pt) for Rhic using AdS/CFT drag and the pQCD- 
based WHDG model. Note how the rapid increase in uq for Rhic as seen in Fig. 
5.4 overcomes the decrease in the fractional energy loss for the pQCD predictions 
for charm and bottom Raa{pt) so that both AdS/CFT drag and pQCD results 
decrease as a function of momentum at Rhic. 



be masked by taking extreme energy loss extrapolations to Lhc. However 
the double ratio R cb revealed very generic behavior, insensitive to the input 
parameters and radically different for the two coupling limits. Of crucial im- 
portance is the momentum range over which pQCD and AdS/CFT drag are 
self-consistent. Certainly pQCD must apply for oo, but the scale below 

which nonperturbative effects become important is not yet well understood. 
Supposing that the AdS/CFT correspondence is relevant for heavy ion col- 
lisions, drag calculation momentum validity is limited from above. Since 7 C 
depends on Mq but the hard pQCD scale does not there is likely a large re- 
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Figure 5.6: The double ratio of R c aa (pt) to R b AA (px) predictions for Rhic using 
Eq. (5.1) for AdS/CFT and WHDG [152] for pQCD with a range of input pa- 
rameters. While the hardening of the production spectrum reduces the dramatic 
bunching at Rhic as compared to Lhc, the lower temperature at Rhic means the 
AdS/CFT drag formalism is applicable to higher momenta. Note that R cb is plotted 
to only 50 GeV for Rhic. 

gion of pt for which both approximations are applicable for bottom quarks. 
The disadvantage of studying the more robust R cb observable is that the AdS 
momentum reach is limited by the much smaller charm quark mass; an over- 
lapping momentum region of validity for both coupling limits may not exist 
for this observable. Further careful study of the "speed limit" and higher order 
corrections to the AdS / CFT result must be done in order to fruitfully compare 
R AA and R cb results to experiment. Or, turning this around, experiment might 
inform theory as to the correct scales for which these theoretical predictions 
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Chapter 6 

Heavy Quark Photon 
Production Radiation 

6.1 Introduction 

The confidence in the application of perturbative QCD methods to jet energy 
loss in heavy ion collisions gained from the early quantitative understanding of 
7r and 7] suppression with null direct-7 control [238, 277] has recently come into 
serious doubt [438]. Evidence from measurements such as high-p^ correlations 
[417, 439] and nonphotonic electrons [423, 424] demonstrates clear disagree- 
ment with perturbative models [152, 357, 374, 435, 440]. Several papers postu- 
late alternative nonperturbative energy loss mechanisms [250, 342, 343, 436], 
and a new measurement, the double ratio of charm to bottom nuclear mod- 
ification factors, has been suggested as a robust observable for testing some 
of these novel ideas [441]. Electromagnetic radiation from quark jets, as it is 
transparent to the partonic medium, holds enormous promise as a new tool for 
investigating jet energy loss mechanisms; naively one expects the spectrum of 
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photons emitted from a jet that underwent a smooth, exponential slowdown 
will differ greatly from one that suffered the emission of a few hard gluons will 
differ from one that experienced a large number of soft scatterings. Nascent 
data of 7-hadron correlations from p + p collisions additionally motivates the 
theoretical exploration of photon bremsstrahlung in heavy ion collisions [442] . 

In this paper we calculate as a warmup problem, and ultimately as an inter- 
esting problem in its own right, the th order in opacity energy loss of a heavy 
quark jet, the radiation associated with the production of a hard parton. We 
will generalize the problem of a zero mass quark emitting a zero mass photon 
to a massive quark emitting a massive photon. The lack of theoretical con- 
sistency in the understanding of light and heavy flavor jet suppression makes 
massive quark calculations of especial interest [152, 374, 440]; moreover heavy 
quark predictions from pQCD will be necessary for comparison to AdS/CFT 
heavy quark drag results [441]. For the case of th order emission QCD and 
QED are identical but for the replacement of gcem with a s and a color Casimir. 
Using a massive photon will allow a comparison to already published results 
on the QCD Ter-Mikayelian effect [305], whose main results were: (1) the Ter- 
Mikayelian effect leads to a large reduction in th order energy loss (~ 30% for 
charm quarks), (2) the full 1-loop HTL gluon propagator can be well approx- 
imated by using a fixed gluon mass m g = = n/y/2, and (3) the small-x, 
soft gluon, number distribution for th order in opacity is 
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where, as usual, a bold variable represents a transverse two-vector. 

6.2 Calculation 

As a first step to compare to previously published results [305] we found the 
number distribution of emitted photons when simply plugging in massive 4- 
vectors into a standard classical E&M calculation: 
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where to get from the first to the second line we used the completeness relation 

Z/A=l 2 e A e A* ~~ — anC ^ wnere we took 
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(6.5) 
(6.6) 
(6.7) 



with brackets indicating lightcone coordinates. 

There are two elements seen in Eq. (6.4) and not in Eq. (6.1): several 
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factors of (1 — x) 2 , and an m 2 in the numerator. The first makes no difference 
in the limit of small x; however energy loss calculations integrate over all x, 
and it turns out that neglecting these factors is a large effect. The second 
simply cannot be reconciled with Eq. (6.1). Interestingly this extra mass term 
in the numerator fills in the "dead cone," the region of small angles with 
respect to the jet axis for which dN g /dx — > in Eq. (6.1) as k — > when 
Mq 7^ 0; this motivates additional study because naively the dead cone leads 
to a reduction in heavy quark energy loss, inconsistent with the observation of 
similar suppression patterns for pions, decay fragments from gluons and light 
quarks, and nonphotonic electrons, decay fragments from heavy charm and 
bottom quarks. 

One may rightly object that the results of Eq. (6.4) were derived using the 
usual massless photon E&M formulae. Surprisingly the only modification of 
Eq. (6.2) when using the Proca Lagrangian is to change the polarization sum 
to include the longitudinal mode. It turns out that the extra terms generated 
by the application of the identity J2\=i 2 3 

e A e A* = ~~ 9pv + k^k u /m 2 exactly cancel, and Eq. (6.4) is also valid for massive 
photon radiation. 

In order to understand the discrepancy from the field theory perspective, 
consider the diagrams contributing to the th order shown in Fig. 6.1. Evalu- 
ation of these leads to 

. , „ , x \2p ■ e* + ft . , , , 2p' ■ e* + ff\ , „ . . 
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Figure 6.1: The two diagrams contributing to the order in opacity pho- 
ton/gluon radiation spectrum. Note the inclusion of the radiation from the away 
side jet, which is usually ignored in pQCD calculations. 

where we have taken Aio(p + k,p') « Mo(p,p' + k) « A^o i n the small x, soft 
radiation limit. Most pQCD calculations ignore the away-side jet; one can 
easily see that the second term in Eq. (6.8), corresponding to the inclusion 
of the second diagram in Fig. 6.1, is crucial for preserving the Ward identity 
[234]. Simultaneously dropping the f, in the numerator and the m 7 in the 
denominator (consistent with the soft photon limit) exactly reproduces the 
classical Proca result. We note that assuming M. commutes with and 
retaining m 7 ^ in the denominator of Eq. (6.8) results in a dN P QCT>/d 3 k 
with leading order identical to Eq. (6.4) but with (1 — x) 2 — > (1 — x/2) 2 as the 
pref actor of the m 7 in the numerator and (1 — x) 2 — > {1 — x) as the pref actor 
of the ?77™ in the denominator. 
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6.3 Size of Effects 



We wish to investigate quantitatively the effect of these extra terms on the th 
order energy loss. To do so we enforce physicality by restricting the x and k 
integration limits so that the emitted photon has _E 7 > m g and leaves the jet 
with Ejet > Mq. For ease of comparison with [305] we set \i = .5 GeV and 
a = .5 fixed. One can see from Fig. 6.2 the large (50-150%) effect on AE/E 
of including the overall prefactor of (1 — x) 2 . Filling in the "dead cone" makes 
only a small difference to the energy lost (5-20%); this is a surprise as the 
"dead cone" is the usual naive justification for heavy quarks having smaller 
radiative energy loss than light quarks. 
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Figure 6.2: (Color online) th order radiative energy loss for (a) charm and (b) 
bottom quarks. All results are to leader order (LO) in \/E + . One sees that the 
largest effect (50-150%) comes from including the (1— x) 2 prefactor and that filling in 
the "dead cone" with the massive photon is a rather small one (5-20%). Comparison 
with m 7 = yields the magnitude of the LO Ter-Mikayelian effect (10-40%). 
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Fig. 6.3 demonstrates the effect of including all the terms generated by Eq. 
(6.4), not just the LO in 1/ ' E + . Of course at higher E and pr the additional 
terms make little difference, but they regulate the otherwise divergent results 
in Apx/pr as px — > 0. The Ter-Mikayelian effect, given by the difference 
between the m 7 ^ and m 7 = plots in Fig. 6.3, varies from 10-40% for 
charm and bottom energy loss. 




(a) (b) 



Figure 6.3: The effect of including all terms from Eq. (6.4) instead of just the 
leading order (LO) terms in 1/E + for (a) AE/E and (b) Apx/pr (the legend 
in (a) applies to both plots). For AE/E, the size of the relative difference in 
magnitude — the Ter-Mikayelian effect — is changed little while the overall normal- 
ization is significantly altered at low energies. For Apx/pT both the relative and 
overall normalizations change quite a bit, with the inclusion of all terms regulating 
the pt —> divergences in the vacuum production radiation spectrum. 
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6.4 Conclusions 

Unfortunately after many years of effort there is still no single satisfactory 
energy loss model for heavy ion collisions at Rhic. This leads to the need 
for basic experimental tests of the gross features of the underlying energy loss 
mechanism, whether it be more like pQCD, AdS/CFT, or some other approxi- 
mation. Medium induced photon bremsstrahlung has the potential to provide 
unprecedented insight into the modes of jet energy loss, and in this paper 
we took the intermediate step of analyzing the th order in opacity produc- 
tion radiation energy loss. While enforcing gauge invariance by not neglecting 
the away side jet fills in the "dead cone," this ultimately has only a small 
effect on the radiation spectrum. On the other hand neglecting the overall 
factor of (1 — x) 2 in the emitted photon distribution makes a surprisingly large 
difference. This prefactor, also neglected in medium-induced gluon radiation 
derivations [274, 278], may significantly alter Raa(pt) calculations, especially 
at smaller momenta. 
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Chapter 7 

Conclusions and Outlook 

Experiment drives physics. Our description of natural laws evolved not from 
aesthetics but falsifiable predictions compared to data. Scientific advance led 
us to the Standard Model of particle physics, with SU (3) the gauge group of 
the strong force coupling to fundamental quarks of fractional charge. Every 
tractable theory of nuclear physics points to a new paradigm at scales on the 
order of Aq CD , at which point quarks and gluons are presumably liberated from 
their hadronic prisons. It is clear that the prerequisite densities are reached 
at Rhic, but many questions remain, theoretically and experimentally. 

If one were to sum up the discoveries at Rhic in two words they would be 
flow and jets. The low-pr azimuthal anisotropy appears to saturate the ideal 
hydrodynamics bound. High-j9 T jets are quenched, and, shockingly, heavy 
quarks are suppressed at a level similar to light quarks and gluons. Addi- 
tionally the intermediate-pT ^2 is surprisingly high for both light hadrons and 
nonphotonic electrons. 

Hydrodynamics results hold the tantalizing possibility for direct evidence 
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of the degrees of freedom and the equation of state of the medium. These are 
of course the golden measurements for demonstrating deconfinement, thermal- 
ization, and a QCD phase transition. Heavy ion collisions may also form the 
most perfect fluid and provide the most stunning example of the application 
of string theoretic techniques to phenomenology. Hydrodynamics is a tough 
business for even qualitative discoveries, though, as so little of the process 
is under theoretical or experimental control. The initial state, so crucial in 
determining the final state, is up for grabs. Viscous corrections during the 
fluid evolution phase are not well understood. And the transition to detected 
particles through low-momentum hadronization is nonperturbative and not 
independently measurable. 

High-pr physics has the advantage of an initial state well under control 
theoretically and experimentally, although current techniques at best have 
significant error for heavy quarks. Similarly the fragmentation into hadrons is 
an independently measured vacuum process that necessarily translates to light 
meson production, although this may not hold for heavy quarks. The chal- 
lenge comes in deriving a confident theoretical description of the in-medium 
energy loss mechanisms. With this in hand one has the possibility of inverting 
the energy loss process using the theory to quantify properties of the medium, 
yielding jet tomography. The original hope that a single uncontroversial the- 
ory of final state energy loss would emerge from the theoretical community 
and immediately agree quantitatively with data has been dashed. There are 
currently four different major pQCD models of radiative energy loss and a 
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number of elastic energy loss approaches. Even more disturbing, none of these 
gives even a consistent qualitative agreement with data. 

In this thesis we confronted the two major areas of disagreement: the si- 
multaneous description of (1) intermediate-p^ Raa and v 2 and (2) nonphotonic 
electron and light meson Raa- For the first we proposed a focusing mechanism 
that provides the necessary increase in v 2 and also decreases as a function of 
momenta. Future v 2 results, should they persist at such high percentages, 
would cast serious doubt on the applicability of perturbative methods to light 
mesons. On the other hand this might point to a medium geometry signifi- 
cantly different from our expectations, with radical consequences for hydrody- 
namics results. 

For the second we included both radiative and collisional energy loss pro- 
cesses with realistic geometrical modeling. This decreased the discrepancy 
between the presumed heavy quark and light quark and gluon energy loss 
as inferred from measurements at Rhic. We also showed that this model is 
sensitive to changes in its input parameter, the density of scattering centers. 
However the significant a s dependence strongly suggests that the irreducible 
theoretical error will be large, which would reduce jets to mere qualitative 
probes of the medium. It is clear that future work in estimating the mag- 
nitude and possible minimization of major sources of theoretical error is of 
critical importance. 

Working under the assumption of jets as qualitative probe we looked for 
a qualitative test of energy loss mechanisms and found one in the momen- 
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turn dependence of heavy quarks. Specifically at Lhc dR c AA jdpx > implies 
perturbative loss whereas dR c AA /dpT < implies AdS/CFT drag loss. While 
the substantial suppression of heavy quarks at Lhc predicted by some models 
could make this difficult to see experimentally the ratio of the two nuclear mod- 
ification factors, R cb (pT) = R c aa (pt) / Raa(.Pt) is extremely robust to changes 
in the theoretical input parameters. Perturbative predictions take R cb — > 1 for 
large momenta while the double ratio is momentum independent and signif- 
icantly below one, ~ .2, from AdS/CFT drag. The only issue is the regions 
of self-consistent applicability of the two approaches. pQCD is valid at large 
momenta, but the momentum below which its approximations begin to break 
down is not quantitatively known. On the other hand AdS/CFT drag is valid 
at small momenta, but the momentum above which its approximations be- 
gin to break down is not quantitatively known. It is likely that for bottom 
quarks there is a momentum range for which both are applicable; for charm it 
is less clear. In this sense the pr behavior of the R cb ratio may indicate either 
the falsification of one or both sets or ideas or may simply give an experi- 
mental indication of these aforementioned momentum cutoffs. Predictions for 
RHIC, where the lower multiplicity will increase the momentum cutoff from 
the AdS/CFT drag derivation, were given. 

Continuing to pursue qualitatively different predictions from energy loss 
formalisms we began to investigate heavy quark photon bremsstrahlung. Even 
the intermediate results from production radiation were quite interesting: in- 
cluding the interference terms from the away-side quark jet filled in the radia- 
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tive 'dead cone.' Not neglecting factors small in x led to a significant reduction 
in the production radiation, although the Ter-Mikayelian effect remained simi- 
lar in magnitude to previous computations. The importance of including these 
factors in gluon bremsstrahlung radiation calculations is not yet known but is 
of course of crucial interest. 

The future of jet tomography in heavy ion collisions is far from rosy. We 
demonstrated in this thesis that very interesting physics can be explored in 
a qualitative way using high-p^ probes. At the same time, making the not 
insignificant assumption of a self-consistent model of energy loss in qualitative 
agreement with data, the current, admittedly rough, estimates of theoretical 
precision are low. Future experimental results will be invaluable in guiding the 
focus of theoretical research on the pertinent approximation regimes in nuclear 
physics; future work in quantifying, then reducing theoretical error may yet 
prove jet tomography possible. 
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Appendix A 



Nuclear Geometry 



For a given nuclear density distribution p(f, normalized such that j d 3 xp = 
A, one defines the nuclear thickness function T A {x = {x, y}) = J dzp(r); see 
Fig. A.l for an illustration of the heavy ion collision geometry. Then the 
nuclear overlap as a function of impact parameter b is T^(x; b) = <jnnTa(x + 
b/2, d)Ta(x — 6/2, y); aNN is the nucleon-nucleon inelastic cross section, taken 
here as 42 mb = 4.2 fm 2 . This gives the density of binary collisions in the 
reaction plane given b. The participant density is given by the probability 
that there is at least one interaction between a nucleon from nucleus A and 
the entire stack of nucleons it passes through in nucleus B, 



PpartO?; b) 



T A (X - b/2, y) (1 - e -°NNT A {x+b/2,y)^ 
+T A {X + 6/2,3/) (1 - e -°NNT A { X -b/2,y)^ _ 



(A.l) 



The total number of participants is N pa _ rt = J d 2 xp pabTt . 

Centrality is an experimental measurement associated with the impact pa- 
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y 




Figure A.l: An illustration of an A + A heavy ion collision at impact parameter 
b. The soft medium particles tend to be distributed similarly to the participant 
density, p par t while the hard collisions scale like the binary density, p co u. Figure 
adapted from [443]. 

rameter of a collision. Theoretically it is found from the inelastic cross section 
as a function of impact parameter, 



The impact parameter 6 cent for a given centrality of C% is defined by 





(A.3) 




(A.4) 



Then the impact parameter for a given centrality class C\ — C%% is found from 
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the average centrality: 



'cent d(J 

o db 



tot- 



(A.5) 
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Figure A. 2: (a) Plot of da/db as a function of b for the Woods-Saxon Au 
nucleus. The shaded region corresponds to 20 — 30% centrality with endpoints at 
b = 6.6 and 8.1 fm; the single representative impact parameter for this centrality 
class, found by properly weighting from Eq. (A.5), is b = 7.4 fm. (b) Plot of the 
nuclear density as a function of radius separately for protons and neutrons in 208 Pb 
[444]. 

The nuclear density profile can take many forms. The simplest numerically 
is the hard cylinder geometry, which is constant in density over the reaction 
plane: 

Phc(x, y) = PoO(V% 2 + V 2 - Rhc). (A.6) 
This of course leads to participant and binary densities that are also constant 
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across the reaction plane. The next simplest is hard spheres: 



PHs(r) = PoO(\f\ - Rhs)- 



(A.7) 



The Woods-Saxon distribution is a commonly used diffuse nuclear density, 
parameterized by 



For 197 Au [445] found R = 6.36 fm and z = .535 fm by comparing the charge 
distribution as measured by high-energy electron scattering to predictions from 
the Dirac equation. Values for 63 Cu can be found in [446]. These experimental 
determinations of the nuclear densities can become nearly arbitrarily complex; 
see Fig. A. 2 for the results of fits of the proton and neutron distributions with 
9 free parameters each [444]. 




(A.8) 
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Appendix B 
Power Law Raa 



B.l Fragmentation Functions 

The fragmentation function D^(z), where z = 2E h /Q, gives the probability of 
producing hadron h of energy E h from parton q in a process of total energy 
Q. If we take the problem as the production of a quark and an antiquark and 
ignore all masses then, in the center of mass frame, each parton carries half 
the total energy, 

p q = E q = Y,E h = Q/2. (B.l) 



Therefore 

z = 2E h = E h = Ph 

Q Pq Pq 

Conservation of energy then demands 



fB.2) 



j:[ 1 Z D h q (z)dz = J2 2 ^ = l. (B.3) 

u JO u V 
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From the definition of a fragmentation function we have that 



dN h {p h ) = dN q {p q )D{z)dz. (B.4) 
dN h , s f dN q 



i "(Pa) = / -^(Pi) D< y z ) dz 
dph J dp h 



where we get to the last line by employing Eq. (B.2) twice. 



B.2 Partonic R\ A 

For in-medium momentum loss in which the final momentum is related to the 
initial momentum by pf = (l — e)pi we have the probability of fractional energy 
loss for parton q with initial momentum pi is P q (e\pi). Then, exactly as for 
fragmentation functions (Eq. (B.4)), 

dN q ( Pf ) = dN 9 (pi)P(e\pi)de, (B.6) 

where we have suppressed the q index on P, and hence (Eq. (B.5)) 




(B.7) 



If we assume that the production spectrum can be approximated by a 
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power law, 



~^p~ {Pi) ~ p»w' (R8) 



where A is some normalization constant, then we may find a simple equation 
for the partonic nuclear modification factor, 

R aa(Pt) = at dPT dN " - (B.9) 

^coll^f (PT) 

N co nJde {pT/i J n{PT/1 _ e) ^-P(e\^- e ) 



N A 

JV Coll n(p T ) 



Mpt) 



tfc p( e |^L)(l _ (b.io) 

FT 

In the case of a slowly varying power law u(pt) is approximately independent 
of pt'i if we further assume a slow variation of P with respect to pt we have 



i^> T ) ~ J deP(e)(l - e) n ^-\ (B.ll) 



We note that if the production spectra are dN^ md / 'd 2 p T or dN^ rod /dp^ then 
the exponent in Eq. (B.ll) goes from u(pt) — 1 to ri(px) ~~ 2. 
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B.3 Hadronic R h AA from Partonic R q AA 

Since 



^(Pt) = ^(pr)^(Rr) (B.13) 

RaaM ~ w^-im • (R14) 
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Appendix C 

Statistical Mechanics Results 

We wish to relate the plasma density to its temperature. In general 

tat t dVd 3 p . _ . 

dN = gn £ dr = gn £ j^-^, (C.l) 

where n £ = [exp ((e — /i)/T) ± 1] _1 is the occupation number (— taken for 
fermions, + taken for bosons), and g is the degeneracy. In an ultrarelativistic 
gas the single-particle energy is e = cp. Therefore d 3 p = Anp 2 dp = ATce 2 de/c 3 
and 

- & T * 2m (a3) 

for a Fermi gas at zero chemical potential, /j, — 0. Then the density is 
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For a Bose gas at /x = we have, as the e integral yields T 3 (3/2)£(3), 



#bC(3) T 3 



V PB An 2 (hcf 



For a plasma of both fermions and bosons N tot = N F + N B and 



Ni 



tot 



V 



P 
T 



C(3) 



7T 2 (7lc) 3 

7T 2 (ft C ) 3 



9f + ^9b 



rp3 

1/3 



C(3) [9f+19b] 

n 2 (hc) 3 



1.202 [2(iV c 2 - 1) + 3N c N f ] 



P 



1/3 



(C.5) 



(C.6) 
(C.7) 
(C.8) 



where g F — 2 ■ (iV 2 — 1) for gluons in SU(N C ) with 2 helicities, g B = 2 ■ 2 ■ 
N c ■ Nf for spin-1/2 quarks and anti-quarks of Nf flavors and N c colors, and 
C(3) ~ 1.202. For a Bj or ken-expanding QGP with a boost-invariant rapidity 
distribution dN/dy of transverse area A T the temperature is then (N c = 3) 



T(r) = 



7T 2 (ft C ) 3 diV/djA 

1.202(16 + 9^] A T r J 



1/3 



(C.9) 
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Appendix D 

Formation Time and Brick 
Estimates 

D.l 7 C for Bethe-Heitler to LPM Regime 
D.l.l Formation Physics 

By the Heisenberg Uncertainty Principle, the formation time, rj., for a gluon 
(or photon) of energy k° separating from its parent quark (or gluon) — without 
loss of generality taken as a quark from now on — in the rest frame of the quark, 
is 

where the primed coordinates are in the rest frame of the quark and unprimed 
coordinates will be in the lab frame. To find ry in the lab frame, one must first 
boost the gluon energy in the lab frame into the rest frame of the quark, then 
boost back to the lab frame. Take a quark with 4-momentum p = (p°,p z ,p) 
and gluon of momentum k = (k°, k z , k). The velocity of the quark is (3 = p z /p°. 
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The energy of the gluon in the quark's rest frame is then k ' = ^y(k° — f3k z ). 
Boosting the formation time, Eq. (D.l), back to the lab frame will dilate the 
time by 7; therefore in the lab 



where for the last line we have taken k° ~ k z so that Eq. (D.3) reproduces 
[290]. 

For a quark of mass M and an emitted gluon of mass m g with momenta 

p = U/M 2 + (1 - x) 2 p 2 z + k 2 , (1 - x)p z , -k) (D.4) 
k = (y\J m 1 + x 2 P 2 z + k 2 , xp z , \tj , (D.5) 

the formation time to lowest order (assuming large p z and small k) is 

2x(l — x) 2 p z 



Tf ~ k 2 [l - 2x(l - x)} + (1 - x) 2 m 2 + x 2 M 2 



(D.6) 



this agrees with [447] in the limit of small-x considered there. Note that using 
Eq. (D.3) one finds (also for large p z and small k) 

2(1 - xfp z 

Tf * x(k 2 + M 2 Y (D - 7) 
which agrees with Eq. (D.6) in the large-x Bethe-Heitler limit. 
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For the same momenta, Eq. (D.4), the leading order lightcone coordinate 
representation, 



M 2 4- k 2 
/If 2 4- k 2 

* [xE + ,^^,-k] (D.9) 



yields exactly the same formation time, Eq. (D.6), but with p z replaced by 
E+. 

D.1.2 Bethe-Heitler to LPM 

The "speed limit" at which incoherent Bethe-Heitler turns over to coherent 
LPM, 7c, is found when the formation length for the produced radiation is on 
the order of the radiation's mean free path: 

Tf ~\ g , (D.10) 

where 77 is given by Eq. (D.2) and 



X g = (D.ll) 

op (V j T s 

~ \- ( D - 12 ) 
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For large-x Bethe-Heitler radiation one sees from Eq. (D.6) that 



T f~ (D.13) 



Eq. (D.10) then gives 



M 

7c ~ (D.14) 



D.2 AdS/CFT Brick Estimate 

We wish to more rigorously derive Eq. (5.5). Recall that Eq. (5.5) justifies Eq. 
(5.6), the back-of-the-envelope explanation for the momentum independences 
and normalization for R cb (p T ) = R c AA (p T ) / R b AA (p T ) for AdS/CFT Drag-type 
heavy quark energy loss. 

From the simple trailing string model of [342, 343] we have by Eq. (5.1) 

d P T m i k\ 

— = - m (D.15) 

where 

P 2M ' v ; 

with T the temperature of the SYM plasma and M the heavy quark mass. 
Solving Eq. (D.15) one finds that 

p(t) = p (0)e-f dt '^ t '\ (D.17) 
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Defining the fractional momentum loss e as 

Pf = {l-e)pi (D.18) 
=> e = l- e -/ d *'M*'). (D.19) 

From Eq. [] we have that for a quark (or gluon) with slowly varying power 
law production index uq(pt) + 1, 

Raa(pt) * ((1 - e) n ), eom , (D.20) 
where (• • • ) ge om indicates averaging over the geometry of the collision. 

D.2.1 Static, Constant Medium 

In a static, constant medium /i(x + vt) = /x. For this case Raa for a single, 
fixed path of length L is 

R AA {p T ) ~ e" n « (pT) ^, (D.21) 

where we have assumed v ~ constant (a good approximation for ultrarela- 
tivistic motion with v ~ 1). From now on we will simply take f ~ 1 and 
suppress it in our derivations. One can see that Eq. (D.21), plotted in Fig. 
D.l, qualitatively reproduces the full numerical results of Fig. 5.2 for the a 
posteriori input parameter choice of L = 1 fm. 

One may also readily find the nuclear modification factor assuming a uni- 
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Figure D.l: Raa{pt) for charm and bottom quarks traversing a single pathlength 
of L = 1 fm in a static medium of temperature T = 300 MeV (from Eq. (C.9) with 
the PHOBOS extrapolation to LHC gluon density dN g /dy = 1750 [87, 403], r = .2 
fm, t = 1 fm, and Nf = 3), with D = 3/2irT from Eq. (D.21). The magnitude and 
momentum dependence are similar to the results of the full numerical simulation 
with Bjorken expansion and realistic medium geometry shown in Fig. 5.2. 

form, static brick of plasma. Then Eq. (D.20) and Eq. (D.19) give 




1 



(D.22) 



The results of Eq. (D.22) are plotted in Fig. D.2 with L = 5 fm, the temper- 
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ature T ~ 210 MeV taken at L/2 — 2.5 fm, and D/2irT = 3; they are again 
in qualitative agreement with the full numerics shown in Fig. 5.2, and the a 
posteriori input parameters are somewhat less arbitrary. 




Pt (GeV) 

Figure D.2: Raa(pt) for charm and bottom quarks created uniformly along a 
static brick of SYM plasma of length L = 5 fm and temperature T ~ 210 MeV (from 
Eq. (C.9) with PHOBOS gluon density dN g /dy = 1750 [87, 403], r = .2 fm, and 
t = L/2 = 2.5 fm) according to Eq. (D.22) with D = 3/2-itT. The magnitude and 
momentum dependence are similar to the results of the full numerical simulation 
with Bjorken expansion and realistic medium geometry shown in Fig. 5.2. 
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D.2.2 Static Nonconstant Medium 

We note in passing the analytically interesting case of a time-independent 
medium of nonconstant density, /i = fi(x). Then 



— = —u(x)dt = —fi-dx. (D.23) 
p ' v 



Since p = ^mv we have that dp = mvd'j + ^mdv = (1 + 7 2 t> 2 )dv . Then Eq. 
(D.23) becomes 

v(l + 7 2 f 2 )'ymdv , _ 

_ = -fjuix, (D.24) 



Moreover 



so we have that 



jmv 



v 2 1 1 + v 

idv = -In - v, (D.25) 



1 - v 2 2 1 - v 



1, l + v 
m 



v f r x f 

fidx, (D.26) 



2 1 -v 

where the integral on the RHS is really a line integral through the medium 
from the initial starting point Xi to the final point of interest, Xf. Therefore 

l±H = l + * e -2f^ (D . 27) 
l-Vf l-Vi 

For the case of ultrarelativistic motion, 1 + Vf ~ 2, and Eq. (D.27) simplifies 
to 

Vf = 1 - (1 - Vi )e- 2 ^ dx . (D.28) 
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The full solution to Eq. (D.27) is 



1+ 



Vf = . (D.29) 

1 1+Vj g-2 / fidx i I 

1-Vi 



D.2.3 Bjorken Expanding Medium 

One might reasonably wonder the consequences of including Bjorken expansion 
on our estimate of Raa from D.2.1. Defining 

7T a/A ,„ 

v = n— —r. D.30 

2M ' y ' 

Eq. (D.20) becomes for a uniform brick of SYM plasma 



»AA 







L d A e -»s^>T\o (D31) 

^-"Wi^) + 1°, (D.32) 
L L 



where Eq. (C.9) gives 



dN g /dy 7T 2 
A T C(3)(16 + 9AM' 



«= f " ' nArV (D.33) 



The integral in the exponential of Eq. (D.31) is 



K 



>:t 



2/3 



3« 2 / 3 
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Changing variables to y = L — £ =^> dy = —di and defining to = 3k 2 ^ 3 u, the 
remaining integral in Eq. (D.31) becomes 

- *f J = *f f ^ [2 + W + ^ 



e 



wr o 1/3 f 3 e -"(^-^) 1/3 



x 



L 1 tu 3 

=L—to 

[2 + 2u(L-l) 1 / 3 + u 2 (L-£) 2 / 3 ] ' 



1=0 



1/3 /■ 1/3 

e WT o 3e _£4,T o 



w 3 



2 + 2wr 1/3 + w 2 t 2/3 



3e- Ll/3 
w 3 



[2 + 2u;L 1 / 3 + u; 2 L 2/3 ] \. (D.35) 



a; is in fact so large that for any reasonable r ~ .2 fm, ur ^> 1. Therefore 
we arrive at, to leading order, 

Raa - ^ + T = (D.36) 

Since a; oc uq/Mq, with the proportionality constant the same for charm and 
bottom quarks, we still find that the ratio of charm to bottom Raa goes as 

RUpt) - (D.37) 
n c M b 

Note the essential importance of the plasma formation time, To, in Eq. (D.36) 
in order to reproduce the linear in M c jM b behavior. The results of the full 
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calculation (not just the leading order terms) for charm and bottom quarks 
are shown in Fig. D.3. The momentum dependence is very similar to Fig. 
5.2. The normalization was made similar by a posteriori setting the Bjorken 
expanding brick to the abbreviated length L = 3 fm. 

i.o 
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0.4 
0.2 
0.0 

50 100 150 

Pt (GeV) 

Figure D.3: Raa{pt) for charm and bottom quarks created uniformly along a 
Bjorken expanding brick of SYM plasma of length L = 3 fm as given by the full 
analytic solution to Eq. (D.31) with D = 3/2ttT. The magnitude and momentum 
dependence are similar to the results of the full numerical simulation with Bjorken 
expansion and realistic medium geometry shown in Fig. 5.2. 
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Appendix E 



Standard pQCD Results 



E.l Light Cone Coordinates 
E.l.l Definitions; Identities 

A four- vector is not bold-faced (e.g. p, k), a three- vector is bold-faced with 
a vector symbol (e.g. p, k), and a transverse two- vector is bold-faced with- 
out a vector symbol (e.g. p, k). Minkowski four- vectors are written with 
parentheses, (); light-cone four- vectors with brackets, []. 



P = (p°,P z ,P) = [P 



P ,Pj- 



(E.l) 



We will use non-symmetrized lightcone coordinates: 



p + = p + p z 



(E.2) 



P 



= p°-p 



z 



(E.3) 



P = 



P 



(E.4) 
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The inverse transformation is then 

P° = \{P + +P-) (E.5) 

P z = \ip + ~P~) (E.6) 

p = p. (E.7) 

The Minkowski dot product in lightcone coordinates is: 

p ■ k = p°k° - p z k z - p • k = \{p + k~ + p~k + ) - p • k. (E.8) 

The length of a vector using lightcone coordinates is then: 

p ■ p = p + p~ — p ■ p (E.9) 



E.1.2 Eikonality in Light Cone Coordinates 

Note that the eikonal approximation requires collinear emission, which in turn 
requires k + ^> k~ (for emitted boson with 4-momentum k). If this were 
otherwise, the k z emission would be in the opposite direction of p z : p z = 
\{p + —p~) — \p + = |(1 — x)E + ; on the other hand k z = \{k + — k~) must be 
approximately \k + or else the sign is wrong. 
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A useful relation comes from eikonality: 



k+ > £r => xE + > 



k 2 



=>- xE + > Ik 



(E.10) 



E.2 Monojet vs. Dijet th Order Radiation 



In order to motivate the work of Appendix F we will derive the order 
bremsstrahlung radiation loss for massive quarks and gluons using two meth- 
ods and find that the results are inconsistent. We will find that the usual 
pQCD approach of calculating the radiation spectrum from a monojet only 
is consistent with a dijet calculation in the massless limit only. For massive 
gluons the longitudinal polarization must be considered, and its contribution 
to the spectrum from the away-side jet will be 0(1) in an expansion in powers 



We will find the following relations useful in evaluating Feynman diagrams 
related to bremsstrahlung, which we will refer to as Peskin's Trick: 



of 1/E+. 



E.2.1 Peskin's Trick 



(f> + m)Ye*AP) = 2P"e£u(P). 



(E.ll) 
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The first relation is proved as follows: 



u{p')l ■ e*(/ + m) = u(j/)'fe*('fti + m ) 



u(p')(2e* -p' - (/ + m)7-e*) 



(E.12) 



The second line comes from using {7^,7^} = 2g ,il '; the third line results from 
the Dirac equation, u{p)(jf> + m) = 0. The proof of the conjugate relation 
follows the exact same line of reasoning. 



E.2.2 Diagrammatic Calculation, Massless Photon and 
Quark Monojet 



Jfc,c 




Qe 




P 



Figure E.l: Zeroth order radiation diagram. The blob on the left represents the 
factorized production of the hard quark. 
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There is no scattering off a medium center. Hence 



P ~ (E, E, 0) = [E + = 2E, 0, 0] (E.13) 

P = K 1 -^ {1 -l)E^ - k] (K14) 

k = ^^ = 2^ = 2^^ (E.15) 

6 = [0,^,e]. (E.16) 



Note that P is off-shell so it cannot be exactly P = (E, E, 0). p~ is found by 
enforcing p ■ p = 0; similarly for k • k — and • e = 0. Some useful relations 
used when deriving the radiation spectrum from this diagram are: 



2ke 

2p-e = p + e~ - 2pe = (1 -x)E + — — + 2ke 

xP+ 



2ke 

x 



2 



(E.17) 



(E.18) 



We assume some kind of factorization so that we don't have to explicitly 
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calculate the production process. Therefore 

iM = u(p)iQeYe^-J^—M (P,P )u(P ) (E.19) 

^ P z + 16 

~ u(p)iQe 7 % J M (P,Po)u(P ) (E.20) 

[P + K) z + 16 

- ^ e o • u(P)M (P,P )u(P ) (E.21) 

2k • 6* 

= ^e-^n(P)>l (P,P ) M (P ) (E.22) 



The first line comes from writing down the result of the diagram, where A^o 
denotes the part of the amplitude corresponding to the hard production. In the 
next line we approximate the numerator of the propagator P with p (dropping 
k because it is small). The third line results from exploiting Peskin's trick, 
Eq. (E.ll), and approximating u(p) by u(P). 

We wish to calculate the distribution of emitted photons given that pro- 
duction has already taken place. Since production is represented by uAiou, 
we will drop this term in our further analysis. 

The quantity of interest is then 

El^l = ^El k '<| 2 - (E.23) 

A A 

We may choose our polarization vector however we may wish; without loss of 
generality take e\= x and e* 2 = y. Then 

^|k-e*| 2 = k 2 (sin(0) 2 + cos(0) 2 ) . (E.24) 

A 
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Therefore 



4(Qe) 
k 2 



(E.25) 



E.2.3 Soft Photon Approximation, Massless Dijet 

Compute electron (mass m — 0) beta decay starting with Eq. (F.84), 



iM=tt(p) [M (p,P))u(P) ■ 



Qe 



P 



p ■ e 



p ■ k P ■ k 



(E.26) 



We will want to exploit the replacement e^e* — > —g^ u so we will have to 
preserve the Ward identity (charge conservation), 



ki, 



p' J 



k ■ p k ■ P 



0. 



(E.27) 



To do so we will have to keep the second term explicitly with Pm — [M, M, 0, 0] 
being the momentum of the "proton" after (as well as before) beta emission. 
Then (essentially Eq. (F.86)), 



{Qef{-g, v ) 



pi- 



k ■ p k ■ Pm ) \ k ■ p k ■ P J 



pu 



M 



(Qef 



2p-P, 



M 



M 2 



m 



(k-p)(k-Pj 



Ml 



{k ■ Pm) 



(k ■ p) r < 



(E.28) 
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We will need the following dot products: 



p-P = [(l-x)E+, ^ -k]-[M,M,0] 

~ \ME + (E.29) 
1 ^ k2 - 2^\ , 1,2 k2 



p.* = - [~ + ^ 2 J+k (E.30) 

1 / k 2 \ 1 

k ■ P = - MxE + + M— — p ~ -MxE + (E.31) 



With these in hand we have that 



V|^| 2 = (Qe) 2 <| / ^) ME+) , - n ^ (E.32) 



4(g e ) 2 4(g e ; 



2 



k 2 (x£+) 2 

4(g e ) 2 



(E.33) 
(E.34) 



where the final relation comes from k <C xE + due to eikonality, Eq. (E.10). 
So 



A 



^iVf = (Qe) 2 _ Q 2 a 
47r 3 wk 2 7r 2 o;k 2 

(E.37) 
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We wish to transform coordinates 



k x = k ± cos((j>) (E.38) 
k y = k ± sm{<f)) (E.39) 



k 2 -k 2 ± ~k (E.40) 



Thus the Jacobian gives d 3 k = k±dkdk±d(j) = k±dcodk±d(j). Hence 

d 3 iVf Q 2 a 



2 



k±dudk±d(f) 7i 2 ujk 
d 2 N\ 0) 2Q 2 a 



k±duidk± nuk 2 
d 2 N^ ] Q 2 a 
dudk 2 nuk 2 

Nicely, we recover the Sudakov double logarithm [251]: 



(E.41) 
(E.42) 
(E.43) 



N (0) = 9^ log (^A h J . (E.44) 

E.2.4 Diagrammatic Calculation, Massive Photon and 
Massive Quark Monojet 

The idea is to reproduce Djordjevic's calculation of zeroth order radiation in 
[305]. The diagram is the same as Fig. E.l except that the quark has a mass 
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Mq and the photon has mass m T The relevant four- vectors are now 



M 2 



p = 


[(1-x) 






k = 


[xE + , 




. ek 


e = 





k 2 + M 2 



Q 



xE+ 
,el. 



l-x)E+' 



-kl 



(E.45) 
(E.46) 

(E.47) 
(E.48) 



The full matrix element is then 



%M = u(p)iQe^el ^ + ^. M (P,PMPo) (E.49) 



i 2 Qe 



2p-e* 



2p ■ k + m 2 + ie 



u(P)M (P,P )u(P ), (E.50) 



where Peskin's trick, Eq. (E.ll), was used. The latter amplitude corresponds 
to factorized production and will be dropped. We need the following relations 



2p ■ k + m. 



[l-x)(k 2 + m 2 ) x(k 2 + M ( 



x 



1-x) 



+2k 2 + ml 



k 2 + (1 - x) 2 m 2 + x 2 M 2 



x(l — x) 



k 2 + (l-x)m 2 +a; 2 M 2 



2p-e 



x(l — x) 
2(1 -x)k- e 



+ 2k-e 



x 



2k e 



x 



(E.51) 
(E.52) 

(E.53) 
(E.54) 
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Therefore 



iMp - Q e k 2 + {1 _ x)m 2 + ^ M 2 /x{1 _ x) (E.55) 

(1 — sc)k • e 

= 2ge k2 + (1 _ x)m 2 +x 2 M 2 (E.56) 

k ■ e 

~ 2ge k2 + ( 1 _ x ) m 2 +x 2 M 2- (E.57) 



Hence the distribution of emitted photons is 
dN\ 0) Q 2 a k 



d 3 k n 2 u [k 2 + (1 - x)m 2 + x 2 M 2 } 2 ' 



(E.58) 



This does not agree with [305], but does agree with the — presumably cor- 
rected version — in [447]. Note that we have not included radiation into the 
longitudinal mode. 



E.2.5 Classical Calculation, Massive Quark Jet and 
Massive Photon 

Note that while the sum over polarizations changes to J2\ e n e t = ~9^u + 
k^ky/m 2 , Eq. (E.28) is unchanged; the additional terms exactly cancel each 
other. We need to calculate the relevant dot products using the four momenta 
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from the previous section. One finds that 

M 



P-Pi 



M 



2(1 -x)E* 



((l-x)E + y + k 2 + M 2 (E.59) 



p-k = 2x{ i_ x) [k 2 + (l-z) 2 m 2 +z 2 M 2 ] (E-60) 
P M -k = ^[(xE + ) 2 + k 2 + m 2 }. (E.61) 

The result for Eq. (E.28) expanded to first order in 1/ E + is 



^ , 4(Qe) 2 (l-x) 2 (k 2 + (l-x) 2 m 2 ) , ~v 
^ 1 1 [k 2 + (1 - a;) 2 m2 + x 2 M^] 2 V / 



There are three differences between Eq. (E.58) and Eq. (E.62): the addi- 
tional overall factor of (1 — x) 2 in Eq. (E.62), (1 — x) 2 in Eq. (E.62) instead 
of 1 — x in Eq. (E.58) as the coefficient of m 2 in the denominator, and the 
additional (1 — x) 2 m 2 in the numerator of Eq. (E.62). In the limit of small x 
for which the calculations are valid, the first two disappear; however the third 
is a true inconsistency. In the following Appendix we will find that the second 
is due to including the mass of the gluon from the k 2 term in the denominator 
of the quark propagator and that the first and third are due to the interference 
terms from the away-side jet. 
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Appendix F 

Heavy Quark Production 
Radiation Derivation 

F.l Classical Computation 

F.l.l Variation of the Determinant of the Metric 



This derivation is due to Simon Judes. 

Let G = Q nit be the metric matrix, and let det G = g. The inverse matrix 
is G- 1 = g^; necessarily detG" 1 = g' 1 (det(M" 1 ) = (detM)" 1 ). 

From the definition of the metric, 



Proof of Identity 



First we must derive a very important relation: 



d^g = gg a ^ d^g a p. 



(F.l) 



ds 2 = g^dx^dx' 



(F.2) 



Appendix F: Heavy Quark Production Radiation Derivation 



223 



it follows that G can be chosen to be symmetric without loss of generality. It 
can therefore be diagonalized by a unitary transformation. Define the diago- 
nalized metric as 

g = UGU~\ (F.3) 

where U is the unitary matrix whose rows are the (normalized) eigenvectors 
of G. Note that 

det Q = detiUGU- 1 ) = det(U^UG) = det G. (FA) 

Consider 

d»g = d^ detG . (F.5) 

Since detG = det Q, from above, and, trivially from the definition, In det Q = 
Tr \nQ, we have that (where Qi is the ith eigenvalue) 

i 

i 

= gTrQ-'d^g. (F.6) 
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Now we must prove that g^d^g^ is equal to this last quantity: 

= Tr UG^U^Ud^GU- 1 

= Tr^T 1 [d^G - (d^GU- 1 - UGid^U- 1 )] 

= TiQ^d^g - UG^U-^lTjGU- 1 - UG^U^UG^d^- 1 ) 

= Trg-^g-u-^d^-u^u- 1 ) 
= Trg-^g-d^u^u) 

= Trg-'d^g. (F.7) 
Use of Identity 

We will use an important example of the identity Eq. (F.l): 

Sg = gg^Sg^ = -gg^gT- (F.8) 
The last equality comes from taking the variation of g^g^ = K = #: 

8{9^gn = g^Sg^ + g^Sg^ 
^{TSg^ = -g^SgT. (F.9) 
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Therefore we find an identity we will use often in the following: 



99»u = 7T9i»>\ ( F - 10 ) 



8g^ v 
similarly 



F.1.2 T^ v from Varying the Metric 

We will use variation with respect to the metric to find the (automatically) 
symmetric, gauge-invariant form of T^ v [319]: 



where g is the determinant of the metric and the coefficient out front gives the 

proper normalization for the usual stress-energy tensor (assuming a H 

type metric). 

E&M 

As a warmup exercise we use Eq. (F.12) to derive the usual E&M energy 
density from its Lagrangian, 



£em = ~\f, u F^ = - l -g«P g ^F ai Fp S . (F.13) 
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Therefore 



t-iEM 
fj,u 



5 



gdg^ 
1 



1 r - 

jV .'/// .'/ ' n- 

^9 



9,uF 2 + ^g-6°6& 5 F aj Fp 



+ V^gg^^lF^Fs 



1 



F ij,\F\ + -g^FapF 1 * 13 . 



(F.14) 



We know that F 2 = 2(B 2 + E 2 ), and since E = -F° l we have that 



rpOO 

EM 



F 0i F 4 + -{B 2 -E 2 ) 

i , - 



= ^(E 2 + B 2 ). 



(F.15) 



Proca 



Let's now find the stress-energy tensor for the Proca Lagrangian, 



--Proc — £em + ^"^At 
1 u 2 



(F.16) 
(F.17) 
(F.18) 
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where /x is the mass of the photon. 

Consider the variation of the new, massive photon part of the Lagrangian: 



2 S 



[Ml/ 



-gg aP A a A p 



-g^A z + V^gA^A 



fi z [ A^A V - —g^AAx 



(F.19) 



Therefore 



Proca 



F» X F X U + -g^F^F^ + /i 2 (a*A v - -g fMU A x A x 
4 V 2 



(F.20) 



And 



rpOO 

Proca 



[E 2 + B 2 + fi 2 (A° 2 + A 2 )]. 



(F.21) 



This agrees with Jackson problem 12.16 (a) [235]. Note that we recover 
for \i — > 0. 

F.1.3 Classical Bremsstrahlung Radiation 

E&M 

The E&M Lagrangian coupled to a classical source is 



C ESz m = -^F iw -^A^ 



(F.22) 



Appendix F: Heavy Quark Production Radiation Derivation 



228 



Using the E-L equations (d a dC/ dd a Ap — dC/dAp = 0), we find that 



d a F aP = f (F.23) 



d n d a AP - d^d a A a = ^ 



d a d a A? = f, (F.24) 

where we get from the second to the third line be choosing to be in the Lorentz 
gauge (d a A a = 0). 

For the Beta decay problem, following Peskin 6.1 (pgs. 177-180) [251], our 
external current takes the form 



f{x) = e j ( x - y( T )) , (F.25) 



where, for the trajectory of our particle, we take 



(p^/m)r for r < 0: 
!/"(r) = < (F.26) 
(p' M /m)r for r > 0. 



We will work with the Fourier transform of the pertinent quantities: 

A»(x) = j ^e-^A^k) (F.27) 
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With this normalization for the transforms we have that 



i"(Jfc) = J d*ke ik - x A»(x) 
~f{k) = [ d 4 ke ik - x f(x). 



(F.29) 
(F.30) 



Solving the EOM (Eq. (F.24)) in momentum space we find that 



A"{k) = ~j"(k). 



(F.31) 



We can find j from the given trajectory (Eq. (F.26)) and the definition of 
the Fourier transform: 



f{x) 



o m 



(x- V -r\ + ef dr P -5^( 

V m J J-oo m V 



p 

x r 

m 



(F.32) 



~f(h) = / d A xj k - x) f{x) 



00 j/v r° n M 

V r i(k-p' /m+ie)T _|_ g / (jr ^{k-p I m-ie)i 

m 



e I dr — e 



le 



p< 



p' J 



k ■ p' + ie k ■ p — ie 



(F.33) 



Therefore we have the field in position space as an integral over momentum: 



A»(x) 



d A k 



-ik-x 



—ie 



p 



P J 



k ■ p' + ie k ■ p — ie 



(F.34) 



The integral on the right hand side (RHS) has the pole structure as shown 
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in Fig. F.l. 



k° = -\k\ 



kp = 



k° = \k\ k ■ p' = 



Figure F.l: Pole structure of Eq. (F.34). 

For t < we close the contour in the upper half of the plane, picking up 
the k° = k ■ p/p° pole; the result is 



AHx) 



d 3 k fa -iffi) t (2VTZ) -ie 
(2tt) 3 2tt k 2 p° 



(F.35) 



where the k° in the k 2 term is given by its value at the pole. Note that the 
l/p° comes from correctly evaluating the contour integral using 



p f ' 



k-p-ie k°p° -k-p 



16 



k° 



k-p 
T7T 



IE 



(F.36) 



We can see that the k ■ p = pole for t < corresponds to the Coulomb 
field of the particle by boosting to its rest frame; we have that k 2 = —\k\ 2 
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(because k ■ p = =>- k° = for p^ = (m, 0) M ), and Eq. (F.35) becomes 



A»(x) 



2tt 



e / rffc / sin^e ifcrcose (l,0f 
Jo Jo 



— \' — 

,7 /' 



4rrr 



(1,5)^. 



(F.37) 



For i > we will ignore the k ■ p' = pole because that merely gives the 
Coulomb field; radiation is produced by the k° = \k\ and k° = —\k\ poles. 
Since k 2 = (k° — \k\)(k° + \k\) and k ■ p = k°p° — k ■ p we have that 



A»(x) 



d 3 k (-2m) 
(2tt) 3 2tt 

1 



-ze 



e -i|fe|* e ife-5 



2|fc| V |/c|p'° — k - ft \k\p°-k-p 



-2\k\ y-|/c|p /0 — k-i? -\k\p°-k-p 



—e 



d 3 k 1 
(27r) 3 2|fc| 



e -i|Jfc|* e ife-5 



V 



\k\p'° — k ■ p! \k\p° — k ■ p 



_|_ ^\t e -ik-s 



P 



ft* 



\k\p'° — k ■ f? \k\p° — k ■ p 



(F.38) 



where, to get to the last line, we were able to freely change k — > —k because 
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of the integration over all space. Therefore we have that 

j^-,A^{k)e- ik - x , (F.39) 
[2-kY 

A»(k) = (F.40) 
\k\ \k-p' k-pj 

Note that we have implicitly taken k° — \k\ and will continue to do so for the 
rest of the calculation. 



If we define 



/H 3 k 
— f*(x)e-** (F.41) 

B\x) = Re J -0^B i (x)e- ihx (F.42) 



and recall that 



E l 
B 



gi A _ gO A i 
*F 

\j jk ($A k - d k A>) = e i jk d^A k 

(v x Ay. 



(F.43) 



(F.44) 
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we then have 



£{k) = ik°A(k) -ikA°(k) 
B(k) = ikxA(k). 



(F.45) 
(F.46) 



Let's now prove some nice relations regarding £ and B that will be useful 

— * — * 

later on: transversality and orthogonality. To prove that 8 is orthogonal to k 
recall that A^ oc (p'^/k ■ p') — (p^/k ■ p). Then 



k ■ £{k) oc 



'3C 



-i\k\ 2 (k-A-A°) 
-i\k\ 



1(1 


p 


k 






p> 


k 


i) 



■oc 



-ilk 



k°p'° k°p° 

k • p' k ■ p 

k ■ p' k ■ p 

k • p' k • p 



= 0. 



Also 



k x S = k x 



ik°A-ikA°j 



(F.47) 



= ik°k x A; 

since k° = ±|fc|, and we won't be interested in the sign, we have that 



(F.48) 



B = k x £. 



(F.49) 
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Let's use all our knowledge to find the energy radiated by this motion. 
Recall that 



EM 



a X1 EM 



+ \B(x 



(F.50) 



Consider first the electric piece: 



d 3 x\E(x)\ 2 

, d k f d k 

(Fx 



8 J ~ ~ J (2tt) 3 J (2tt) 3 
S(k')e- ik '- x + £*(k')e ik '- x 
d 3 k 



S\k)e- ihx + £*(k)e lk - : 



II 



£(k) ■ £{-k)e' 2tkH + £*(k) ■ £*(-k)e 



(2tt) 3 
+2£{k) -£*{k) 



2ik°t 



(F.51) 
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Now consider the magnetic piece: 



E-a — 



1 

2 

lf**J (2 



d 3 x\B(x)\ 2 

d 3 k f d 3 k' 



' 3 _k_ r d 3 k' 

tt) 5 J (2tt) 3 



k x £(k)e~ + k x £*(k)e 



ik-x 



k' x £{k u 
d 3 k 



k x £(k)} ■ \(-k) x S(—k) 

2ik°t 



(2tt) 3 

x £*(£)] • |"(-fc) x £*(-*) 
fc x • [jfe x £*(£) 



+ 

+2 

We will now prove that (k x £) • (k x £) — £ • £. 

e ijk k j £ k e ilm k l £ m = (5 jl 5 km - 5 jm 5 kl )k j £ k k l £ r 



(F.52) 



= {%■%){£■£)-{%■£){%■£) 

= £■£. (F.53) 



Therefore we find that the energy radiated is 



E: 



1 f d 3 k 



EM 



2 J (2tt) : 



£(k)-£*(k). 



(F.54) 
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With two transverse polarizations e\(k) we have that 



£(k) ■ S*(k) 



A=l,2 A=l,2 



k° 2 m)-A(k)\ 2 = i*i 2 £ m)-A(k)\ 2 . 



A=l,2 A=l,2 



(F.55) 



This leads us finally to 




p p 



2 



k ■ p' k ■ p 



(F.56) 



We would now like to show that one may replace the three vector quantities 
with their corresponding four vectors, then take J2 e v e t ~ * ~9^v- If we orient 
k^ = (A;, 0,0, k) 11 along the z-direction, then we may take our two transverse 
polarizations as = (0,1,0,0)^ and = (0,0,1,0)^. It is now clear that 
we may freely take the three vectors to four, as the zeroth components of the 
polarization vectors is zero. To turn the sum over polarizations into a metric 
we need that the Ward identity be satisfied, k^j^ = k^A^ = 0. Then, with 
our k oriented along z, 



k^A^ = kA + kA 3 = 



A = -A 3 . 



(F.57) 
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Then 



= -\A Q \ 2 + \A^ 2 + \A 2 \ 2 + \A Z \ 2 

= - gpv A»A w . (F.58) 

We immediately arrive at our final result: 



f d 3 k e 2 ( p'>* p^ 



p' u p v 



k ■ p' k ■ p 



f d 3 k e 2 ( 2p-p' m 2 m 2 \ 

J (2tt) 3 2 \(k-p')(k-p) (k-p'y {k-p) 2 )' 1 ' 

Proca 

When we include an external current in the Proca Lagrangian we have 



^Proca = -\F» V F» V + - j„A». (F.60) 



The EOM are then 

d a F a/5 + fx 2 A 13 = j 13 . (F.61) 
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For a conserved current, dpj 13 = 0, we have a consistency equation 

dp (d a d a A p - d a d?A a + fi 2 A p ) = dpf 

=> d a d a dpAP - d^d a A a + ^dpA? = 

=>dpA fi = 0. (F.62) 

Our EOM, Eq. (F.61), become 

d 2 A p + p, 2 A^f. (F.63) 
Fourier analysis immediately leads one to 



d 4 k _ ik . x -ie ( p'^ p^ 



A»(x) = / — e'* 1 - . (F.64) 

' (2vr) 4 k 2 -/i 2 \k -p' + ie k-p-iej 



The poles in this occur at k ■ p' — 0, k ■ p — 0, and k 2 — fi 2 = 0; the pole 
structure is exactly the same as for the massless case, Fig. F.l, but with the 
radiation poles giving a dispersion relation of 



k° ± = ±Jk 2 + fi 2 . (F.65) 



Let's examine the t < case and make sure we recover the correct Yukawa 
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potential in the rest frame. For t < we pick up only the k ■ p = pole. Then 



J (2jt)3 V 2» J -B - f P° 

= e(l,0f— ^ / k 2 dk / Sin^- jocose 



Arcr 



neT/r 

^ r (l,0) M . (F.66) 



Let's now examine the t > solution. Again we'll ignore the uninteresting 
k ■ p' = Yukawa pole contribution. We find that 



At,/ \ . -27ri f d 3 k 



2tt J (2tt) 3 

1 / j^'A* 



2/c+ V k° + p>° -k-f k° + p° -k-p 



^2 j ^^zA/"«j 



2&° \ p ,o _ £ . ^ fc o p o _ £ . ^ 
We can therefore write ^4p roca in the form of Eq. (F.39) with 



where, from now on, we will implicitly take A; = However because /c^fc^ = 
li 2 instead of as in the E&M case we cannot immediately write | E \ 2 + | B | 2 as 
£■£* as we did before. Additionally, since we do not have gauge invariance like 
in the E&M case, we do not have transversality of the fields; there are three 
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polarization vectors (whose sum in fact obeys Y^\=i e \ e *\ = ~ 5V + k^k" / fi 2 ). 
As a result we won't use any of the tricks exploited in the E&M calculation; 
rather we will plug our solution of -4.p r0 ca into the formula for -Ep roC a, Eq. 
(F.21), and brute force manipulate it. 

We are interested in finding quantities of the form J d 3 xQ 2 (x) where Q(x) = 
Re / d 3 kQ(k)exp(-ik- x)/(2tt) 3 . Then 

x (Q{k')e- ik '- x + Q*{k')e ik '- X 
-Q(k)Q(-k)e- 2ik0t + Q*{k)Q*{-k)e 2m 



4 J (2tt) 3 

+ 2Q(k)Q*(k) 

1 d 3 k 



2 (2ny 



Q(k)Q(-k) cos(2k u t) + Q(k)Q*(k), (F.69) 



where we are taking k = (k°,k) and — k = (k°,—k). We get the last line 
because, as we will see below, A* 11 = A^ =>- Q(k)Q(—k) = Q*(k)Q*(—k) for 
all the quantities of interest. 
From Eq. (F.45) we see that 



£ k -£_ k = [i(k°A k - kA° k )] [i(k°A-k + kA°_ kJ 

= k°A° k k ■ A- k + k 2 A° k A°_ k - k° 2 A k ■ A- k - k°A°_ k k ■ A k 

= k 2 A° k A°_ k -k° 2 A k - A- k (F.70) 

= £*k-£-ki (F.71) 
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and 



E k ■ El = k° 2 A + k 2 Af - 2k°A° k k ■ A k 



(F.72) 



From Eq. (F.46) we find 



Bk-B. k = [i(kxA k )][-i(kxA- k )] 
= k 2 Ak ■ A-k — k ■ Akk ■ A-k 



(F.73) 
(F.74) 



where we have used the vector identity (axb)- (cx d) = (a-c)(b-d) — (a-d)(b-c), 
and 

Bk-Bl = k 2 A\-{k-Akf. (F.75) 
Plugging these pieces into Eq. (F.21) we find that 



Proca 



d 3 k 



k 2 A° k A _k -k° 2 A k - A- k + k 2 A k ■ A^ 



(2tt) 3 

+ (k ■ A k )(k ■ A- k ) + ^ 2 AlA°_ k + V 2 A k ■ A- 



cos2k°t 



+ k o2 At + k 2 A° k 2 
V Af + tfAl 



2k {J A k ] k ■ A k + k z A{ -(k-A 



(F.76) 



First let's prove that the coefficient of the cos term is 0. It's clear that the 
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A k • A^k terms all cancel as k 2 + /z 2 = 
k° 2 A k A°_ k -(k-A k )(k-A- k ) = 



k° . The remaining part is 

(k°A° k -k- A k )(k°A°_ k + k-A- k ) 
+ k°A°_ k k ■ A k - k°A° k k ■ A- k 
(k»A,)(k G A\ + k-A- k ) 
0, (F.77) 



where we note that the pieces added to the end of the first line sum to zero. 
Using the definition of k° to combine terms we have the Ep mca integrand as 

2k° 2 rf k +k o2 A k 2 

-2k°A° k k -A k -(k- A k ) 2 = 2 k° 2 A\ - 2{k ■ A k f 

= 2k« 2 (Al-Al 2 ) 

= —2 k° 2 A^A^, (F.78) 

where we have repeatedly exploited k^A^ = 0, which implies that 

(k^A^f = k° 2 Af- 2k°A° k k -A k +(k-A k y = 0. (F.79) 

With EA=i, 2 ,3 e A e r = -9^ + k^k u /n 2 , k^A^ = 0, and Eq. (F.68) we finally 
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arrive at 



E 



Proca 



2k° 2 \ex-Af 



A=l,2,3 



d 3 k 



(2tt) 3 2 



A=l,2,3 



P 



P 



k ■ p' k ■ p 



d 3 k 



2p ■ p' 



m 



m 



(2tt) 3 2 \(k -p')(k -p) (k-p') 2 (k-pf 



(F.80) 
. (F.81) 



F.2 Field Theory Calculation 
F.2.1 Bremsstrahlung Radiation 

We wish to calculate the field theory analog of the previous classical calcula- 
tions. Following Peskin pgs. 182-183 [251], consider the diagrams in Fig. F.2 
contributing to bremsstrahlung radiation. 




Figure F.2: Diagrams contributing to the probability of emitting a bremsstrahlung 
photon. 

Then, taking Ai to denote the part of the amplitude that comes from the 
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external field "tickling" the charge, we find that 



iM 



ieu(p') ( Mo(p',p — k) 



+ m) 



(p — k) 2 — nv 



7 •£*(*) 



+ k + m) 

+ 7 - e ^ ) (/+fe)"-m^ ° (P+fe ' p)) ' /(/,) - 



(F.82) 



In order to make contact with the classical limit we assume a soft photon: 
| A; | <C \p — p\. Then we can take 



M (p' + k,p) w M (p',p- k) ~M (p',p), 



(F.83) 



and we can ignore the k* in the numerators. 

Using Eq. (E.ll) we have that Eq. (F.82) becomes 



iM. = u(p') [M. (p' ,p)]u(p) 



p' ■ e* p ■ e" 



p 1 ■ k p ■ k 



(F.84) 



Summing over polarization states, integrating over phase space, and taking 
the conditional probability yields 



(iV 7 ) 



d 3 k 1 
(27r) 3 2^' 



E 



p p 



k ■ p' k ■ p 



(F.85) 



For massless QED we have that in a diagram that satisfies the Ward iden- 



tity e^e* = —g^v Then we have that the energy radiated, to lowest order, 
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is 

j-,qft f d 3 k e 2 / 2p ■ p' m 2 m 2 \ 

EM = J (27)3 ~2 \{k-p'){k-p) " {k-p') 2 ~ {k-p) 2 ) ' ( ' 86) 

in exact agreement with the classical E&M result, Eq. (F.59). 

For massive Proca we have the additional longitudinal polarization that 



the photon can radiate into; the polarization sum is modified to e*e v 



A _*A 



— Q^iu + k^k v /fi 2 . Additionally k 2 = m 2 ^ in the propagators. The soft 
photon limit corresponds to the limit of small k°,k, hence m 2 <^2p-k. To the 
same level of approximation, then, we ignore both the photon momentum in 
the numerators and the photon masses in the denominators of the propagators. 
Then we have that 



Hit. 



d 3 k e 2 ( 2p ■ p' m 2 m 2 



l>1 " ,;> ' (2tt) 3 2 \{k-p'){k-p) {k-p') 2 {k-p) 2 

{k-p') 2 {k-p) 2 2{k-p'){k-p) 



{k-p') 2 ^ 2 {k-p) 2 iJ? {k ■ p'){k ■ p)ji 2 

2 ™2 \ 

(F.87) 



d 3 k e 2 ( 2p ■ p' nr in 



(2vr) 3 2 \{k-p'){k-p) {k-p') 2 {k-p) 2 



since the extra terms all cancel, we may use the equation interchangeably for 
either the massive or massless case. 
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F.3 Comparison 

F.3.1 Full Result: Massive Field, Charge, with Recoil 

We wish to make a connection to the zeroth order radiation emitted from the 
production of a hard parton at midrapidity. To that end, let's consider the 
massive quark and field configuration first. Using Ivan's notation: a four- 
vector is not bold-faced (e.g. p, k), a three- vector is bold-faced with a vector 
symbol (e.g. p, k), and a transverse two- vector is bold-faced without a vector 
symbol (e.g. p, k). Minkowski four- vectors are written with parentheses, 
(); light-cone four- vectors with brackets, []. In this notation our momenta of 
interest are 



p = [M Q ,M Q ,0] (F.88) 

P = (1 _ x)jg V k ] (F.89) 

_i_ 2 

k = [xE+, xE+ \ K (F-90) 



In order to compute the Proca result we must convert these into three 
vectors. Doing so faithfully and using these in Eq. (F.80) yields to lowest 
order, when expanding in powers of 1/E + , 



dE 4(l-;r) 2 (k 2 + (l-:r) 2 m 2 ) 



d?k 



C— '—± - - (F91) 

(k 2 + (l-x) 2 m2+x 2 M 2 ) 2 ' 



where C = (2(2tt) 3 ) is a constant that will be have the same value for 
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the rest of this section. There are two interesting features in this result that 
are different from Magda's zeroth order massive QCD result (ignoring the 
additional difference from the error in her published result that she corrected 
in her thesis). The first is the overall factor of (1 — x) 2 that does not reproduce 
Ivan's massless result, which is C /k 2 . While this may be swept under the rug as 
1 — x ~ 1 in the small- x limit, we will actually recover this overall factor below. 
The other difference is the (1 — x) 2 m 2 term in the numerator: this cannot be 
ignored. Nevertheless this filling in of the "dead cone" in the massive field 
case is not so far-fetched; the additional longitudinal polarization could easily 
create a mode to emit radiation collinearly. 

Plugging the above momenta into Eq. (F.87) and expanding in powers of 
1/E + , the lowest order QFT calculation gives 



exactly reproducing the classical limit. This gives us quite a bit of confidence 
in our results. 

F.3.2 Massive Field and Charge, No Recoil 

Interestingly the result when there is no recoil, 



dE 



C 



4(1 -x) 2 (k 2 + (1 -x) 2 m 2 ' 
(k 2 + (1 - x) 2 m 2 + x 2 M 2 ) 



) 



2 ! 



(F.92) 



p' = [E 



<+ 



k 2 + M 2 

E+ 




(F.93) 
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with the other momenta the same, is 

dE = c A((l + x)*V + m*) 
<Pk ((l + a; ) 2 k 2 + m2 + x 2 M 2 ) 2 

F.3.3 Massive Charge and Massless Field with Recoil 

In the case of a massless field the QFT and classical results are necessarily the 
same; we have 

™ = C lii^g (F 95) 

d 3 k (k 2 + x 2 M 2 ) 2 ' 

Everything is consistent with taking m 7 — > from our previous result. Again 

we see the overall factor of (1 — x) 2 out front. We're tempted to simply take 

Mq — > to recover the completely massless limit. However in the massless 

limit we must actually solve a completely different problem: for a massless 

charge we cannot take it to be at rest for t < 0. Perhaps of anecdotal interest 

the totally massless with recoil bremsstrahlung result, using 

p = [E + ,0,0] (F.96) 
P' = [(l-x)E + , {1 _ x) £ + ,-k] (F.97) 

1^2 _i_ m 2 

k = [xE + ,—^,k], (F.98) 



is 



dE x 2 



this is clearly not what we are interested in. 
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What we will do is to solve a problem more relevant to finding the zeroth 
order contribution than bremsstrahlung radiation: beta decay radiation (pair 
production) . 



F.4 Pair Production 



F.4.1 Classical 



We now take as our current a (charge) pair production 



j»(x) = e dr 



d ^5^ ((x - yi (r)) - ((x - y 2 (r)) 



, (F.100) 



where we have 



< 2 (r) = 



0, 

"11,2 



T, 



T < 



T > 



(F.101) 



we have allowed ourselves the freedom to have both difference masses and 
different momenta for the two (opposite) charges. 
The Fourier transform of our current is 



f(k) = J d 4 xe ik - x j»(x) 



= e lim / dr 



mi 



l>'l ( :( ■ 

m 2 



m 2 J 



ie 



Pi 



P2 



k ■ pi + ie k ■ pi + ie 



(F.102) 



The only difference classically from the bremsstrahlung result is that the 
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"k-p = 0" pole is shoved below the real axis. Since we throw away these poles 
anyway we know that we may safely use our previously derived equations for 
the radiated energy lost, Eq. (F.59) and Eq. (F.80). All we need to do is alter 
the p momentum appropriately before plugging in. 

F.4.2 QFT 

We'd like to know whether consideration of pair production, Fig. F.3, results 
in an amplitude different from the usual bremsstrahlung diagrams, Fig. F.2. 
The amplitude from Fig. F.3 is 



iM. ~ — ieu(p')u(p)j • e*(k) 
i{p" + ^ + m) 



{p' + k) 2 — m 2 



(p + k) 2 — m 2 

M u---u, (F.103) 



where we see that the sign difference between the terms no longer comes from 
the switch in sign of the (p — k) 2 vs. the (p' + k) 2 term in the denominator but 
rather from the opposite sign of charge. We may safely employ our previous 
derivation, Eqs. (F.86,F.87), using momenta associated with pair production. 
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Pi Q e 




p, Qe 



ft 



k 



+ 



p', -Qe 



p\ -Qe 



Figure F.3: Diagrams contributing to the probability of emitting a photon in a 
pair production event. 

F.4.3 Comparison 

Massive Field, Charge, with Recoil 

It turns out with 



P = Wl/E 




,0] 

k 2 + ilf2 



(F.104) 



p' = [(l-x)E 



•+ 



(F.105) 



(l-x)E+ 



[xE + , 



k 2 + ml 




(F.106) 



k 



xE+ 



the Proca result, Eq. (F.91), is unchanged. 
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Massless Field and Charge with Recoil 

With 



the E&M result is 



p = [0,E + ,0] (F.107) 
P' = l(l-x)E\ (F.108) 

k = [xE + ,-^,k], (F.109) 



where we now find that the limiting form of the massive result is now recov- 
ered. 



